ex2.cpp File Reference

Main file for 2nd programming exercise. More...

#include <iostream>
#include <Eigen/Dense>
#include <Eigen/Sparse>
#include <Eigen/SparseCholesky>
#include <adaptiveTriangMesh.h>
#include <configuratorsShellFE.h>
#include <EigenDataContainer.h>
#include <quadratureShellFE.h>
#include <legacyVtkWriter.h>
#include <shellHandler.h>
#include <triangleShellFE.h>
#include "errorEstimator.h"

Go to the source code of this file.

Typedefs
typedef double	RealType
typedef shellFE::ShellElementWithTangentSpaceAtVertex< DataTypeContainerShellFE >	TriangleType
typedef shellFE::AdaptiveTriangMesh< DataTypeContainerShellFE, TriangleType >	MeshType
typedef shellFE::CenterOfMassQuadrature< RealType, typename DataTypeContainerShellFE::DomVecType >	QuadType
typedef shellFE::UnitTriangMeshConfiguratorP1< DataTypeContainerShellFE, MeshType, QuadType >	ConfiguratorType
typedef ConfiguratorType::VectorType	VectorType
typedef ConfiguratorType::SparseMatrixType	SparseMatrixType

Functions
void	printHelp ()
int	main (int argc, char **argv)
	Main function containing the basic structure of the code. More...

Detailed Description

Main file for 2nd programming exercise.

Definition in file ex2.cpp.

Typedef Documentation

typedef shellFE::UnitTriangMeshConfiguratorP1<DataTypeContainerShellFE, MeshType, QuadType> ConfiguratorType

Definition at line 25 of file ex2.cpp.

typedef shellFE::AdaptiveTriangMesh<DataTypeContainerShellFE, TriangleType > MeshType

Definition at line 23 of file ex2.cpp.

typedef shellFE::CenterOfMassQuadrature<RealType, typename DataTypeContainerShellFE::DomVecType> QuadType

Definition at line 24 of file ex2.cpp.

typedef double RealType

Definition at line 21 of file ex2.cpp.

typedef ConfiguratorType::SparseMatrixType SparseMatrixType

Definition at line 27 of file ex2.cpp.

typedef shellFE::ShellElementWithTangentSpaceAtVertex<DataTypeContainerShellFE> TriangleType

Definition at line 22 of file ex2.cpp.

typedef ConfiguratorType::VectorType VectorType

Definition at line 26 of file ex2.cpp.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Main function containing the basic structure of the code.

Definition at line 40 of file ex2.cpp.

{
    unsigned int numRefinements = 10;
    switch (argc) {
    case 3: {
        // First two arguments are handled below
        break;
    }
    case 4: {
        numRefinements = atoi(argv[3]);
        break;
    }
    default: {
        printHelp();
        return 1;
    }
    }

    // Load mesh
    cout << "* Loading grid" << endl;
    MeshType mesh(argv[1]);

    cout << "* Refining " << numRefinements << " times ";
    cout << "  [";
    cout.flush();
    for (unsigned int r = 0; r < numRefinements; ++r) {
        // Estimate error and refine mesh
        ErrorEstimator<ConfiguratorType> errorEstimator(mesh);
        errorEstimator.markAndRefineTriangles();

        // Save the mesh and solution as legacy VTK
        shellFE::LegacyVtkWriter<MeshType> vtkWriter(mesh);
        std::stringstream outputFileName;
        outputFileName << argv[2] << r << ".vtk";
        vtkWriter.save(outputFileName.str().c_str());

        cout << '|';
        cout.flush();
    }
    cout << "]\n\n";

    return 0;
}

void printHelp

(

)

Definition at line 29 of file ex2.cpp.

{
    cout << "Usage: ex2 meshFile resultFile [numRefinements]\n";
    cout << "\tmeshFile\tLegacy VTK file defining the Finite Element mesh "
            << "to use\n";
    cout << "\tresultFile\tLegacy VTK file to write the result to\n";
    cout << "\tnumRefinements\tNumber of times the mesh is refined"
            << " (default: 10)\n\n";
}