#include <dg.h>

Inheritance diagram for oofem::DGProblem:

Collaboration diagram for oofem::DGProblem:

Classes
class	DGBoundaryEntity

Public Member Functions
	DGProblem (int i, EngngModel *_master)
void	constructBoundaryEntities ()
void	solveYourselfAt (TimeStep *tStep) override
void	updateComponent (TimeStep tStep, NumericalCmpn cmpn, Domain d) override
bool	newDofHandling () override
void	updateSolution (FloatArray &solutionVector, TimeStep tStep, Domain d) override
void	updateInternalRHS (FloatArray &answer, TimeStep tStep, Domain d, FloatArray *eNorm) override
void	updateMatrix (SparseMtrx &mat, TimeStep tStep, Domain d) override
double	giveUnknownComponent (ValueModeType mode, TimeStep tStep, Domain d, Dof *dof) override
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
void	postInitialize () override
virtual void	applyIC ()
int	requiresUnknownsDictionaryUpdate () override
int	giveUnknownDictHashIndx (ValueModeType mode, TimeStep *tStep) override
void	updateDomainLinks () override
Function *	giveDtFunction ()
double	giveDeltaT (int n)
double	giveDiscreteTime (int iStep)
TimeStep *	giveNextStep () override
	Returns next time step (next to current step) of receiver.
TimeStep *	giveSolutionStepWhenIcApply (bool force=false) override
NumericalMethod *	giveNumericalMethod (MetaStep *mStep) override
	Returns reference to receiver's numerical method.
void	initializeFrom (InputRecord &ir) override
bool	requiresEquationRenumbering (TimeStep *tStep) override
int	forceEquationNumbering () override
void	printOutputAt (FILE file, TimeStep tStep) override
void	updateYourself (TimeStep *tStep) override
int	checkConsistency () override
FieldPtr	giveField (FieldType key, TimeStep *tStep) override
const char *	giveInputRecordName () const
const char *	giveClassName () const override
	Returns class name of the receiver.
fMode	giveFormulation () override
double	giveFinalTime () override
	EngngModel (int i, EngngModel *_master=NULL)
virtual	~EngngModel ()
	Destructor.
	EngngModel (const EngngModel &)=delete
EngngModel &	operator= (const EngngModel &)=delete
Domain *	giveDomain (int n)
void	setDomain (int i, Domain *ptr, bool iDeallocateOld=true)
int	giveNumberOfDomains ()
	Returns number of domains in problem.
const std::string &	giveDescription () const
const time_t &	giveStartTime ()
bool	giveSuppressOutput () const
virtual ErrorEstimator *	giveDomainErrorEstimator (int n)
virtual MaterialInterface *	giveMaterialInterface (int n)
void	setNumberOfEquations (int id, int neq)
FILE *	giveOutputStream ()
	Returns file descriptor of output file.
std::string	giveOutputBaseFileName ()
std::string	giveReferenceFileName ()
void	letOutputBaseFileNameBe (const std ::string &src)
ContextOutputMode	giveContextOutputMode () const
int	giveContextOutputStep () const
void	setContextOutputMode (ContextOutputMode contextMode)
void	setUDContextOutputMode (int cStep)
double	giveDeltaT ()
	Returns time step size from the time step controlelr.
void	setDeltaT (double dT)
	Returns time step size through the time step controlelr.
void	setProblemMode (problemMode pmode)
void	setParallelMode (bool newParallelFlag)
problemMode	giveProblemMode () const
	Returns domain mode.
void	setProblemScale (problemScale pscale)
problemScale	giveProblemScale () const
	Returns scale in multiscale simulation.
virtual void	setRenumberFlag ()
	Sets the renumber flag to true.
virtual void	resetRenumberFlag ()
	Sets the renumber flag to false.
double	giveSolutionStepTime ()
void	giveAnalysisTime (int &rhrs, int &rmin, int &rsec, int &uhrs, int &umin, int &usec)
void	terminateAnalysis ()
virtual void	solveYourself ()
virtual void	restartYourself (TimeStep *tS)
virtual void	terminate (TimeStep *tStep)
virtual void	doStepOutput (TimeStep *tStep)
void	saveStepContext (TimeStep *tStep, ContextMode mode)
virtual void	initializeYourself (TimeStep *tStep)
virtual int	initializeAdaptive (int tStepNumber)
virtual int	giveNumberOfDomainEquations (int di, const UnknownNumberingScheme &num)
virtual FieldPtr	giveField (InternalStateType key, TimeStep *)
EngngModel *	giveMasterEngngModel ()
	Returns the master engnmodel.
virtual double	giveLoadLevel ()
	Returns the current load level.
virtual double	giveEigenValue (int eigNum)
	Only relevant for eigen value analysis. Otherwise returns zero.
virtual void	setActiveVector (int i)
	Only relevant for eigen value analysis. Otherwise does noting.
int	updateSharedDofManagers (FloatArray &answer, const UnknownNumberingScheme &s, int ExchangeTag)
int	exchangeRemoteElementData (int ExchangeTag)
virtual int	giveCurrentNumberOfIterations ()
MPI_Comm	giveParallelComm ()
	Returns the communication object of reciever.
int	packRemoteElementData (ProcessCommunicator &processComm)
int	unpackRemoteElementData (ProcessCommunicator &processComm)
int	packDofManagers (ArrayWithNumbering *src, ProcessCommunicator &processComm)
int	unpackDofManagers (ArrayWithNumbering *dest, ProcessCommunicator &processComm)
ProblemCommunicator *	giveProblemCommunicator (EngngModelCommType t)
void	initializeCommMaps (bool forceInit=false)
virtual int	instanciateYourself (DataReader &dr, InputRecord &ir, const char outFileName, const char desc)
void	Instanciate_init ()
int	instanciateDomains (DataReader &dr)
	Instanciate problem domains by calling their instanciateYourself() service.
int	instanciateMetaSteps (DataReader &dr)
	Instanciate problem meta steps by calling their instanciateYourself() service.
virtual int	instanciateDefaultMetaStep (InputRecord &ir)
	Instanciate default metastep, if nmsteps is zero.
virtual void	updateAttributes (MetaStep *mStep)
void	initMetaStepAttributes (MetaStep *mStep)
MetaStep *	giveCurrentMetaStep ()
	Returns current meta step.
virtual TimeStep *	giveCurrentStep (bool force=false)
virtual void	adaptTimeStep (double nIter)
virtual TimeStep *	givePreviousStep (bool force=false)
virtual void	preInitializeNextStep ()
	Does a pre-initialization of the next time step (implement if necessarry).
virtual int	giveNumberOfFirstStep (bool force=false)
int	giveNumberOfMetaSteps ()
	Return number of meta steps.
MetaStep *	giveMetaStep (int i)
	Returns the i-th meta step.
int	giveNumberOfSteps (bool force=false)
virtual double	giveEndOfTimeOfInterest ()
	Returns end of time interest (time corresponding to end of time integration).
int	giveNumberOfTimeStepWhenIcApply ()
	Returns the time step number, when initial conditions should apply.
ExportModuleManager *	giveExportModuleManager ()
	Returns receiver's export module manager.
EngngModelTimer *	giveTimer ()
	Returns reference to receiver timer (EngngModelTimer).
virtual double	giveInitialTime ()
	return time at the begining of analysis
virtual int	giveNewEquationNumber (int domain, DofIDItem)
virtual int	giveNewPrescribedEquationNumber (int domain, DofIDItem)
std::string	giveContextFileName (int tStepNumber, int stepVersion) const
std::string	giveDomainFileName (int domainNum, int domainSerNum) const
virtual void	initForNewIteration (Domain d, TimeStep tStep, int iterationNumber, const FloatArray &solution)
virtual void	initStepIncrements ()
virtual int	forceEquationNumbering (int i)
virtual void	updateDofUnknownsDictionary (DofManager , TimeStep )
virtual ParallelContext *	giveParallelContext (int n)
virtual void	initParallelContexts ()
virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &s, Domain domain)
virtual void	assemble (SparseMtrx &answer, TimeStep tStep, const MatrixAssembler &ma, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s, Domain domain)
void	assembleVector (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromDofManagers (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromElements (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleVectorFromBC (FloatArray &answer, TimeStep tStep, const VectorAssembler &va, ValueModeType mode, const UnknownNumberingScheme &s, Domain domain, FloatArray *eNorms=NULL)
void	assembleExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)
void	assemblePrescribedExtrapolatedForces (FloatArray &answer, TimeStep tStep, CharType type, Domain domain)
virtual int	checkProblemConsistency ()
virtual void	init ()
virtual void	printOutputAt (FILE file, TimeStep tStep, const IntArray &nodeSets, const IntArray &elementSets)
void	outputNodes (FILE file, Domain &domain, TimeStep tStep, int setNum)
void	outputElements (FILE file, Domain &domain, TimeStep tStep, int setNum)
void	printYourself ()
	Prints state of receiver. Useful for debugging.
virtual void	printDofOutputAt (FILE stream, Dof iDof, TimeStep *tStep)
virtual int	useNonlocalStiffnessOption ()
	Returns nonzero if nonlocal stiffness option activated.
bool	isParallel () const
	Returns true if receiver in parallel mode.
int	giveRank () const
	Returns domain rank in a group of collaborating processes (0..groupSize-1).
int	giveNumberOfProcesses () const
	Returns the number of collaborating processes.
EngngModelContext *	giveContext ()
	Context requesting service.
virtual int	giveNumberOfSlaveProblems ()
	Returns number of slave problems.
virtual EngngModel *	giveSlaveProblem (int i)
	Returns i-th slave problem.
virtual bool	giveEquationScalingFlag ()
	Returns the Equation scaling flag, which is used to indicate that governing equation(s) are scaled, or non-dimensionalized.
virtual double	giveVariableScale (VarScaleType varId)
	Returns the scale factor for given variable type.
virtual int	estimateMaxPackSize (IntArray &commMap, DataStream &buff, int packUnpackType)
virtual void	balanceLoad (TimeStep *tStep)
virtual LoadBalancer *	giveLoadBalancer ()
virtual LoadBalancerMonitor *	giveLoadBalancerMonitor ()
void	initParallel ()
	Request domain rank and problem size.
EngngModel *	giveEngngModel ()
	Returns reference to itself -> required by communicator.h.
virtual bool	isElementActivated (int elemNum)
virtual bool	isElementActivated (Element *e)
TimeStepController *	giveTimeStepController ()
	Returns the time step controller.
virtual void	drawYourself (oofegGraphicContext &gc)
virtual void	drawElements (oofegGraphicContext &gc)
virtual void	drawNodes (oofegGraphicContext &gc)
virtual void	showSparseMtrxStructure (int type, oofegGraphicContext &gc, TimeStep *tStep)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Member Functions
std::unique_ptr< Element >	CreateBoundaryElement (Element_Geometry_Type egt, int elemNum, Domain *domain, IntArray &bentityNodes) const
void	assembleDirichletBcRhsVector (FloatArray &answer, TimeStep tStep, ValueModeType mode, const UnknownNumberingScheme &ns, Domain d) const
virtual void	packMigratingData (TimeStep *tStep)
virtual void	unpackMigratingData (TimeStep *tStep)

Protected Attributes
Variable::VariableQuantity	unknownQuantity = Variable::VariableQuantity::VolumeFraction
LinSystSolverType	solverType = ST_Direct
SparseMtrxType	sparseMtrxType = SMT_Skyline
std ::unique_ptr< DofDistributedPrimaryField >	field
std ::unique_ptr< SparseMtrx >	lhsMatrix
std ::unique_ptr< SparseMtrx >	rhsMatrix
FloatArray	solution
FloatArray	internalForces
FloatArray	eNorm
std ::unique_ptr< SparseLinearSystemNM >	nMethod
	Numerical method used to solve the problem.
double	initT = 0.
	Initial time from which the computation runs. Default is zero.
double	deltaT = 0.
	Length of time step.
double	alpha = 0.
int	dtFunction = 0
	Associated time function for time step increment.
FloatArray	prescribedTimes
	Specified times where the problem is solved.
bool	keepTangent = false
bool	hasTangent = false
IntArray	exportFields
std::string	problemType
	identifies what problem to solve (UP, UPV, etc)
bool	preprocessFEM2DG = false
	flag indicating whether the FEM input should be preprocessed to DG problem input (boundary entity generation).
IntArray	sets2preprocess
	array of sets to process to generate boundary node sets
IntArray	targetBoundaryNodeSets
	array of target sets (same size as sets2preprocess) to store the generated boundary node sets
int	targetAllNodeSet = 0
	int target sets for all nodes and all elements
int	targetAllElementSet = 0
std ::vector< std::unique_ptr< DGBoundaryEntity > >	boundaryEntities
Protected Attributes inherited from oofem::EngngModel
int	ndomains
	Number of receiver domains.
std ::vector< std ::unique_ptr< Domain > >	domainList
	List of problem domains.
int	numberOfSteps
	Total number of time steps.
int	numberOfEquations
	Total number of equation in current time step.
int	numberOfPrescribedEquations
	Total number or prescribed equations in current time step.
IntArray	domainNeqs
	Number of equations per domain.
IntArray	domainPrescribedNeqs
	Number of prescribed equations per domain.
bool	renumberFlag
	Renumbering flag (renumbers equations after each step, necessary if Dirichlet BCs change).
bool	profileOpt
	Profile optimized numbering flag (using Sloan's algorithm).
int	equationNumberingCompleted
	Equation numbering completed flag.
int	nMetaSteps
	Number of meta steps.
std ::vector< MetaStep >	metaStepList
	List of problem metasteps.
std ::unique_ptr< TimeStep >	stepWhenIcApply
	Solution step when IC (initial conditions) apply.
std ::unique_ptr< TimeStep >	currentStep
	Current time step.
std ::unique_ptr< TimeStep >	previousStep
	Previous time step.
int	number
	Receivers id.
std::string	dataOutputFileName
	Path to output stream.
std::string	coreOutputFileName
	String with core output file name.
FILE *	outputStream
	Output stream.
std::string	referenceFileName
	String with reference file name.
ContextOutputMode	contextOutputMode
	Domain context output mode.
int	contextOutputStep
ExportModuleManager	exportModuleManager
	Export module manager.
InitModuleManager	initModuleManager
	Initialization module manager.
MonitorManager	monitorManager
	Monitor manager.
problemMode	pMode
	Domain mode.
problemScale	pScale
	Multiscale mode.
time_t	startTime
	Solution start time.
EngngModel *	master
	Master e-model; if defined receiver is in maintained (slave) mode.
EngngModelContext *	context
	Context.
EngngModelTimer	timer
	E-model timer.
int	parallelFlag
	Flag indicating that the receiver runs in parallel.
enum fMode	nonLinFormulation
	Type of non linear formulation (total or updated formulation).
std::unique_ptr< ErrorEstimator >	defaultErrEstimator
	Error estimator. Useful for adaptivity, or simply printing errors output.
std::unique_ptr< TimeStepController >	timeStepController
	Time Step controller is responsible for collecting data from analysis, elements, and materials, and select the appropriate timestep size for the next step, or reduce the step in case of convergence problems.
int	rank
	Domain rank in a group of collaborating processes (0..groupSize-1).
int	numProcs
	Total number of collaborating processes.
int	nonlocalExt
	Flag indicating if nonlocal extension active, which will cause data to be sent between shared elements before computing the internal forces.
char	processor_name [PROCESSOR_NAME_LENGTH]
	Processor name.
MPI_Comm	comm
	Communication object for this engineering model.
std::unique_ptr< LoadBalancer >	lb
	Load Balancer.
std::unique_ptr< LoadBalancerMonitor >	lbm
bool	loadBalancingFlag
	If set to true, load balancing is active.
bool	force_load_rebalance_in_first_step
	Debug flag forcing load balancing after first step.
CommunicatorBuff *	commBuff
	Common Communicator buffer.
ProblemCommunicator *	communicator
	Communicator.
ProblemCommunicator *	nonlocCommunicator
	NonLocal Communicator. Necessary when nonlocal constitutive models are used.
std ::vector< ParallelContext >	parallelContextList
	List where parallel contexts are stored.
bool	suppressOutput
	Flag for suppressing output to file.
std::string	simulationDescription

Additional Inherited Members
Public Types inherited from oofem::EngngModel
enum	EngngModel_UpdateMode { EngngModel_Unknown_Mode , EngngModel_SUMM_Mode , EngngModel_SET_Mode }
enum	EngngModelCommType { PC_default , PC_nonlocal }
enum	InitialGuess { IG_None = 0 , IG_Tangent = 1 }
Protected Types inherited from oofem::EngngModel
enum	{ InternalForcesExchangeTag , MassExchangeTag , LoadExchangeTag , ReactionExchangeTag , RemoteElementExchangeTag }
	Message tags. More...

Detailed Description

This class represents generic Discontinuous Galerkin solver. The problem can be growing/decreasing, signalized by flag "changingproblemsize" in the problem description. The solution is stored in UnknownsField, which can obtain/ project solution from/to DOFs (nodes). If the problem keeps the same equation numbers, solution is taken from UnknownsField without any projection, which is more efficient. See the matlibmanual for solution strategy of balance equations and the solution algorithm.

Todo: Documentation errors (there is no "UnknownsField" used here).

Definition at line 130 of file dg.h.

Constructor & Destructor Documentation

◆ DGProblem()

oofem::DGProblem::DGProblem	(	int	i,
		EngngModel *	_master = nullptr )

Definition at line 153 of file dg.C.

References oofem::EngngModel::EngngModel(), and oofem::EngngModel::ndomains.

Member Function Documentation

◆ applyIC()

void oofem::DGProblem::applyIC ( )

virtual

Definition at line 650 of file dg.C.

References field, oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), giveSolutionStepWhenIcApply(), OOFEM_LOG_INFO, oofem::Element::updateInternalState(), and oofem::Element::updateYourself().

Referenced by solveYourselfAt().

◆ assembleDirichletBcRhsVector()

void oofem::DGProblem::assembleDirichletBcRhsVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		ValueModeType	mode,
		const UnknownNumberingScheme &	ns,
		Domain *	d ) const

protected

@BP: This method is used to assemble the Dirichlet boundary conditions into the RHS vector. Note that it assembles the contributions from Lhs and Rhs! This is needed, as the Rhs is at present evaluated as Matrix * PreviousSolutionVector, and the Dirichlet BCs are not included in previous solution vector!

Definition at line 787 of file dg.C.

References oofem::FloatMatrix::add(), oofem::FloatArray::assemble(), oofem::FloatArray::beProductOf(), oofem::Element::computeVectorOfPrescribed(), oofem::FloatArray::containsOnlyZeroes(), deltaT, oofem::Element::giveCharacteristicMatrix(), oofem::Domain::giveElement(), oofem::Element::giveElementDofIDMask(), oofem::Element::giveLocationArray(), oofem::Domain::giveNumberOfElements(), oofem::TimeStep::givePreviousStep(), oofem::FloatMatrix::isNotEmpty(), oofem::FloatArray::negated(), and oofem::FloatMatrix::times().

Referenced by solveYourselfAt().

◆ checkConsistency()

int oofem::DGProblem::checkConsistency ( )

overridevirtual

Allows programmer to test some receiver's internal data, before computation begins.

Returns: Nonzero if receiver check is o.k.

Reimplemented from oofem::EngngModel.

Definition at line 746 of file dg.C.

◆ constructBoundaryEntities()

void oofem::DGProblem::constructBoundaryEntities ( )

Definition at line 305 of file dg.C.

References oofem::IntArray::at(), boundaryEntities, oofem::Set::clear(), CreateBoundaryElement(), oofem::IntArray::findFirstIndexOf(), oofem::IntArray::findSorted(), oofem::IntArray::followedBy(), oofem::FEInterpolation::giveBoundaryGeometryType(), oofem::Set::giveBoundaryList(), oofem::Element::giveBoundaryNodes(), oofem::DofManager::giveCoordinates(), oofem::Domain::giveDofManager(), oofem::Element::giveDofManager(), oofem::Element::giveDofManagerNumber(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElement(), oofem::ConnectivityTable::giveElementsWithNodes(), oofem::Element::giveInterpolation(), oofem::FEMComponent::giveNumber(), oofem::Element::giveNumberOfBoundarySides(), oofem::Domain::giveNumberOfDofManagers(), oofem::Element::giveNumberOfDofManagers(), oofem::Domain::giveNumberOfElements(), oofem::Domain::giveSet(), oofem::IntArray::giveSize(), OOFEM_ERROR, OOFEM_LOG_INFO, oofem::IntArray::resize(), oofem::Domain::resizeDofManagers(), oofem::Domain::resizeElements(), oofem::Domain::setDofManager(), oofem::Element::setDofManagers(), oofem::Domain::setElement(), oofem::Set::setElementList(), oofem::Set::setNodeList(), sets2preprocess, oofem::IntArray::sort(), targetAllElementSet, targetAllNodeSet, targetBoundaryNodeSets, and oofem::EngngModel::timer.

Referenced by postInitialize().

◆ CreateBoundaryElement()

std::unique_ptr< Element > oofem::DGProblem::CreateBoundaryElement	(	Element_Geometry_Type	egt,
		int	elemNum,
		Domain *	domain,
		IntArray &	bentityNodes ) const

protected

Definition at line 500 of file dg.C.

References oofem::classFactory, and OOFEM_ERROR.

Referenced by constructBoundaryEntities().

◆ forceEquationNumbering()

int oofem::DGProblem::forceEquationNumbering ( )

overridevirtual

Forces equation renumbering on all domains associated to engng model. All equation numbers in all domains for all dofManagers are invalidated, and new equation numbers are generated starting from 1 on each domain. It will update numberOfEquations variable accordingly. Should be used at startup to force equation numbering and therefore sets numberOfEquations. Must be used if model supports changes of static system to assign new valid equation numbers to dofManagers.

Reimplemented from oofem::EngngModel.

Definition at line 691 of file dg.C.

References lhsMatrix, and rhsMatrix.

◆ giveClassName()

const char * oofem::DGProblem::giveClassName ( ) const

inlineoverridevirtual

Returns class name of the receiver.

Implements oofem::EngngModel.

Definition at line 232 of file dg.h.

◆ giveDeltaT()

double oofem::DGProblem::giveDeltaT ( int n )

Definition at line 248 of file dg.C.

References deltaT, dtFunction, giveDiscreteTime(), oofem::EngngModel::giveDomain(), and prescribedTimes.

Referenced by giveNextStep(), and giveSolutionStepWhenIcApply().

◆ giveDiscreteTime()

double oofem::DGProblem::giveDiscreteTime ( int iStep )

Definition at line 260 of file dg.C.

References initT, OOFEM_ERROR, and prescribedTimes.

Referenced by giveDeltaT().

◆ giveDtFunction()

Function * oofem::DGProblem::giveDtFunction ( )

◆ giveField()

FieldPtr oofem::DGProblem::giveField	(	FieldType	key,
		TimeStep *	)

overridevirtual

Returns the smart pointer to requested field, Null otherwise. The return value uses shared_ptr, as some registered fields may be owned (and maintained) by emodel, while some may be created on demand and thus reliable reference counting mechanism is essential.

Reimplemented from oofem::EngngModel.

Definition at line 760 of file dg.C.

References field, oofem::EngngModel::giveContext(), oofem::EngngModel::giveCurrentStep(), oofem::FieldManager::giveField(), oofem::EngngModelContext::giveFieldManager(), and OOFEM_ERROR.

◆ giveFinalTime()

double oofem::DGProblem::giveFinalTime ( )

inlineoverridevirtual

Returns final time of the simulation

Reimplemented from oofem::EngngModel.

Definition at line 236 of file dg.h.

References deltaT, oofem::EngngModel::numberOfSteps, and prescribedTimes.

◆ giveFormulation()

fMode oofem::DGProblem::giveFormulation ( )

inlineoverridevirtual

Indicates type of non linear computation (total or updated formulation). This is used for example on Nodal level to update coordinates if updated formulation is done, or on element level, when non linear contributions are computed.

Reimplemented from oofem::EngngModel.

Definition at line 233 of file dg.h.

References oofem::TL.

◆ giveInputRecordName()

const char * oofem::DGProblem::giveInputRecordName ( ) const

inline

Definition at line 231 of file dg.h.

References _IFT_DGProblem_Name.

◆ giveNextStep()

TimeStep * oofem::DGProblem::giveNextStep ( )

overridevirtual

Returns next time step (next to current step) of receiver.

Reimplemented from oofem::EngngModel.

Definition at line 273 of file dg.C.

References alpha, oofem::EngngModel::currentStep, giveDeltaT(), giveSolutionStepWhenIcApply(), and oofem::EngngModel::previousStep.

◆ giveNumericalMethod()

NumericalMethod * oofem::DGProblem::giveNumericalMethod ( MetaStep * mStep )

overridevirtual

Returns reference to receiver's numerical method.

Reimplemented from oofem::EngngModel.

Definition at line 158 of file dg.C.

References oofem::classFactory, oofem::EngngModel::giveDomain(), oofem::EngngModel::isParallel(), nMethod, OOFEM_ERROR, solverType, oofem::ST_Feti, and oofem::ST_Petsc.

Referenced by solveYourselfAt(), and updateDomainLinks().

◆ giveSolutionStepWhenIcApply()

TimeStep * oofem::DGProblem::giveSolutionStepWhenIcApply ( bool force = false )

overridevirtual

Returns the solution step when Initial Conditions (IC) apply.

Parameters

force when set to true then receiver reply is returned instead of master (default)

Reimplemented from oofem::EngngModel.

Definition at line 288 of file dg.C.

References alpha, giveDeltaT(), oofem::EngngModel::giveNumberOfTimeStepWhenIcApply(), initT, oofem::EngngModel::master, and oofem::EngngModel::stepWhenIcApply.

Referenced by applyIC(), and giveNextStep().

◆ giveUnknownComponent()

double oofem::DGProblem::giveUnknownComponent	(	ValueModeType	,
		TimeStep *	,
		Domain *	,
		Dof *	)

overridevirtual

Returns requested unknown. Unknown at give time step is characterized by its type and mode and by its equation number. This function is used by Dofs, when they are requested for their associated unknowns.

See also: Dof::giveUnknown

Reimplemented from oofem::EngngModel.

Definition at line 241 of file dg.C.

References field.

◆ giveUnknownDictHashIndx()

int oofem::DGProblem::giveUnknownDictHashIndx	(	ValueModeType	mode,
		TimeStep *	tStep )

overridevirtual

This method is responsible for computing unique dictionary id (ie hash value) from given valueModeType and time step. This function is used by particular dofs to access unknown identified by given parameters from its dictionary using computed index. Usually the hash algorithm should produce index that depend on time step relatively to actual one to avoid storage of complete history.

Reimplemented from oofem::EngngModel.

Definition at line 733 of file dg.C.

References oofem::TimeStep::giveNumber().

◆ initializeFrom()

void oofem::DGProblem::initializeFrom ( InputRecord & ir )

overridevirtual

Initializes receiver according to object description in input reader. InitString can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.

Reimplemented from oofem::EngngModel.

Definition at line 179 of file dg.C.

References _IFT_DGProblem_alpha, _IFT_DGProblem_deltat, _IFT_DGProblem_deltatfunction, _IFT_DGProblem_exportFields, _IFT_DGProblem_initt, _IFT_DGProblem_keepTangent, _IFT_DGProblem_preprocessFEM2DG, _IFT_DGProblem_prescribedtimes, _IFT_DGProblem_problemType, _IFT_DGProblem_sets2preprocess, _IFT_DGProblem_targetAllElementSet, _IFT_DGProblem_targetAllNodeSet, _IFT_DGProblem_targetBoundaryNodeSets, _IFT_EngngModel_smtype, alpha, deltaT, dtFunction, exportFields, field, oofem::EngngModel::giveContext(), oofem::InputRecord::hasField(), initT, IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, keepTangent, OOFEM_ERROR, OOFEM_LOG_RELEVANT, preprocessFEM2DG, prescribedTimes, problemType, oofem::FieldManager::registerField(), sets2preprocess, oofem::SMT_Skyline, sparseMtrxType, targetAllElementSet, targetAllNodeSet, and targetBoundaryNodeSets.

◆ newDofHandling()

bool oofem::DGProblem::newDofHandling ( )

inlineoverridevirtual

Temporary method for allowing code to seamlessly convert from the old to new way of handling DOF values. (the new way expects the field to store all values, regardless of if they are computed, from BC, or IC.) This is used by MasterDof

Todo: When all models have converted to using a field, this should be removed.

Reimplemented from oofem::EngngModel.

Definition at line 196 of file dg.h.

◆ postInitialize()

void oofem::DGProblem::postInitialize ( )

overridevirtual

Performs post-initialization for all the problem contents (which is called after initializeFrom). Currently, it calls Domain::postInitialize for all problem domains.

Reimplemented from oofem::EngngModel.

Definition at line 518 of file dg.C.

References constructBoundaryEntities(), oofem::EngngModel::domainList, oofem::Element::giveClassName(), oofem::Element::giveDofManager(), oofem::EngngModel::giveDomain(), oofem::ClonedDofManager::giveMasterNumber(), oofem::Set::giveNodeList(), oofem::FEMComponent::giveNumber(), oofem::Element::giveNumberOfDofManagers(), oofem::EngngModel::giveNumberOfFirstStep(), oofem::EngngModel::master, oofem::EngngModel::metaStepList, preprocessFEM2DG, and targetBoundaryNodeSets.

◆ printOutputAt()

void oofem::DGProblem::printOutputAt	(	FILE *	file,
		TimeStep *	tStep )

overridevirtual

Prints output of receiver to output domain stream, for given time step. Corresponding function for element gauss points is invoked (gaussPoint::printOutputAt).

Reimplemented from oofem::EngngModel.

Definition at line 700 of file dg.C.

References oofem::EngngModel::giveDomain().

◆ requiresEquationRenumbering()

bool oofem::DGProblem::requiresEquationRenumbering ( TimeStep * tStep )

overridevirtual

Returns true if equation renumbering is required for given solution step. This may of course change the number of equation and in general there is no guarantee that for a certain dof the same equation will be assigned. So the use of DOF unknowns dictionaries is generally recommended.

Todo: This method should be set as the default behavior instead of relying on a user specified flag. Then this function should be removed.

Reimplemented from oofem::EngngModel.

Definition at line 668 of file dg.C.

References oofem::Domain::giveBcs(), oofem::EngngModel::giveDomain(), oofem::GeneralBoundaryCondition::giveNumberOfInternalDofManagers(), oofem::TimeStep::givePreviousStep(), oofem::TimeStep::isTheFirstStep(), and oofem::ActiveBoundaryCondition::requiresActiveDofs().

◆ requiresUnknownsDictionaryUpdate()

int oofem::DGProblem::requiresUnknownsDictionaryUpdate ( )

overridevirtual

Indicates if EngngModel requires Dofs dictionaries to be updated. If EngngModel does not support changes of static system, the dof forwards the requests for its unknowns to EngngModel, where unknowns are naturally kept. This is possible, because dof equation number is same during whole solution. But when changes of static system are allowed, several problem arise. For example by solving simple incremental static with allowed static changes, the incremental displacement vector of structure can not be added to total displacement vector of structure, because equation numbers may have changed, and one can not simply add these vector to obtain new total displacement vector, because incompatible displacement will be added. To solve this problem, unknown dictionary at dof level has been assumed. Dof then keeps its unknowns in its own private dictionary. After computing increment of solution, engngModel updates for each dof its unknowns in its dictionary (using updateUnknownsDictionary function). For aforementioned example engngModel updates incremental values but also total value by asking dof for previous total value (dof will use its dictionary, does not asks back EngngModel) adds corresponding increment and updates total value in dictionary.

Reimplemented from oofem::EngngModel.

Definition at line 740 of file dg.C.

◆ restoreContext()

void oofem::DGProblem::restoreContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

Restores the state of model from output stream. Restores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by restoring element context also contexts of all associated integration points (and material statuses) are restored. Each context is associated with unique time step. Only one context per time step is allowed. Restore context function will restore such context, which is related (through its step number) to time step number and version given in obj parameter. Restoring context will change current time step in order to correspond to newly restored context.

Parameters

stream	Context file.
mode	Determines amount of info in stream.

Exceptions

ContextIOERR exception if error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 725 of file dg.C.

References field.

◆ saveContext()

void oofem::DGProblem::saveContext	(	DataStream &	stream,
		ContextMode	mode )

overridevirtual

Stores the state of model to output stream. Stores not only the receiver state, but also same function is invoked for all DofManagers and Elements in associated domain. Note that by storing element context also contexts of all associated integration points (and material statuses) are stored.

Parameters

stream	Context stream.
mode	Determines amount of info in stream.

Exceptions

ContextIOERR If error encountered.

Reimplemented from oofem::EngngModel.

Definition at line 717 of file dg.C.

References field.

◆ solveYourselfAt()

void oofem::DGProblem::solveYourselfAt ( TimeStep * tStep )

overridevirtual

Solves problem for given time step. Should assemble characteristic matrices and vectors if necessary and solve problem using appropriate numerical method. After finishing solution, this->updateYourself function for updating solution state and then this->terminate function (for updating nodal and element values) should be called.

Reimplemented from oofem::EngngModel.

Definition at line 578 of file dg.C.

References alpha, applyIC(), oofem::assemble(), assembleDirichletBcRhsVector(), oofem::EngngModel::assembleVector(), oofem::classFactory, oofem::TimeStep::convergedReason, field, oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), giveNumericalMethod(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::isTheFirstStep(), lhsMatrix, nMethod, OOFEM_LOG_INFO, solution, sparseMtrxType, unknownQuantity, and updateSolution().

◆ updateComponent()

void oofem::DGProblem::updateComponent	(	TimeStep *	tStep,
		NumericalCmpn	cmpn,
		Domain *	d )

overridevirtual

Reimplemented from oofem::EngngModel.

Definition at line 643 of file dg.C.

References OOFEM_ERROR.

◆ updateDomainLinks()

void oofem::DGProblem::updateDomainLinks ( )

overridevirtual

Updates domain links after the domains of receiver have changed. Used mainly after restoring context - the domains may change and this service is then used to update domain variables in all components belonging to receiver like error estimators, solvers, etc, having domains as attributes.

Reimplemented from oofem::EngngModel.

Definition at line 753 of file dg.C.

References oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), and giveNumericalMethod().

◆ updateInternalRHS()

void oofem::DGProblem::updateInternalRHS	(	FloatArray &	answer,
		TimeStep *	tStep,
		Domain *	d,
		FloatArray *	eNorm )

overridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson).

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.
eNorm	Optional per-element norm (for normalization).

Reimplemented from oofem::EngngModel.

Definition at line 633 of file dg.C.

References eNorm.

◆ updateMatrix()

void oofem::DGProblem::updateMatrix	(	SparseMtrx &	mat,
		TimeStep *	tStep,
		Domain *	d )

overridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson).

Note: For performance, the matrix should keep it's non-zero structure between calls, so the caller should make sure not to clear the matrix object before called.

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.

Reimplemented from oofem::EngngModel.

Definition at line 638 of file dg.C.

◆ updateSolution()

void oofem::DGProblem::updateSolution	(	FloatArray &	solutionVector,
		TimeStep *	tStep,
		Domain *	d )

overridevirtual

Updates the solution (guess) according to the new values. Callback for nonlinear solvers (e.g. Newton-Raphson), and are called before new internal forces are computed.

Parameters

solutionVector	New solution.
tStep	Time when component is updated.
d	Domain.

Todo: NRSolver should report when the solution changes instead of doing it this way.

Todo: Need to reset the boundary conditions properly since some "update" is doing strange things such as applying the (wrong) boundary conditions. This call will be removed when that code can be removed.

Reimplemented from oofem::EngngModel.

Definition at line 622 of file dg.C.

References field.

Referenced by solveYourselfAt().

◆ updateYourself()

void oofem::DGProblem::updateYourself ( TimeStep * tStep )

overridevirtual

Updates internal state after finishing time step. (for example total values may be updated according to previously solved increments). Then element values are also updated (together with related integration points and material statuses).

Reimplemented from oofem::EngngModel.

Definition at line 711 of file dg.C.

Member Data Documentation

◆ alpha

double oofem::DGProblem::alpha = 0.

protected

Definition at line 154 of file dg.h.

Referenced by giveNextStep(), giveSolutionStepWhenIcApply(), initializeFrom(), and solveYourselfAt().

◆ boundaryEntities

std :: vector< std::unique_ptr<DGBoundaryEntity> > oofem::DGProblem::boundaryEntities

protected

Definition at line 185 of file dg.h.

Referenced by constructBoundaryEntities().

◆ deltaT

double oofem::DGProblem::deltaT = 0.

protected

Length of time step.

Definition at line 153 of file dg.h.

Referenced by assembleDirichletBcRhsVector(), giveDeltaT(), giveFinalTime(), and initializeFrom().

◆ dtFunction

int oofem::DGProblem::dtFunction = 0

protected

Associated time function for time step increment.

Definition at line 156 of file dg.h.

Referenced by giveDeltaT(), and initializeFrom().

◆ eNorm

FloatArray oofem::DGProblem::eNorm

protected

Definition at line 145 of file dg.h.

Referenced by updateInternalRHS().

◆ exportFields

IntArray oofem::DGProblem::exportFields

protected

Definition at line 160 of file dg.h.

Referenced by initializeFrom().

◆ field

std :: unique_ptr< DofDistributedPrimaryField > oofem::DGProblem::field

protected

Definition at line 138 of file dg.h.

Referenced by applyIC(), giveField(), giveUnknownComponent(), initializeFrom(), restoreContext(), saveContext(), solveYourselfAt(), and updateSolution().

◆ hasTangent

bool oofem::DGProblem::hasTangent = false

protected

Definition at line 159 of file dg.h.

◆ initT

double oofem::DGProblem::initT = 0.

protected

Initial time from which the computation runs. Default is zero.

Definition at line 151 of file dg.h.

Referenced by giveDiscreteTime(), giveSolutionStepWhenIcApply(), and initializeFrom().

◆ internalForces

FloatArray oofem::DGProblem::internalForces

protected

Definition at line 144 of file dg.h.

◆ keepTangent

bool oofem::DGProblem::keepTangent = false

protected

Definition at line 159 of file dg.h.

Referenced by initializeFrom().

◆ lhsMatrix

std :: unique_ptr< SparseMtrx > oofem::DGProblem::lhsMatrix

protected

Definition at line 140 of file dg.h.

Referenced by forceEquationNumbering(), and solveYourselfAt().

◆ nMethod

std :: unique_ptr< SparseLinearSystemNM > oofem::DGProblem::nMethod

protected

Numerical method used to solve the problem.

Definition at line 148 of file dg.h.

Referenced by giveNumericalMethod(), and solveYourselfAt().

◆ preprocessFEM2DG

bool oofem::DGProblem::preprocessFEM2DG = false

protected

flag indicating whether the FEM input should be preprocessed to DG problem input (boundary entity generation).

Definition at line 164 of file dg.h.

Referenced by initializeFrom(), and postInitialize().

◆ prescribedTimes

FloatArray oofem::DGProblem::prescribedTimes

protected

Specified times where the problem is solved.

Definition at line 158 of file dg.h.

Referenced by giveDeltaT(), giveDiscreteTime(), giveFinalTime(), and initializeFrom().

◆ problemType

std::string oofem::DGProblem::problemType

protected

identifies what problem to solve (UP, UPV, etc)

Definition at line 162 of file dg.h.

Referenced by initializeFrom().

◆ rhsMatrix

std :: unique_ptr< SparseMtrx > oofem::DGProblem::rhsMatrix

protected

Definition at line 141 of file dg.h.

Referenced by forceEquationNumbering().

◆ sets2preprocess

IntArray oofem::DGProblem::sets2preprocess

protected

array of sets to process to generate boundary node sets

Definition at line 166 of file dg.h.

Referenced by constructBoundaryEntities(), and initializeFrom().

◆ solution

FloatArray oofem::DGProblem::solution

protected

Definition at line 143 of file dg.h.

Referenced by solveYourselfAt().

◆ solverType

LinSystSolverType oofem::DGProblem::solverType = ST_Direct

protected

Definition at line 136 of file dg.h.

Referenced by giveNumericalMethod().

◆ sparseMtrxType

SparseMtrxType oofem::DGProblem::sparseMtrxType = SMT_Skyline

protected

Definition at line 137 of file dg.h.

Referenced by initializeFrom(), and solveYourselfAt().

◆ targetAllElementSet

int oofem::DGProblem::targetAllElementSet = 0

protected

Definition at line 170 of file dg.h.

Referenced by constructBoundaryEntities(), and initializeFrom().

◆ targetAllNodeSet

int oofem::DGProblem::targetAllNodeSet = 0

protected

int target sets for all nodes and all elements

Definition at line 170 of file dg.h.

Referenced by constructBoundaryEntities(), and initializeFrom().

◆ targetBoundaryNodeSets

IntArray oofem::DGProblem::targetBoundaryNodeSets

protected

array of target sets (same size as sets2preprocess) to store the generated boundary node sets

Definition at line 168 of file dg.h.

Referenced by constructBoundaryEntities(), initializeFrom(), and postInitialize().

◆ unknownQuantity

Variable::VariableQuantity oofem::DGProblem::unknownQuantity = Variable::VariableQuantity::VolumeFraction

protected

Definition at line 134 of file dg.h.

Referenced by solveYourselfAt().

The documentation for this class was generated from the following files:

dg.h
dg.C

Classes

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DGProblem()

Member Function Documentation

◆ applyIC()

◆ assembleDirichletBcRhsVector()

◆ checkConsistency()

◆ constructBoundaryEntities()

◆ CreateBoundaryElement()

◆ forceEquationNumbering()

◆ giveClassName()

◆ giveDeltaT()

◆ giveDiscreteTime()

◆ giveDtFunction()

◆ giveField()

◆ giveFinalTime()

◆ giveFormulation()

◆ giveInputRecordName()

◆ giveNextStep()

◆ giveNumericalMethod()

◆ giveSolutionStepWhenIcApply()

◆ giveUnknownComponent()

◆ giveUnknownDictHashIndx()

◆ initializeFrom()

◆ newDofHandling()

◆ postInitialize()

◆ printOutputAt()

◆ requiresEquationRenumbering()

◆ requiresUnknownsDictionaryUpdate()

◆ restoreContext()

◆ saveContext()

◆ solveYourselfAt()

◆ updateComponent()

◆ updateDomainLinks()

◆ updateInternalRHS()

◆ updateMatrix()

◆ updateSolution()

◆ updateYourself()

Member Data Documentation

◆ alpha

◆ boundaryEntities

◆ deltaT

◆ dtFunction

◆ eNorm

◆ exportFields

◆ field

◆ hasTangent

◆ initT

◆ internalForces

◆ keepTangent

◆ lhsMatrix

◆ nMethod

◆ preprocessFEM2DG

◆ prescribedTimes

◆ problemType

◆ rhsMatrix

◆ sets2preprocess

◆ solution

◆ solverType

◆ sparseMtrxType

◆ targetAllElementSet

◆ targetAllNodeSet

◆ targetBoundaryNodeSets

◆ unknownQuantity