petsc-3.9.4 2018-09-11
Report Typos and Errors

PetscFEIntegrateResidual

Produce the element residual vector for a chunk of elements by quadrature integration

Synopsis

#include "petscfe.h" 
PetscErrorCode PetscFEIntegrateResidual(PetscFE fem, PetscDS prob, PetscInt field, PetscInt Ne, PetscFECellGeom *cgeom,
                                        const PetscScalar coefficients[], const PetscScalar coefficients_t[], PetscDS probAux, const PetscScalar coefficientsAux[], PetscReal t, PetscScalar elemVec[])
Not collective

Input Parameters

fem - The PetscFE object for the field being integrated
prob - The PetscDS specifying the discretizations and continuum functions
field - The field being integrated
Ne - The number of elements in the chunk
cgeom - The cell geometry for each cell in the chunk
coefficients - The array of FEM basis coefficients for the elements
coefficients_t - The array of FEM basis time derivative coefficients for the elements
probAux - The PetscDS specifying the auxiliary discretizations
coefficientsAux - The array of FEM auxiliary basis coefficients for the elements
t - The time

Output Parameter

elemVec -the element residual vectors from each element

Note

Loop over batch of elements (e):
  Loop over quadrature points (q):
    Make u_q and gradU_q (loops over fields,Nb,Ncomp) and x_q
    Call f_0 and f_1
  Loop over element vector entries (f,fc --> i):
    elemVec[i] += \psi^{fc}_f(q) f0_{fc}(u, \nabla u) + \nabla\psi^{fc}_f(q) \cdot f1_{fc,df}(u, \nabla u)

See Also

PetscFEIntegrateResidual()

Level

developer

Location

src/dm/dt/interface/dtfe.c
Index of all DT routines
Table of Contents for all manual pages
Index of all manual pages