-
Notifications
You must be signed in to change notification settings - Fork 96
/
IsotropicLinearMaxwell.i
112 lines (94 loc) · 3.97 KB
/
IsotropicLinearMaxwell.i
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
// -*- C++ -*-
//
// ----------------------------------------------------------------------
//
// Brad T. Aagaard, U.S. Geological Survey
// Charles A. Williams, GNS Science
// Matthew G. Knepley, University at Buffalo
//
// This code was developed as part of the Computational Infrastructure
// for Geodynamics (http://geodynamics.org).
//
// Copyright (c) 2010-2022 University of California, Davis
//
// See LICENSE.md for license information.
//
// ----------------------------------------------------------------------
//
/** @file modulesrc/materials/IsotropicLinearMaxwell.i
*
* Python interface to C++ IsotropicLinearMaxwell.
*/
namespace pylith {
namespace materials {
class IsotropicLinearMaxwell : public pylith::materials::RheologyElasticity {
// PUBLIC METHODS //////////////////////////////////////////////////////////////////////////////////////////
public:
/// Default constructor.
IsotropicLinearMaxwell(void);
/// Destructor.
~IsotropicLinearMaxwell(void);
/// Deallocate PETSc and local data structures.
void deallocate(void);
/** Use reference stress and strain in computation of stress and
* strain?
*
* @param[in] value Flag indicating to include reference stress and strain.
*/
void useReferenceState(const bool value);
/** Use reference stress and strain in computation of stress and
* strain?
*
* @returns True if using reference stress and strain, false otherwise.
*/
bool useReferenceState(void) const;
/** Get auxiliary factory associated with physics.
*
* @return Auxiliary factory for physics object.
*/
pylith::materials::AuxiliaryFactoryElasticity* getAuxiliaryFactory(void);
/** Add rheology subfields to auxiliary field.
*
* @param[inout] auxiliaryField Auxiliary field.
*/
void addAuxiliarySubfields(void);
/** Get stress kernel for RHS residual, G(t,s).
*
* @param[in] coordsys Coordinate system.
*
* @return RHS residual kernel for stress.
*/
PetscPointFunc getKernelResidualStress(const spatialdata::geocoords::CoordSys* coordsys) const;
/** Get elastic constants kernel for RHS Jacobian G(t,s).
*
* @param[in] coordsys Coordinate system.
*
* @return RHS Jacobian kernel for elastic constants.
*/
PetscPointJac getKernelJacobianElasticConstants(const spatialdata::geocoords::CoordSys* coordsys) const;
/** Get stress kernel for derived field.
*
* @param[in] coordsys Coordinate system.
*
* @return Project kernel for computing stress subfield in derived field.
*/
PetscPointFunc getKernelDerivedCauchyStress(const spatialdata::geocoords::CoordSys* coordsys) const;
/** Add kernels for updating state variables.
*
* @param[inout] kernels Array of kernels for updating state variables.
* @param[in] coordsys Coordinate system.
*/
void addKernelsUpdateStateVars(std::vector<pylith::feassemble::IntegratorDomain::ProjectKernels>* kernels,
const spatialdata::geocoords::CoordSys* coordsys) const;
/** Update kernel constants.
*
* @param[inout] kernelConstants Array of constants used in integration kernels.
* @param[in] dt Current time step.
*/
void updateKernelConstants(pylith::real_array* kernelConstants,
const PylithReal dt) const;
};
// class IsotropicLinearMaxwell
} // materials
} // pylith
// End of file