![]() |
programmer's documentation
|
This function builds the matrix of advection/diffusion for a vector field with a tensorial diffusivity. More...
Functions/Subroutines | |
subroutine | matrvv (ncelet, ncel, nfac, nfabor, iconvp, idiffp, ndircp, isym, nfecra, thetap, ifacel, ifabor, coefbu, cofbfu, fimp, flumas, flumab, viscf, viscb, da, xa) |
This function builds the matrix of advection/diffusion for a vector field with a tensorial diffusivity.
The advection is upwind, the diffusion is not reconstructed. The matrix is splitted into a diagonal block (3x3 times number of cells) and an extra diagonal part (of dimension 2 times 3x3 the number of internal faces).
subroutine matrvv | ( | integer | ncelet, |
integer | ncel, | ||
integer | nfac, | ||
integer | nfabor, | ||
integer | iconvp, | ||
integer | idiffp, | ||
integer | ndircp, | ||
integer | isym, | ||
integer | nfecra, | ||
double precision | thetap, | ||
integer, dimension(2,nfac) | ifacel, | ||
integer, dimension(nfabor) | ifabor, | ||
double precision, dimension(3,3,nfabor) | coefbu, | ||
double precision, dimension(3,3,nfabor) | cofbfu, | ||
double precision, dimension(3,3,ncelet) | fimp, | ||
double precision, dimension(nfac) | flumas, | ||
double precision, dimension(nfabor) | flumab, | ||
double precision, dimension(3,3,nfac) | viscf, | ||
double precision, dimension(nfabor) | viscb, | ||
double precision, dimension(3,3,ncelet) | da, | ||
double precision, dimension(3,3,isym,nfac) | xa | ||
) |
[in] | ncelet | number of extended (real + ghost) cells |
[in] | ncel | number of cells |
[in] | nfac | number of interior faces |
[in] | nfabor | number of boundary faces |
[in] | iconvp | indicator
|
[in] | idiffp | indicator
|
[in] | ndircp | indicator
|
[in] | isym | indicator
|
[in] | thetap | weightening coefficient for the theta-schema,
|
[in] | ifacel | cell indexes of interior faces |
[in] | ifabor | no de l'elt voisin d'une face de bord |
[in] | coefbu | boundary condition array for the variable (Impplicit part - 3x3 tensor array) |
[in] | cofbfu | boundary condition array for the variable flux (Impplicit part - 3x3 tensor array) |
[in] | flumas | mass flux at interior faces |
[in] | flumab | mass flux at border faces |
[in] | viscf | ![]() |
[in] | viscb | ![]() |
[out] | da | diagonal part of the matrix |
[out] | xa | extra interleaved diagonal part of the matrix |