This routine builds up the matrix for least square fit used in linear and quadratic interpolation.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(tem_intpMatrixLSF_type), | intent(inout) | :: | me | intpMatrix for LSF fill |
||
integer, | intent(inout) | :: | order | interpolation order calculated for current element depending on nSources if quadratic LSF matrix is singular fall back to linear |
||
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(out) | :: | pos | Pointer to position of interpolation matrix in growing array of matrix |
||
logical, | intent(out) | :: | success | success if false if matrix is singular reduce interpolation order |