#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <strings.h>
#include <assert.h>
#include "postprocessing.h"

Include dependency graph for postprocessing.c:

Functions

void	invertLocal (int nSpace, double A[3][3], double AI[3][3])

void	postProcessRT0velocityFromP1nc (int nElements_global, int nQuadraturePoints_element, int nDOF_test_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nFreeDOF_element, int freeLocal_element, double detJ, double sqrt_det_g, double n, double elementBarycenters, double quad_a, double quad_f, double w_dV_r, double w_dV_m, double u, double gradu, double a, double f, double r, double mt, double *rt0vdofs)

void	postProcessRT0velocityFromP1nc_sd (int nElements_global, int nQuadraturePoints_element, int nDOF_test_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int nFreeDOF_element, int freeLocal_element, double detJ, double sqrt_det_g, double n, double elementBarycenters, double quad_a, double quad_f, double w_dV_r, double w_dV_m, double u, double gradu, double a, double f, double r, double mt, double *rt0vdofs)

void	postProcessRT0velocityFromP1ncNoMass (int nElements_global, int nQuadraturePoints_element, int nDOF_test_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nFreeDOF_element, int freeLocal_element, double detJ, double sqrt_det_g, double n, double elementBarycenters, double quad_a, double quad_f, double w_dV_r, double u, double gradu, double a, double f, double r, double *rt0vdofs)

void	postProcessRT0velocityFromP1ncNoMass_sd (int nElements_global, int nQuadraturePoints_element, int nDOF_test_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int nFreeDOF_element, int freeLocal_element, double detJ, double sqrt_det_g, double n, double elementBarycenters, double quad_a, double quad_f, double w_dV_r, double u, double gradu, double a, double f, double r, double *rt0vdofs)

void	updateRT0velocityWithAveragedPotentialP1nc (int nElements_global, int nQuadraturePoints_element, int nSpace, double detJ, double quad_a, double phi, double gradphi, double a, double rt0vdofs)

void	updateRT0velocityWithAveragedPotentialP1nc_sd (int nElements_global, int nQuadraturePoints_element, int nSpace, int rowptr, int colind, double detJ, double quad_a, double phi, double gradphi, double a, double rt0vdofs)

void	getElementRT0velocityValues (int nElements_global, int nPoints_element, int nSpace, double x_element, double rt0vdofs_element, double *v_element)

void	getElementBoundaryRT0velocityValues (int nElements_global, int nElementBoundaries_element, int nPoints_elementBoundary, int nSpace, double x_elementBoundary, double rt0vdofs_element, double *v_elementBoundary)

void	getGlobalElementBoundaryRT0velocityValues (int nElementBoundaries_global, int nPoints_elementBoundary, int nSpace, int elementBoundaryElementsArray, double x_elementBoundary_global, double rt0vdofs_element, double v_elementBoundary_global)

void	getGlobalExteriorElementBoundaryRT0velocityValues (int nExteriorElementBoundaries_global, int nPoints_elementBoundary, int nSpace, int elementBoundaryElementsArray, int exteriorElementBoundariesArray, double x_elementBoundary_global, double rt0vdofs_element, double *v_elementBoundary_global)

void	postProcessRT0potentialFromP1nc (int nElements_global, int nQuadraturePoints_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, double uQuadratureWeights_element, double elementBarycenters, double aElementQuadratureWeights, double detJ, double uQuadratureWeights_elementBoundary, double x, double u, double gradu, double x_elementBoundary, double u_elementBoundary, double n, double a, double f, double r, double rt0vdofs, double rt0potential)

void	postProcessRT0potentialFromP1nc_sd (int nElements_global, int nQuadraturePoints_element, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, double uQuadratureWeights_element, double elementBarycenters, double aElementQuadratureWeights, double detJ, double uQuadratureWeights_elementBoundary, double x, double u, double gradu, double x_elementBoundary, double u_elementBoundary, double n, double a, double f, double r, double rt0vdofs, double rt0potential)

void	projectElementBoundaryVelocityToRT0fluxRep (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, double elementBoundaryQuadratureWeights, double n, double v_elementBoundary, double rt0vdofs_element)

void	projectElementBoundaryFluxToRT0fluxRep (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_RT0V_element, int elementBoundaryElementsArray, int elementBoundariesArray, double elementBoundaryQuadratureWeights, double flux_elementBoundary, double *rt0vdofs_element)

void	getElementRT0velocityValuesFluxRep (int nElements_global, int nElementBoundaries_element, int nPoints_element, int nSpace, int nDetVals_element, double nodeArray, int elementNodesArray, double abs_det_J, double x_element, double rt0vdofs_element, double v_element)

void	getElementBoundaryRT0velocityValuesFluxRep (int nElements_global, int nElementBoundaries_element, int nPoints_elementBoundary, int nSpace, int nDetVals_element, double nodeArray, int elementNodesArray, double abs_det_J, double x_elementBoundary, double rt0vdofs_element, double v_elementBoundary)

void	getGlobalElementBoundaryRT0velocityValuesFluxRep (int nElementBoundaries_global, int nPoints_elementBoundary_global, int nSpace, int nDetVals_element, double nodeArray, int elementNodesArray, int elementBoundaryElementsArray, double abs_det_J, double x_elementBoundary_global, double rt0vdofs_element, double *v_elementBoundary_global)

void	getGlobalExteriorElementBoundaryRT0velocityValuesFluxRep (int nExteriorElementBoundaries_global, int nPoints_elementBoundary_global, int nSpace, int nDetVals_element, double nodeArray, int elementNodesArray, int elementBoundaryElementsArray, int exteriorElementBoundariesArray, double abs_det_J, double x_ebqe, double rt0vdofs_element, double v_ebqe)

void	getRT0velocityValuesFluxRep_arbitraryElementMembership (int nElements_global, int nElementBoundaries_element, int nPoints, int nSpace, int nDetVals_element, const double nodeArray, const int elementNodesArray, const double abs_det_J, const double x, const int element_locations, const double rt0vdofs_element, double *v_element)

void	buildLocalBDM1projectionMatrices_orig (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nDOFs_test_element, int nDOFs_trial_element, int nVDOFs_element, double w_dS_f, double ebq_n, double ebq_v, double BDMprojectionMat_element)

void	buildLocalBDM1projectionMatrices (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nDOFs_test_element, int nDOFs_trial_element, int nVDOFs_element, double w_dS_f, double ebq_n, double ebq_v, double BDMprojectionMat_element)

void	buildLocalBDM2projectionMatrices (int degree, int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_elementInterior, int nSpace, int nDOFs_test_element, int nDOFs_trial_boundary_element, int nDOFs_trial_interior_element, int nVDOFs_element, int edgeFlags, double w_dS_f, double ebq_n, double ebq_v, double BDMprojectionMat_element, double q_basis_vals, double w_int_test_grads, double w_int_div_free, double *piola_trial_fun)

void	factorLocalBDM1projectionMatrices (int nElements_global, int nVDOFs_element, double BDMprojectionMat_element, int BDMprojectionMatPivots_element)

void	factorLocalBDM2projectionMatrices (int nElements_global, int nVDOFs_element, double BDMprojectionMat_element, int BDMprojectionMatPivots_element)

void	solveLocalBDM1projection (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nDOFs_test_element, int nVDOFs_element, double BDMprojectionMatFact_element, int BDMprojectionMatPivots_element, double w_dS_f, double ebq_n, double ebq_velocity, double p1_velocity_dofs)

void	buildBDM2rhs (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_elementInterior, int nSpace, int nDOFs_test_element, int nVDOFs_element, int nDOFs_trial_interior_element, double BDMprojectionMatFact_element, int BDMprojectionMatPivots_element, int edgeFlags, double w_dS_f, double ebq_n, double w_interior_grads, double w_interior_divfree, double ebq_velocity, double q_velocity, double p1_velocity_dofs)

void	solveLocalBDM2projection (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int nDOFs_test_element, int nVDOFs_element, double BDMprojectionMatFact_element, int BDMprojectionMatPivots_element, double w_dS_f, double ebq_n, double w_interior_gradients, double q_velocity, double ebq_velocity, double p1_velocity_dofs)

void	solveLocalBDM1projectionFromFlux (int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOFs_test_element, int nVDOFs_element, double BDMprojectionMatFact_element, int BDMprojectionMatPivots_element, int elementBoundaryElementsArray, int elementBoundariesArray, double w_dS_f, double ebq_global_flux, double *p1_velocity_dofs)

void	getElementBDM1velocityValuesLagrangeRep_orig (int nElements_global, int nQuadraturePoints_element, int nSpace, int nDOF_trial_element, int nVDOF_element, double q_v, double p1_velocity_dofs, double *q_velocity)

void	getElementBDM1velocityValuesLagrangeRep (int nElements_global, int nQuadraturePoints_element, int nSpace, int nDOF_trial_element, int nVDOF_element, double q_v, double p1_velocity_dofs, double *q_velocity)

void	getElementBDM2velocityValuesLagrangeRep (int nElements_global, int nQuadraturePoints_element, int nSpace, int nDOF_trial_element, int nVDOF_element, double q_v, double p1_velocity_dofs, double *q_velocity)

void	getElementLDGvelocityValuesLagrangeRep (int nElements_global, int nQuadraturePoints_element, int nSpace, int nDOF_trial_element, int nVDOF_element, double q_v, double velocity_dofs, double *q_velocity)

void	getGlobalExteriorElementBoundaryBDM1velocityValuesLagrangeRep (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int nVDOF_element, int elementBoundaryElementsArray, int exteriorElementBoundariesArray, double ebqe_v, double p1_velocity_dofs, double *ebqe_velocity)

void	getGlobalElementBoundaryBDM1velocityValuesLagrangeRep (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int nVDOF_element, int elementBoundaryElementsArray, int exteriorElementBoundariesArray, double ebqe_v, double p1_velocity_dofs, double *ebq_global_velocity)

void	getElementBoundaryBDM1velocityValuesLagrangeRep (int nElements_global, int nBoundaries_Element, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int nVDOF_element, int elementBoundaryElementsArray, int exteriorElementBoundariesArray, double ebq_v, double p1_velocity_dofs, double *ebq_velocity)

void	calculateConservationResidualGlobalBoundaries (int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int exteriorElementBoundariesToSkip, double dS, double normal, double elementResidual, double velocity, double conservationResidual)

void	sunWheelerGSsweep (int nElements_global, int nElementBoundaries_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double dS, double normal, double sqrt_det_g, double alpha, double fluxCorrection, double conservationResidual)

void	fluxCorrectionVelocityUpdate (int nElements_global, int nElementBoundaries_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double dS, double normal, double fluxCorrection, double vConservative, double *vConservative_element)

void	computeFluxCorrectionPWC (int nElementBoundaries_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, double pwcW, double pwcV, double fluxCorrection)

int	nodeStar_init (int nElements_global, int nNodes_element, int nNodes_global, int nElements_node, int nodeStarElementsArray, int nodeStarElementNeighborsArray, int N_p, int subdomain_dim_p, double subdomain_L_p, double subdomain_R_p, double subdomain_U_p, PROTEUS_LAPACK_INTEGER subdomain_pivots_p, PROTEUS_LAPACK_INTEGER *subdomain_column_pivots_p)

int	nodeStar_free (int N, int subdomain_dim, double subdomain_L, double subdomain_R, double subdomain_U, PROTEUS_LAPACK_INTEGER subdomain_pivots, PROTEUS_LAPACK_INTEGER *subdomain_column_pivots)

int	nodeStar_setU (NodeStarFactorStruct *nodeStarFactor, double val)

int	nodeStar_copy (int other_N, int other_subdomain_dim, double other_subdomain_L, double other_subdomain_R, double other_subdomain_U, PROTEUS_LAPACK_INTEGER other_subdomain_pivots, PROTEUS_LAPACK_INTEGER other_subdomain_column_pivots, int N_p, int subdomain_dim_p, double subdomain_L_p, double subdomain_R_p, double subdomain_U_p, PROTEUS_LAPACK_INTEGER subdomain_pivots_p, PROTEUS_LAPACK_INTEGER *subdomain_column_pivots_p)

void	calculateConservationResidualPWL (int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nDOF_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int dofMapl2g, int nodeStarElements, int nodeStarElementNeighbors, int nElements_node, int fluxElementBoundaries, double elementResidual, double vAverage, double dX, double w, double normal, NodeStarFactorStruct nodeStarFactor, double conservationResidual, double vConservative, double *vConservative_element)

void	calculateConservationJacobianPWL (int nDOF_global, int nNodes_internal, int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nDOF_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int dofMapl2g, int dofStarElements, int dofStarElementNeighbors, int nElements_node, int internalNodes, int fluxElementBoundaries, int fluxBoundaryNodes, double w, double normal, NodeStarFactorStruct *nodeStarFactor)

void	calculateConservationFluxPWL (int nNodes_global, int nNodes_internal, int nElements_node, int internalNodes, int fluxBoundaryNodes, NodeStarFactorStruct nodeStarFactor)

void	calculateConservationFluxPWL_noNeumannFix (int nNodes_global, int nElements_node, NodeStarFactorStruct nodeStarFactor)

void	calculateConservationResidualPWL_opt (int nNodes_owned, int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int nodeStarElements, int nodeStarElementNeighbors, int nElements_node, int fluxElementBoundaries, double elementResidual, double vAverage, double dX, double w, double normal, NodeStarFactorStruct nodeStarFactor, double conservationResidual, double vConservative, double vConservative_element)

void	calculateConservationResidualPWL_primative (int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int skipflag_elementBoundaries, double elementResidual, double dX, double normal, double conservationResidual, double vConservative)

void	calculateConservationJacobianPWL_opt (int nNodes_owned, int nNodes_global, int nNodes_internal, int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int nodeStarElements, int nodeStarElementNeighbors, int nElements_node, int internalNodes, int fluxElementBoundaries, int fluxBoundaryNodes, double w, double normal, NodeStarFactorStruct nodeStarFactor)

void	calculateConservationFluxPWL_opt (int nNodes_owned, int nNodes_global, int nNodes_internal, int nElements_node, int internalNodes, int fluxBoundaryNodes, NodeStarFactorStruct nodeStarFactor)

void	subdomain_U_copy_global2local (int max_nN_owned, int nElements_global, int nNodes_element, int elementNodes, int nodeStarElements, NodeStarFactorStruct nodeStarFactor, double subdomain_U)

void	subdomain_U_copy_local2global (int max_nN_owned, int nElements_global, int nNodes_element, int elementNodes, int nodeStarElements, NodeStarFactorStruct nodeStarFactor, double subdomain_U)

void	updateSelectedExteriorElementBoundaryFlux (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_test_element, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int skipflag_elementBoundaries, double flux, double w_dS, double *residual)
	Update the element boundary flux on exterior element boundaries. More...

void	postprocessAdvectiveVelocityPointEval (int nPoints, int nSpace, double updateCoef, const double f, double velocity)

void	postprocessDiffusiveVelocityPointEval (int nPoints, int nSpace, double updateCoef, const double a, const double grad_phi, double *velocity)

void	postprocessDiffusiveVelocityPointEval_sd (int nPoints, int nSpace, double updateCoef, int rowptr, int colind, const double a, const double grad_phi, double *velocity)

void	calculateElementResidualPWL (int nElements, int nDOF_element_res, int nDOF_element_resPWL, double alpha, double elementResidual, double *elementResidualPWL)

void	calculateConservationResidualPWL_interiorBoundaries (int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int nodeStarElements, int nodeStarElementNeighbors, int nElements_node, int fluxElementBoundaries, double elementResidual, double vAverage, double dX, double w, double normal, NodeStarFactorStruct nodeStarFactor, double conservationResidual, double vConservative, double vConservative_element)

void	calculateConservationJacobianPWL_interiorBoundaries (int nNodes_global, int nElements_global, int nInteriorElementBoundaries_global, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nNodes_element, int nSpace, int interiorElementBoundaries, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int elementNodes, int nodeStarElements, int nodeStarElementNeighbors, int nElements_node, int fluxElementBoundaries, double w, double normal, NodeStarFactorStruct nodeStarFactor)

Function Documentation

◆ buildLocalBDM1projectionMatrices_orig()

void buildLocalBDM1projectionMatrices_orig	(	int	nElements_global,
		int	nElementBoundaries_element,
		int	nQuadraturePoints_elementBoundary,
		int	nSpace,
		int	nDOFs_test_element,
		int	nDOFs_trial_element,
		int	nVDOFs_element,
		double *	w_dS_f,
		double *	ebq_n,
		double *	ebq_v,
		double *	BDMprojectionMat_element
	)

Definition at line 2469 of file postprocessing.c.

◆ calculateConservationFluxPWL_noNeumannFix()

void calculateConservationFluxPWL_noNeumannFix	(	int	nNodes_global,
		int *	nElements_node,
		NodeStarFactorStruct *	nodeStarFactor
	)

Definition at line 4801 of file postprocessing.c.

◆ getElementBDM1velocityValuesLagrangeRep_orig()

void getElementBDM1velocityValuesLagrangeRep_orig	(	int	nElements_global,
		int	nQuadraturePoints_element,
		int	nSpace,
		int	nDOF_trial_element,
		int	nVDOF_element,
		double *	q_v,
		double *	p1_velocity_dofs,
		double *	q_velocity
	)

Definition at line 3345 of file postprocessing.c.

◆ getGlobalElementBoundaryBDM1velocityValuesLagrangeRep()

void getGlobalElementBoundaryBDM1velocityValuesLagrangeRep	(	int	nExteriorElementBoundaries_global,
		int	nQuadraturePoints_elementBoundary,
		int	nSpace,
		int	nDOF_trial_element,
		int	nVDOF_element,
		int *	elementBoundaryElementsArray,
		int *	exteriorElementBoundariesArray,
		double *	ebqe_v,
		double *	p1_velocity_dofs,
		double *	ebq_global_velocity
	)

Definition at line 3551 of file postprocessing.c.

◆ invertLocal()

void invertLocal	(	int	nSpace,
		double	A[3][3],
		double	AI[3][3]
	)