MFEM
v3.4
Finite element discretization library
|
#include <maxwell_solver.hpp>
Public Member Functions | |
MaxwellSolver (ParMesh &pmesh, int sOrder, double(*eps)(const Vector &), double(*muInv)(const Vector &), double(*sigma)(const Vector &), void(*j_src)(const Vector &, double, Vector &), Array< int > &abcs, Array< int > &dbcs, void(*dEdt_bc)(const Vector &, double, Vector &)) | |
~MaxwellSolver () | |
int | GetLogging () const |
void | SetLogging (int logging) |
HYPRE_Int | GetProblemSize () |
void | PrintSizes () |
void | SetInitialEField (VectorCoefficient &EFieldCoef) |
void | SetInitialBField (VectorCoefficient &BFieldCoef) |
void | Mult (const Vector &B, Vector &dEdt) const |
Perform the action of the operator: y = k = f(x, t), where k solves the algebraic equation F(x, k, t) = G(x, t) and t is the current time. More... | |
void | ImplicitSolve (const double dt, const Vector &x, Vector &k) |
Solve the equation: k = f(x + dt k, t), for the unknown k at the current time t. More... | |
double | GetMaximumTimeStep () const |
double | GetEnergy () const |
Operator & | GetNegCurl () |
Vector & | GetEField () |
Vector & | GetBField () |
void | SyncGridFuncs () |
void | RegisterVisItFields (VisItDataCollection &visit_dc) |
void | WriteVisItFields (int it=0) |
void | InitializeGLVis () |
void | DisplayToGLVis () |
Public Member Functions inherited from mfem::TimeDependentOperator | |
TimeDependentOperator (int n=0, double t_=0.0, Type type_=EXPLICIT) | |
Construct a "square" TimeDependentOperator y = f(x,t), where x and y have the same dimension n. More... | |
TimeDependentOperator (int h, int w, double t_=0.0, Type type_=EXPLICIT) | |
Construct a TimeDependentOperator y = f(x,t), where x and y have dimensions w and h, respectively. More... | |
virtual double | GetTime () const |
Read the currently set time. More... | |
virtual void | SetTime (const double _t) |
Set the current time. More... | |
bool | isExplicit () const |
True if type is EXPLICIT. More... | |
bool | isImplicit () const |
True if type is IMPLICIT or HOMOGENEOUS. More... | |
bool | isHomogeneous () const |
True if type is HOMOGENEOUS. More... | |
virtual void | ExplicitMult (const Vector &x, Vector &y) const |
Perform the action of the explicit part of the operator, G: y = G(x, t) where t is the current time. More... | |
virtual void | ImplicitMult (const Vector &x, const Vector &k, Vector &y) const |
Perform the action of the implicit part of the operator, F: y = F(x, k, t) where t is the current time. More... | |
virtual Operator & | GetImplicitGradient (const Vector &x, const Vector &k, double shift) const |
Return an Operator representing (dF/dk shift + dF/dx) at the given x, k, and the currently set time. More... | |
virtual Operator & | GetExplicitGradient (const Vector &x) const |
Return an Operator representing dG/dx at the given point x and the currently set time. More... | |
virtual | ~TimeDependentOperator () |
Public Member Functions inherited from mfem::Operator | |
Operator (int s=0) | |
Construct a square Operator with given size s (default 0). More... | |
Operator (int h, int w) | |
Construct an Operator with the given height (output size) and width (input size). More... | |
int | Height () const |
Get the height (size of output) of the Operator. Synonym with NumRows(). More... | |
int | NumRows () const |
Get the number of rows (size of output) of the Operator. Synonym with Height(). More... | |
int | Width () const |
Get the width (size of input) of the Operator. Synonym with NumCols(). More... | |
int | NumCols () const |
Get the number of columns (size of input) of the Operator. Synonym with Width(). More... | |
virtual void | MultTranspose (const Vector &x, Vector &y) const |
Action of the transpose operator: y=A^t(x) . The default behavior in class Operator is to generate an error. More... | |
virtual Operator & | GetGradient (const Vector &x) const |
Evaluate the gradient operator at the point x. The default behavior in class Operator is to generate an error. More... | |
virtual const Operator * | GetProlongation () const |
Prolongation operator from linear algebra (linear system) vectors, to input vectors for the operator. NULL means identity. More... | |
virtual const Operator * | GetRestriction () const |
Restriction operator from input vectors for the operator to linear algebra (linear system) vectors. NULL means identity. More... | |
void | FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, Vector &b, Operator *&A, Vector &X, Vector &B, int copy_interior=0) |
Form a constrained linear system using a matrix-free approach. More... | |
virtual void | RecoverFEMSolution (const Vector &X, const Vector &b, Vector &x) |
Reconstruct a solution vector x (e.g. a GridFunction) from the solution X of a constrained linear system obtained from Operator::FormLinearSystem(). More... | |
void | PrintMatlab (std::ostream &out, int n=0, int m=0) const |
Prints operator with input size n and output size m in Matlab format. More... | |
virtual | ~Operator () |
Virtual destructor. More... | |
Type | GetType () const |
Return the type ID of the Operator class. More... | |
Additional Inherited Members | |
Public Types inherited from mfem::TimeDependentOperator | |
enum | Type { EXPLICIT, IMPLICIT, HOMOGENEOUS } |
Public Types inherited from mfem::Operator | |
enum | Type { ANY_TYPE, MFEM_SPARSEMAT, Hypre_ParCSR, PETSC_MATAIJ, PETSC_MATIS, PETSC_MATSHELL, PETSC_MATNEST, PETSC_MATHYPRE, PETSC_MATGENERIC } |
Enumeration defining IDs for some classes derived from Operator. More... | |
Protected Attributes inherited from mfem::TimeDependentOperator | |
double | t |
Current time. More... | |
Type | type |
Describes the form of the TimeDependentOperator. More... | |
Protected Attributes inherited from mfem::Operator | |
int | height |
Dimension of the output / number of rows in the matrix. More... | |
int | width |
Dimension of the input / number of columns in the matrix. More... | |
Definition at line 36 of file maxwell_solver.hpp.
mfem::electromagnetics::MaxwellSolver::MaxwellSolver | ( | ParMesh & | pmesh, |
int | sOrder, | ||
double(*)(const Vector &) | eps, | ||
double(*)(const Vector &) | muInv, | ||
double(*)(const Vector &) | sigma, | ||
void(*)(const Vector &, double, Vector &) | j_src, | ||
Array< int > & | abcs, | ||
Array< int > & | dbcs, | ||
void(*)(const Vector &, double, Vector &) | dEdt_bc | ||
) |
Definition at line 29 of file maxwell_solver.cpp.
mfem::electromagnetics::MaxwellSolver::~MaxwellSolver | ( | ) |
Definition at line 300 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::DisplayToGLVis | ( | ) |
Definition at line 632 of file maxwell_solver.cpp.
|
inline |
Definition at line 70 of file maxwell_solver.hpp.
|
inline |
Definition at line 69 of file maxwell_solver.hpp.
double mfem::electromagnetics::MaxwellSolver::GetEnergy | ( | ) | const |
Definition at line 562 of file maxwell_solver.cpp.
|
inline |
Definition at line 49 of file maxwell_solver.hpp.
double mfem::electromagnetics::MaxwellSolver::GetMaximumTimeStep | ( | ) | const |
Definition at line 507 of file maxwell_solver.cpp.
|
inline |
Definition at line 67 of file maxwell_solver.hpp.
HYPRE_Int mfem::electromagnetics::MaxwellSolver::GetProblemSize | ( | ) |
Definition at line 364 of file maxwell_solver.cpp.
|
virtual |
Solve the equation: k = f(x + dt k, t), for the unknown k at the current time t.
For general F and G, the equation for k becomes: F(x + dt k, k, t) = G(x + dt k, t).
The input vector x corresponds to time index (or cycle) n, while the currently set time, t, and the result vector k correspond to time index n+1. The time step dt corresponds to the time interval between cycles n and n+1.
This method allows for the abstract implementation of some time integration methods, including diagonal implicit Runge-Kutta (DIRK) methods and the backward Euler method in particular.
If not re-implemented, this method simply generates an error.
Reimplemented from mfem::TimeDependentOperator.
Definition at line 404 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::InitializeGLVis | ( | ) |
Definition at line 610 of file maxwell_solver.cpp.
Perform the action of the operator: y = k = f(x, t), where k solves the algebraic equation F(x, k, t) = G(x, t) and t is the current time.
Reimplemented from mfem::TimeDependentOperator.
Definition at line 398 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::PrintSizes | ( | ) |
Definition at line 370 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::RegisterVisItFields | ( | VisItDataCollection & | visit_dc | ) |
Definition at line 575 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::SetInitialBField | ( | VectorCoefficient & | BFieldCoef | ) |
Definition at line 390 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::SetInitialEField | ( | VectorCoefficient & | EFieldCoef | ) |
Definition at line 382 of file maxwell_solver.cpp.
|
inline |
Definition at line 50 of file maxwell_solver.hpp.
void mfem::electromagnetics::MaxwellSolver::SyncGridFuncs | ( | ) |
Definition at line 500 of file maxwell_solver.cpp.
void mfem::electromagnetics::MaxwellSolver::WriteVisItFields | ( | int | it = 0 | ) |
Definition at line 588 of file maxwell_solver.cpp.