solve_axb Subroutine

public subroutine solve_axb(a, b, d, ierr, yes_verbose)

solves AxB=D for D

Arguments

Type IntentOptional Attributes Name
real(kind=long), intent(in) :: a(:,:)

matrix A of shape NxN

real(kind=long), intent(in) :: b(:)

vector B of len N

real(kind=long), intent(out) :: d(size(b,1))

result vector of len N

integer, intent(out) :: ierr

output status -1 is solution could not be found (A is singular)

logical, intent(in), optional :: yes_verbose

print message if A is singular