build_matrixLSF_linearIntp Subroutine

private subroutine build_matrixLSF_linearIntp(me, QQ, nDims, nSources, cxDirRK, neighDir, nCoeffs, success)

This routine builds up the matrix for least square fit used in linear interpolation.

Arguments

TypeIntentOptionalAttributesName
type(tem_matrix_type), intent(out) :: me

Matrix to fill

integer, intent(in) :: QQ

Number of stencil directions

integer, intent(in) :: nDims

Number of dimensions

integer, intent(in) :: nSources

Number of sources from coarser found

real(kind=rk), intent(in) :: cxDirRK(3,QQ)

Stencil directions

integer, intent(in) :: neighDir(nSources)

direction in which sources are found

integer, intent(in) :: nCoeffs

nUnknown coeffs

logical, intent(out) :: success

success if false if matrix is singular reduce interpolation order


Calls

proc~~build_matrixlsf_linearintp~~CallsGraph proc~build_matrixlsf_linearintp build_matrixLSF_linearIntp proc~alloc_matrix alloc_matrix proc~build_matrixlsf_linearintp->proc~alloc_matrix proc~polylinear_1d polyLinear_1D proc~build_matrixlsf_linearintp->proc~polylinear_1d proc~polylinear_3d polyLinear_3D proc~build_matrixlsf_linearintp->proc~polylinear_3d proc~tem_matrix_dump tem_matrix_dump proc~build_matrixlsf_linearintp->proc~tem_matrix_dump proc~tem_abort tem_abort proc~build_matrixlsf_linearintp->proc~tem_abort proc~invert_matrix invert_matrix proc~build_matrixlsf_linearintp->proc~invert_matrix proc~polylinear_2d polyLinear_2D proc~build_matrixlsf_linearintp->proc~polylinear_2d proc~alloc_matrix->proc~tem_abort mpi_abort mpi_abort proc~tem_abort->mpi_abort proc~invert_matrix->proc~tem_abort dgetrf dgetrf proc~invert_matrix->dgetrf dgetri dgetri proc~invert_matrix->dgetri

Called by

proc~~build_matrixlsf_linearintp~~CalledByGraph proc~build_matrixlsf_linearintp build_matrixLSF_linearIntp proc~append_intpmatrixlsf append_intpMatrixLSF proc~append_intpmatrixlsf->proc~build_matrixlsf_linearintp interface~append~38 append interface~append~38->proc~append_intpmatrixlsf

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iDir
integer, private :: iSrc
type(tem_matrix_type), private :: tmp_matrix

Each row represents a polynomial evaluated at coord of elements in

real(kind=rk), private :: inv_AtA(nCoeffs,nCoeffs)
real(kind=rk), private :: AtA(nCoeffs,nCoeffs)
integer, private :: errCode