http://trilinos.sandia.gov/packages/docs/r10.4/packages/intrepid/example/Drivers/example_10.cpp File Reference

Construction of Laplace operator on a uniform hexahedral mesh using arbitrary-degree elements. This is the most naive implementation wherein we form the stiffness matrix on each cell by quadrature and do not preallocate the global matrix graph before assembling. 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"`

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 is the most naive implementation wherein we form the stiffness matrix on each cell by quadrature and do not preallocate the global matrix graph before assembling.

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

Discrete linear system for nodal coefficients(x):

Kx = b

b - right hand side vector

```
Remarks:
Usage
```
./Intrepid_example_Drivers_Example_10.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_10.exe 2 10 10 10
```

Definition in file example_10.cpp.

Generated on Mon Jan 31 09:55:49 2011 for Intrepid by  1.6.3