example_15.cpp File Reference

Construction of Laplace operator on a uniform hexahedral mesh using arbitrary-degree elements. This assembles the matrix into STL data structures vector<map<int,double> > to get logarithmic access to columns. More...

#include "Intrepid_FunctionSpaceTools.hpp"
#include "Intrepid_FieldContainer.hpp"
#include "Intrepid_CellTools.hpp"
#include "Intrepid_HGRAD_HEX_Cn_FEM.hpp"
#include "Intrepid_DefaultCubatureFactory.hpp"
#include "Intrepid_Utils.hpp"
#include "Epetra_Time.h"
#include "Epetra_Map.h"
#include "Epetra_FEVector.h"
#include "Epetra_FECrsMatrix.h"
#include "Epetra_SerialComm.h"
#include "Teuchos_oblackholestream.hpp"
#include "Teuchos_RCP.hpp"
#include "Shards_CellTopology.hpp"
#include "EpetraExt_MultiVectorOut.h"
#include <vector>
#include <map>

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Detailed Description

Construction of Laplace operator on a uniform hexahedral mesh using arbitrary-degree elements. This assembles the matrix into STL data structures vector<map<int,double> > to get logarithmic access to columns.

div grad u = f in Omega
u = 0 on Gamma 

Discrete linear system for nodal coefficients(x):
    
Kx = b

K - HGrad stiffness matrix
b - right hand side vector 

Author

Created by P. Bochev, R. Kirby, D. Ridzal and K. Peterson.

Remarks

Usage

./Intrepid_example_Drivers_Example_15.exe N verbose
int degree          - polynomial degree
int NX              - num intervals in x direction (assumed box domain, 0,1)
int NY              - num intervals in x direction (assumed box domain, 0,1)
int NZ              - num intervals in x direction (assumed box domain, 0,1)
verbose (optional)  - any character, indicates verbose output

Sample command line

./Intrepid_example_Drivers_Example_15.exe 2 10 10 10 

Definition in file example_15.cpp.