pdd_plast_flow_mixed/weak_formulation/pds_plast_flow_mixed_weakform.c File Reference

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <assert.h>
#include "pdh_intf.h"
#include "pdh_control_intf.h"
#include "mmh_intf.h"
#include "aph_intf.h"
#include "uth_intf.h"
#include "../include/pdh_plast_flow_mixed.h"
#include "../include/pdh_plast_flow_mixed_problem.h"
#include "../include/pdh_plast_flow_mixed_weakform.h"
#include "uth_log.h"

Include dependency graph for pds_plast_flow_mixed_weakform.c:

Functions
int	pdr_plast_flow_select_el_coeff_vect (int Problem_id, int *Coeff_vect_ind)
double *	pdr_plast_flow_select_el_coeff (int Problem_id, double Mval, double Axx, double Axy, double Axz, double Ayx, double Ayy, double Ayz, double Azx, double Azy, double Azz, double Bx, double By, double Bz, double Tx, double Ty, double Tz, double Cval, double Lval, double Qx, double Qy, double Qz, double Sval)
int	pdr_plast_flow_el_coeff (int Problem_id, int El_id, int Mat_num, double Hsize, int Pdeg, double X_loc, double Base_phi, double Base_dphix, double Base_dphiy, double Base_dphiz, double Xcoor, double Uk_val, double Uk_x, double Uk_y, double Uk_z, double Un_val, double Un_x, double Un_y, double Un_z, double Mval, double Axx, double Axy, double Axz, double Ayx, double Ayy, double Ayz, double Azx, double Azy, double Azz, double Bx, double By, double Bz, double Tx, double Ty, double Tz, double Cval, double Lval, double Qx, double Qy, double Qz, double Sval, double Tk, double *Dynamic_viscosity, double Density, double Reference_density, double Reference_velocity, double Delta_t, double Implicitness_coeff, double Body_force_x, double Body_force_y, double Body_force_z)
int	pdr_plast_flow_comp_el_stiff_mat (int Problem_id, int El_id, int Comp_sm, int Pdeg_in, int Nrdofs_loc, double Stiff_mat, double Rhs_vect, char Rewr_dofs)
int	pdr_plast_flow_comp_fa_stiff_mat (int Problem_id, int Fa_id, int Comp_sm, int Pdeg_in, int Nrdofs_loc, double Stiff_mat, double Rhs_vect, char Rewr_dofs)
int	pdr_plast_flow_get_velocity_at_point (int Problem_id, int El_id, double X_loc, double Base_phi, double Base_dphix, double Base_dphiy, double Base_dphiz, double Velocity, double DVel_dx, double DVel_dy, double *DVel_dz)
int	pdr_plast_flow_compute_CFL (int Problem_id, double CFL_min_p, double CFL_max_p, double *CFL_ave_p)

Function Documentation

int pdr_plast_flow_comp_el_stiff_mat	(	int	Problem_id,
		int	El_id,
		int	Comp_sm,
		int	Pdeg_in,
		int *	Nrdofs_loc,
		double *	Stiff_mat,
		double *	Rhs_vect,
		char *	Rewr_dofs
	)

----------------------------------------------------------- pdr_plast_flow_comp_el_stiff_mat - to construct UNCONSTRAINED stiffness matrix and load vector for an element ------------------------------------------------------------

!! coarse element dofs should be supplied by calling routine !!!

----------------------------------------------------------------------!

--------------------- LOOP OVER GAUSS POINTS -------------------------!

----------------------------------------------------------------------!

-------------------- MATERIAL DATA AT GAUSS POINT ----------------- --!

----------------------------------------------------------------------!

------------ CALCULATE BOUSINESSQUE FORCE AT GAUSS POINT -------------!

----------------------------------------------------------------------!

Here is the call graph for this function:

int pdr_plast_flow_comp_fa_stiff_mat	(	int	Problem_id,
		int	Fa_id,
		int	Comp_sm,
		int	Pdeg_in,
		int *	Nrdofs_loc,
		double *	Stiff_mat,
		double *	Rhs_vect,
		char *	Rewr_dofs
	)

----------------------------------------------------------- pdr_plast_flow_comp_fa_stiff_mat - to construct UNCONSTRAINED stiffness matrix and load vector for a face ------------------------------------------------------------

!! need some trick for coarse solve with different Pdeg !!!

DIAGNOSTICS - uncomment !

-------------------------------------------------------------------!

------------------ FLOW BOUNDARY CONDITIONS -----------------------!

-------------------------------------------------------------------!

-------------------------- STRESS PINS --------------------------!

DIAGNOSTICS - uncomment !

-------------------------------------------------------------------!

-------------------------- VELOCITY PINS --------------------------!

DIAGNOSTICS - uncomment !

-------------------------------------------------------------------!

--------------------------- BC_PLAST_FLOW_SYMM ----------------------!

DIAGNOSTICS - uncomment !

-------------------------------------------------------------------!

--------------------------- BC_PLAST_FLOW_FREE_SLIP ----------------------!

DIAGNOSTICS - uncomment !

------------------------------------------------------------------!

-------------------------- BC_PLAST_FLOW_VEL -----------------------!

------------------------ BC_PLAST_FLOW_NOSLIP ----------------------!

DIAGNOSTICS - uncomment !

-------------------------------------------------------------------!

----------------------- BC_PLAST_FLOW_TRACTION -----------------------!

----------------------------------------------------------------------!

------------ CALCULATE EQUIVALENT VISCOSITY AT GAUSS POINT -----------!

----------------------------------------------------------------------!

-------------------------------------------------------------------!

-------------------------- FRICTION ------------------------------!

Here is the call graph for this function:

int pdr_plast_flow_compute_CFL	(	int	Problem_id,
		double *	CFL_min_p,
		double *	CFL_max_p,
		double *	CFL_ave_p
	)

----------------------------------------------------------- pdr_plast_flow_compute_CFL - to compute global CFL numbers (for a subdomain) ------------------------------------------------------------

Here is the call graph for this function:

int pdr_plast_flow_el_coeff	(	int	Problem_id,
		int	Elem,
		int	Mat_num,
		double	Hsize,
		int	Pdeg,
		double *	X_loc,
		double *	Base_phi,
		double *	Base_dphix,
		double *	Base_dphiy,
		double *	Base_dphiz,
		double *	Xcoor,
		double *	Uk_val,
		double *	Uk_x,
		double *	Uk_y,
		double *	Uk_z,
		double *	Un_val,
		double *	Un_x,
		double *	Un_y,
		double *	Un_z,
		double *	Mval,
		double *	Axx,
		double *	Axy,
		double *	Axz,
		double *	Ayx,
		double *	Ayy,
		double *	Ayz,
		double *	Azx,
		double *	Azy,
		double *	Azz,
		double *	Bx,
		double *	By,
		double *	Bz,
		double *	Tx,
		double *	Ty,
		double *	Tz,
		double *	Cval,
		double *	Lval,
		double *	Qx,
		double *	Qy,
		double *	Qz,
		double *	Sval,
		double	Tk,
		double *	Dynamic_viscosity,
		double	Density,
		double	Reference_density,
		double	Reference_velocity,
		double	Delta_t,
		double	Implicitness_coeff,
		double	Body_force_x,
		double	Body_force_y,
		double	Body_force_z
	)

-------------------------------------------------------- pdr_plast_flow_el_coeff - to return coefficients for element integrals ----------------------------------------------------------

----------------------------------------------------------------------!

------------ CALCULATE EQUIVALENT VISCOSITY AT GAUSS POINT -----------!

----------------------------------------------------------------------!

SPLITTING OF TERMS INTO LHS AND RHS FOR ALPHA TIME INTEGRATION SCHEME !

----------------------------------------------------------------------!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!*/

----------------------------------------------------------------------!

------ CALCULATE MATRICES OF WEAK FORM COEFFICIENTS FOR LHS ----------!

----------------------------------------------------------------------!

!! ---- EACH matrix has to have proper dimension so that the corrersponding --- !!!

!! ---- term in the weak formulation has the final dimension [kg/(m*s^3)] ----- !!!

----------------------------------------------------------------------!

------ CALCULATE TIME DERIVATIVE AND STABILIZATION TERMS FOR RHS -----!

----------------------------------------------------------------------!

---------------- CALCULATE BODY FORCES TERMS FOR RHS -----------------!

----------------------------------------------------------------------!

Here is the call graph for this function:

int pdr_plast_flow_get_velocity_at_point	(	int	Problem_id,
		int	El_id,
		double *	X_loc,
		double *	Base_phi,
		double *	Base_dphix,
		double *	Base_dphiy,
		double *	Base_dphiz,
		double *	Velocity,
		double *	DVel_dx,
		double *	DVel_dy,
		double *	DVel_dz
	)

----------------------------------------------------------- pdr_plast_flow_get_velocity_at_point - to provide the velocity and its gradient at a particular point given its local coordinates within an element MODULE PROVIDES IMPLEMENTATION FOR ALL OTHER MODULES (in pds_plast_flow_weakform.c) ------------------------------------------------------------

Here is the call graph for this function:

double* pdr_plast_flow_select_el_coeff	(	int	Problem_id,
		double **	Mval,
		double **	Axx,
		double **	Axy,
		double **	Axz,
		double **	Ayx,
		double **	Ayy,
		double **	Ayz,
		double **	Azx,
		double **	Azy,
		double **	Azz,
		double **	Bx,
		double **	By,
		double **	Bz,
		double **	Tx,
		double **	Ty,
		double **	Tz,
		double **	Cval,
		double **	Lval,
		double **	Qx,
		double **	Qy,
		double **	Qz,
		double **	Sval
	)

!!!!! OLD OBSOLETE VERSION !!!!!!

Here is the call graph for this function:

int pdr_plast_flow_select_el_coeff_vect	(	int	Problem_id,
		int *	Coeff_vect_ind
	)

----------------------------------------------------------- pdr_plast_flow_select_el_coeff_vect - to select coefficients returned to approximation routines for element integrals ------------------------------------------------------------

pdd_plast_flow_mixed/weak_formulation/pds_plast_flow_mixed_weakform.c File Reference

Functions

Function Documentation