oft_petsc_la Module Reference

Detailed Description

PETSc vector implementation.

Authors

Chris Hansen

Date

January 2013

Data Types
type	oft_petsc_matrix
	PETSc matrix class. More...
type	oft_petsc_vector
	PETSc vector class. More...

Functions/Subroutines
subroutine	mat_add_values (self, i_inds, j_inds, b, n, m, iblock, jblock, loc_cache)
	Add values to a PETSc matrix.
subroutine	mat_apply_cvec (self, a, b)
	Apply the matrix to a field.
subroutine	mat_apply_vec (self, a, b)
	Apply the matrix to a field.
subroutine	mat_applyt_cvec (self, a, b)
	Apply the matrix to a field.
subroutine	mat_applyt_vec (self, a, b)
	Apply the matrix to a field.
subroutine	mat_assemble (self, diag)
	Finish assembly of matrix.
subroutine	mat_delete (self)
	Delete matrix.
subroutine	mat_set_values (self, i_inds, j_inds, b, n, m, iblock, jblock)
	Set values of a PETSc matrix.
subroutine	mat_zero (self)
	Zero all entries in matrix.
subroutine	mat_zero_rows (self, nrows, irows, iblock, keep_diag)
	Zero all entries in the specified rows.
integer(i4) function	oft_petsc_matrix_cast (self, source)
	Attempt to cast a matrix object to a oft_petsc_matrix.
integer(i4) function	oft_petsc_vector_cast (self, source)
	Cast oft_vector to oft_petsc_la::oft_petsc_vector.
subroutine	vec_add (self, gamma, alpha, a, beta, b)
	Add vectors.
subroutine	vec_delete (self)
	Finalize vector.
function	vec_dot_cvec (self, a)
	Dot product with a second vector.
real(r8) function	vec_dot_vec (self, a)
	Dot product with a second vector.
subroutine	vec_get_local (self, array, iblock)
	Set all elements to a scalar.
subroutine	vec_get_slice (self, array, iblock)
	Set all elements to a scalar.
real(r8) function, dimension(n)	vec_mdot_cvec (self, a, n)
	Dot product with a second vector.
real(r8) function, dimension(n)	vec_mdot_vec (self, a, n)
	Dot product with a second vector.
subroutine	vec_mult (self, a, div_flag)
	Multiply fields element by element.
subroutine	vec_new_cvec (self, new)
	Create a new vector as a bare copy.
subroutine	vec_new_vec (self, new)
	Create a new vector as a bare copy.
real(r8) function	vec_norm (self, itype)
	Norm of a vector.
subroutine	vec_restore_local (self, array, iblock, add, wait)
	Set all elements to a scalar.
subroutine	vec_restore_slice (self, array, iblock, wait)
	Set all elements to a scalar.
subroutine	vec_scale (self, alpha)
	Scale vector by a scalar.
subroutine	vec_set (self, alpha, iblock, random)
	Set all elements to a scalar or random number.
subroutine	vec_stitch (self, up_method)
	Perform global stitching.
real(r8) function	vec_sum (self)
	Sum reduction over a vector.

Function/Subroutine Documentation

mat_add_values()

subroutine mat_add_values ( class(oft_petsc_matrix), intent(inout)  self, integer(i4), dimension(n), intent(in)  i_inds, integer(i4), dimension(m), intent(in)  j_inds, real(r8), dimension(n,m), intent(in)  b, integer(i4), intent(in)  n, integer(i4), intent(in)  m, integer(i4), intent(in), optional  iblock, integer(i4), intent(in), optional  jblock, integer(i4), dimension(n,m), intent(inout), optional  loc_cache  )

private

Add values to a PETSc matrix.

Parameters

[in]	i_inds	Row indices of entries to set [n]
[in]	j_inds	Column indices of entries to set [m]
[in]	b	Values to add [n,m]
[in]	n	Number of rows in local matrix
[in]	m	Number of columns in local matrix
[in]	iblock	Row block (optional)
[in]	jblock	Column block (optional)

mat_apply_cvec()

subroutine mat_apply_cvec ( class(oft_petsc_matrix), intent(inout)  self, class(oft_cvector), intent(inout), target  a, class(oft_cvector), intent(inout)  b  )

private

Apply the matrix to a field.

b = self * a

Parameters

[in]	a	Source field
[out]	b	Result of matrix product

mat_apply_vec()

subroutine mat_apply_vec ( class(oft_petsc_matrix), intent(inout)  self, class(oft_vector), intent(inout), target  a, class(oft_vector), intent(inout)  b  )

private

Apply the matrix to a field.

b = self * a

Parameters

[in]	a	Source field
[out]	b	Result of matrix product

mat_applyt_cvec()

subroutine mat_applyt_cvec ( class(oft_petsc_matrix), intent(inout)  self, class(oft_cvector), intent(inout)  a, class(oft_cvector), intent(inout)  b  )

private

Apply the matrix to a field.

b = self * a

Parameters

[in]	a	Source field
[out]	b	Result of matrix product

mat_applyt_vec()

subroutine mat_applyt_vec ( class(oft_petsc_matrix), intent(inout)  self, class(oft_vector), intent(inout)  a, class(oft_vector), intent(inout)  b  )

private

Apply the matrix to a field.

b = self * a

Parameters

[in]	a	Source field
[out]	b	Result of matrix product

mat_assemble()

subroutine mat_assemble ( class(oft_petsc_matrix), intent(inout)  self, class(oft_vector), intent(inout), optional, target  diag  )

private

Finish assembly of matrix.

Warning

Currently diagonal extraction is not supported for PETSc matrices

Parameters

[in,out]

diag

Diagonal entries of matrix [nr] (optional)

mat_delete()

subroutine mat_delete ( class(oft_petsc_matrix), intent(inout)  self )

private

Delete matrix.

mat_set_values()

subroutine mat_set_values ( class(oft_petsc_matrix), intent(inout)  self, integer(i4), dimension(n), intent(in)  i_inds, integer(i4), dimension(m), intent(in)  j_inds, real(r8), dimension(n,m), intent(in)  b, integer(i4), intent(in)  n, integer(i4), intent(in)  m, integer(i4), intent(in), optional  iblock, integer(i4), intent(in), optional  jblock  )

private

Set values of a PETSc matrix.

Parameters

[in]	i_inds	Row indices of entries to set [n]
[in]	j_inds	Column indices of entries to set [m]
[in]	b	Values to set [n,m]
[in]	n	Number of rows in local matrix
[in]	m	Number of columns in local matrix
[in]	iblock	Row block (optional)
[in]	jblock	Column block (optional)

mat_zero()

subroutine mat_zero ( class(oft_petsc_matrix), intent(inout)  self )

private

Zero all entries in matrix.

mat_zero_rows()

subroutine mat_zero_rows ( class(oft_petsc_matrix), intent(inout)  self, integer(i4), intent(in)  nrows, integer(i4), dimension(nrows), intent(in)  irows, integer(i4), intent(in), optional  iblock, logical, intent(in), optional  keep_diag  )

private

Zero all entries in the specified rows.

Parameters

[in]	nrows	Number of rows to zero
[in]	irows	Indices of rows to zero [nrows]
[in]	iblock	Row block (optional)
[in]	keep_diag	Keep diagonal entries

oft_petsc_matrix_cast()

integer(i4) function oft_petsc_matrix_cast ( class(oft_petsc_matrix), intent(out), pointer  self, class(oft_matrix), intent(in), target  source  )

private

Attempt to cast a matrix object to a oft_petsc_matrix.

The source matrix must be oft_petsc_matrix or a child class, otherwise an error will be thrown

Parameters

[out]	self	Pointer to cast oft_petsc_matrix
[in]	source	Source matrix to cast

oft_petsc_vector_cast()

integer(i4) function oft_petsc_vector_cast	(	class(oft_petsc_vector), intent(out), pointer	self,
		class(oft_vector), intent(in), target	source
	)

Cast oft_vector to oft_petsc_la::oft_petsc_vector.

Parameters

[out]	self	Object of desired type, unassociated if cast fails
[in]	source	Source object to cast

Returns

Error flag

vec_add()

subroutine vec_add ( class(oft_petsc_vector), intent(inout)  self, real(r8), intent(in)  gamma, real(r8), intent(in)  alpha, class(oft_vector), intent(in), target  a, real(r8), intent(in), optional  beta, class(oft_vector), intent(in), optional, target  b  )

private

Add vectors.

self = \( \gamma \) self + \( \alpha \) a + \( \beta \) b

Parameters

[in]	gamma	Scale of source vector
[in]	alpha	Scale of first vector
[in]	a	First vector to add
[in]	beta	Scale of second vector (optional)
[in]	b	Second vector to add (optional)

vec_delete()

subroutine vec_delete ( class(oft_petsc_vector), intent(inout)  self )

private

Finalize vector.

vec_dot_cvec()

function vec_dot_cvec ( class(oft_petsc_vector), intent(inout)  self, class(oft_cvector), intent(in), target  a  )

private

Dot product with a second vector.

Parameters

[in]

a

Second vector for dot product

Returns

\( \sum_i self_i a_i \)

vec_dot_vec()

real(r8) function vec_dot_vec ( class(oft_petsc_vector), intent(inout)  self, class(oft_vector), intent(in), target  a  )

private

Dot product with a second vector.

Parameters

[in]

a

Second vector for dot product

Returns

\( \sum_i self_i a_i \)

vec_get_local()

subroutine vec_get_local ( class(oft_petsc_vector), intent(inout)  self, real(r8), dimension(:), intent(inout), pointer  array, integer(i4), intent(in), optional  iblock  )

private

Set all elements to a scalar.

Parameters

[in]

alpha

Updated field value

Note

This subroutine is a dummy routine used to specify the interface of the member function and catch errors in uninitialized fields.

vec_get_slice()

subroutine vec_get_slice ( class(oft_petsc_vector), intent(inout)  self, real(r8), dimension(:), intent(inout), pointer  array, integer(i4), intent(in), optional  iblock  )

private

Set all elements to a scalar.

Parameters

[in]

alpha

Updated field value

Note

This subroutine is a dummy routine used to specify the interface of the member function and catch errors in uninitialized fields.

vec_mdot_cvec()

real(r8) function, dimension(n) vec_mdot_cvec ( class(oft_petsc_vector), intent(inout)  self, type(oft_vector_ptr), dimension(n), intent(in)  a, integer(i4), intent(in)  n  )

private

Dot product with a second vector.

Parameters

[in]

a

Second vector for dot product

Returns

\( \sum_i self_i a_i \)

vec_mdot_vec()

real(r8) function, dimension(n) vec_mdot_vec ( class(oft_petsc_vector), intent(inout)  self, type(oft_vector_ptr), dimension(n), intent(in)  a, integer(i4), intent(in)  n  )

private

Dot product with a second vector.

Parameters

[in]

a

Second vector for dot product

Returns

\( \sum_i self_i a_i \)

vec_mult()

subroutine vec_mult ( class(oft_petsc_vector), intent(inout)  self, class(oft_vector), intent(in), target  a, logical, intent(in), optional  div_flag  )

private

Multiply fields element by element.

\( self_i = self_i * a_i \)

Parameters

[in]

a

First field to add

vec_new_cvec()

subroutine vec_new_cvec ( class(oft_petsc_vector), intent(in)  self, class(oft_cvector), intent(out), pointer  new  )

private

Create a new vector as a bare copy.

Creates a new vector of the same size and with the same gobal mapping. The new vector is also initialized to 0.

Parameters

[out]

new

New vector

vec_new_vec()

subroutine vec_new_vec ( class(oft_petsc_vector), intent(in)  self, class(oft_vector), intent(out), pointer  new  )

private

Create a new vector as a bare copy.

Creates a new vector of the same size and with the same gobal mapping. The new vector is also initialized to 0.

Parameters

[out]

new

New vector

vec_norm()

real(r8) function vec_norm ( class(oft_petsc_vector), intent(inout)  self, integer(i4), intent(in)  itype  )

private

Norm of a vector.

Parameters

[in]

itype

Type of norm (1-> 1-norm, 2-> 2-norm, 3-> Inf-norm)

Returns

norm(self)

vec_restore_local()

subroutine vec_restore_local ( class(oft_petsc_vector), intent(inout)  self, real(r8), dimension(:), intent(inout)  array, integer(i4), intent(in), optional  iblock, logical, intent(in), optional  add, logical, intent(in), optional  wait  )

private

Set all elements to a scalar.

Parameters

[in]

alpha

Updated field value

Note

This subroutine is a dummy routine used to specify the interface of the member function and catch errors in uninitialized fields.

vec_restore_slice()

subroutine vec_restore_slice ( class(oft_petsc_vector), intent(inout)  self, real(r8), dimension(:), intent(in)  array, integer(i4), intent(in), optional  iblock, logical, intent(in), optional  wait  )

private

Set all elements to a scalar.

Parameters

[in]

alpha

Updated field value

Note

This subroutine is a dummy routine used to specify the interface of the member function and catch errors in uninitialized fields.

vec_scale()

subroutine vec_scale ( class(oft_petsc_vector), intent(inout)  self, real(r8), intent(in)  alpha  )

private

Scale vector by a scalar.

self = \( \alpha \) self

Parameters

[in]

alpha

Factor to scale field

vec_set()

subroutine vec_set ( class(oft_petsc_vector), intent(inout)  self, real(r8), intent(in)  alpha, integer(i4), intent(in), optional  iblock, logical, intent(in), optional  random  )

private

Set all elements to a scalar or random number.

Parameters

[in]	alpha	Updated vector value
[in]	random	Set to random number, if true alpha is ignored (optional)

vec_stitch()

subroutine vec_stitch ( class(oft_petsc_vector), intent(inout)  self, integer(i4), intent(in)  up_method  )

private

Perform global stitching.

Parameters

[in]

up_method

Type of stitching to perform

vec_sum()

real(r8) function vec_sum ( class(oft_petsc_vector), intent(inout)  self )

private

Sum reduction over a vector.

Returns

\( \sum_i self_i \)