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) | :: | 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 |

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 |