matrix_sum Function

public function matrix_sum(aij, akj, undeflag)

Calculate product for a component in a matrix

Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: aij(:,:)

Matrix A1

double precision, intent(in) :: akj(size(aij,1),size(aij,2))

Matrix A2

logical, intent(in), optional :: undeflag(size(aij,1),size(aij,2))

Undefined flag at each sampling point

Return Value doubleprecision