This routine builds up the matrix for least square fit used in linear interpolation.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(tem_matrix_type), | intent(out) | :: | me | Matrix to fill |
||
| integer, | intent(in) | :: | 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 |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| 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 |