6.8 psb_spins — Insert a set of coefficients into a sparse matrix

call psb_spins(nz, ia, ja, val, a, desc_a, info [,local])
call psb_spins(nr, irw, irp, ja, val, a, desc_a, info [,local])

Type:

Asynchronous.

On Entry

nz

the number of coefficients to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer scalar.

nr

the number of rows to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer scalar.

irw

the first row to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer scalar.

ia

the row indices of the coefficients to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer array of size nz.

irp

the row pointers of the coefficients to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer array of size nr + 1.

ja

the column indices of the coefficients to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an integer array of size nz.

val

the coefficients to be inserted.
Scope:local.
Type:required.
Intent: in.
Specified as: an array of size nz. Must be of the same type and kind of the coefficients of the sparse matrix a.

desc_a

The communication descriptor.
Scope: local.
Type: required.
Intent: inout.
Specified as: a variable of type psb_desc_type.

local

Whether the entries in the indices vectors ia, ja are already in local numbering.
Scope:local.
Type:optional.
Specified as: a logical value; default: .false..

On Return

a

the matrix into which coefficients will be inserted.
Scope:local
Type:required
Intent: inout.
Specified as: a structured data of type psb_Tspmat_type.

desc_a

The communication descriptor.
Scope: local.
Type: required.
Intent: inout.
Specified as: a variable of type psb_desc_type.

info

Error code.
Scope: local
Type: required
Intent: out.
An integer value; 0 means no error has been detected.

Notes

1. On entry to this routine the descriptor may be in either the build or assembled state.
2. On entry to this routine the sparse matrix may be in either the build or update state.
3. If the descriptor is in the build state, then the sparse matrix must also be in the build state; the action of the routine is to (implicitly) call psb_cdins to add entries to the sparsity pattern; each sparse matrix entry implicitly defines a graph edge, that is passed to the descriptor routine for the appropriate processing;
4. The input data can be passed in either COO or CSR formats;
5. In COO format the coefficients to be inserted are represented by the ordered triples ia(i),ja(i),val(i), for i = 1,…,nz; these triples are arbitrary;
6. In CSR format the coefficients to be inserted for each input row i = 1,nr are represented by the ordered triples (i + irw - 1),ja(j),val(j), for j = irp(i),…,irp(i + 1) - 1; these triples should belong to the current process, i.e. i+irw-1 should be one of the local indices, but are otherwise arbitrary;
7. There is no requirement that a given row must be passed in its entirety to a single call to this routine: the buildup of a row may be split into as many calls as desired (even in the CSR format);
8. Coefficients from different rows may also be mixed up freely in a single call, according to the application needs;
9. Coefficients from matrix rows not owned by the calling process are treated according to the value of bldmode specified at allocation time; if bldmode was chosen as psb_matbld_remote_ the library will keep track of them, otherwise they are silently ignored;
10. If the descriptor is in the assembled state, then any entries in the sparse matrix that would generate additional communication requirements are ignored;
11. If the matrix is in the update state, any entries in positions that were not present in the original matrix are ignored.