linalg.block

Classes

BlockLinearOperator(V1, V2[, blocks])	Linear operator that can be written as blocks of other Linear Operators.
BlockVector(V[, blocks])	Block of Vectors, which is an element of a BlockVectorSpace.
BlockVectorSpace(*spaces[, connectivity])	Product Vector Space V of two Vector Spaces (V1,V2) or more.

Details

class BlockVectorSpace(*spaces, connectivity=None)

Bases: VectorSpace

Product Vector Space V of two Vector Spaces (V1,V2) or more.

Parameters:

*spacespsydac.linalg.basic.VectorSpace

A list of Vector Spaces.

property dimension

The dimension of a product space V = (V1, V2, …] is the cardinality (i.e. the number of vectors) of a basis of V over its base field.

property dtype

The data type of the field over which the space is built.

zeros()

Get a copy of the null element of the product space V = [V1, V2, …]

Returns:

nullBlockVector

A new vector object with all components equal to zero.

axpy(a, x, y)

Increment the vector y with the a-scaled vector x, i.e. y = a * x + y, provided that x and y belong to the same vector space V (self). The scalar value a may be real or complex, depending on the field of V.

Parameters:

ascalar

The scaling coefficient needed for the operation.

xBlockVector

The vector which is not modified by this function.

yBlockVector

The vector modified by this function (incremented by a * x).

property spaces

property parallel

Returns True if the memory is distributed.

property starts

property ends

property pads

property n_blocks

property connectivity

class BlockVector(V, blocks=None)

Bases: Vector

Block of Vectors, which is an element of a BlockVectorSpace.

Parameters:

Vpsydac.linalg.block.BlockVectorSpace

Space to which the new vector belongs.

blockslist or tuple (psydac.linalg.basic.Vector)

List of Vector objects, belonging to the correct spaces (optional).

property space

Vector space to which this vector belongs.

property dtype

The data type of the vector field V this vector belongs to.

dot(v)

Evaluate the scalar product with the vector v of the same space.

Parameters:

vVector

Vector belonging to the same space as self.

copy(out=None)

Ensure x.copy(out=x) returns x and not a new object.

conjugate(out=None)

Compute the complex conjugate vector.

If the field is real (i.e. self.dtype in (np.float32, np.float64)) this method is equivalent to copy. If the field is complex (i.e. self.dtype in (np.complex64, np.complex128)) this method returns the complex conjugate of self, element-wise.

The behavior of this function is similar to numpy.conjugate(self, out=None).

property ghost_regions_in_sync

update_ghost_regions()

start_update_interface_ghost_regions()

end_update_interface_ghost_regions(req)

exchange_assembly_data()

property n_blocks

property blocks

toarray(order='C')

Convert to Numpy 1D array.

toarray_local(order='C')

Convert to petsc Nest vector.

topetsc()

Convert to petsc data structure.

class BlockLinearOperator(V1, V2, blocks=None)

Bases: LinearOperator

Linear operator that can be written as blocks of other Linear Operators. Either the domain or the codomain of this operator, or both, should be of class BlockVectorSpace.

Parameters:

V1psydac.linalg.block.VectorSpace

Domain of the new linear operator.

V2psydac.linalg.block.VectorSpace

Codomain of the new linear operator.

blocksdict | (list of lists) | (tuple of tuples)

LinearOperator objects (optional).

1. ‘blocks’ can be dictionary with

   . key = tuple (i, j), where i and j are two integers >= 0 . value = corresponding LinearOperator Lij

2. ‘blocks’ can be list of lists (or tuple of tuples) where blocks[i][j]

   is the LinearOperator Lij (if None, we assume null operator)

property domain

The domain of the linear operator - an element of Vectorspace

property codomain

The codomain of the linear operator - an element of Vectorspace

property dtype

conjugate(out=None)

conj(out=None)

tosparse(**kwargs)

Convert to any Scipy sparse matrix format.

toarray(**kwargs)

Convert to Numpy 2D array.

dot(v, out=None)

Apply linear operator to Vector v. Result is written to Vector out, if provided.

transpose(conjugate=False, out=None)

” Return the transposed BlockLinearOperator, or the Hermitian Transpose if conjugate==True

Parameters:

conjugateBool(optional)

True to get the Hermitian adjoint.

outBlockLinearOperator(optional)

Optional out for the transpose to avoid temporaries

diagonal(*, inverse=False, sqrt=False, out=None)

Get the coefficients on the main diagonal as another BlockLinearOperator object.

Parameters:

inversebool

If True, get the inverse of the diagonal. (Default: False). Can be combined with sqrt to get the inverse square root.

sqrtbool

If True, get the square root of the diagonal. (Default: False). Can be combined with inverse to get the inverse square root.

outBlockLinearOperator

If provided, write the diagonal entries into this matrix. (Default: None).

Returns:

BlockLinearOperator

The matrix which contains the main diagonal of self (or its inverse).

property blocks

Immutable 2D view (tuple of tuples) of the linear operator, including the empty blocks as ‘None’ objects.

property n_block_rows

property n_block_cols

property nonzero_block_indices

Tuple of (i, j) pairs which identify the non-zero blocks: i is the row index, j is the column index.

update_ghost_regions()

exchange_assembly_data()

remove_spurious_entries()

property ghost_regions_in_sync

transform(operation)

Applies an operation on each block in this BlockLinearOperator.

Parameters:

operationLinearOperator -> LinearOperator

The operation which transforms each block.

backend()

copy(out=None)

Create a copy of self, that can potentially be stored in a given BlockLinearOperator.

Parameters:

outBlockLinearOperator(optional)

The existing BlockLinearOperator in which we want to copy self.

Returns:

BlockLinearOperator

The copy of self, either stored in the given BlockLinearOperator out (if provided) or in a new one. In the corner case where out=self the self object is immediately returned.

topetsc()

Convert to petsc data structure.

compute_interface_matrices_transpose()

set_backend(backend)