WriteWallSurfaceAreaVTKUnitTest.cpp File Reference

#include <Mercury3D.h>
#include <Species/LinearViscoelasticSpecies.h>
#include <Walls/InfiniteWall.h>
#include <errno.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <locale.h>

Classes
class	AreaVTK

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

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

{
    setlocale(LC_ALL, "C");
 
    {
        AreaVTK problem("InfiniteWall");
        problem.solve(argc, argv);
 
        double A0 = problem.getTotalsurfaceAreaFromFile(problem.getName() + "Wall_0.vtu");
        double A1 = problem.getTotalsurfaceAreaFromFile(problem.getName() + "Wall_1.vtu");
 
        double checkedArea = 4.0;   // calculated from domain size
        helpers::check(A0, A1, 0, "Surface areas of first object are equal");
        helpers::check(A0, checkedArea, 1e-3, "Surface area of first object");
    }
 
    {
        AreaVTK problem("../SelfTests/Walls/Casing.stl");
        problem.solve(argc, argv);
 
        double A0 = problem.getTotalsurfaceAreaFromFile(problem.getName() + "Wall_0.vtu");
        double A1 = problem.getTotalsurfaceAreaFromFile(problem.getName() + "Wall_1.vtu");
 
        double checkedArea = 675075.0;  // checked via MeshLab
        helpers::check(A0, A1, 0, "Surface areas of second object are equal");
        helpers::check(A0, checkedArea, 1e-1, "Surface area of second object");
    }
 
    return 0;
}

References helpers::check(), DPMBase::getName(), AreaVTK::getTotalsurfaceAreaFromFile(), and DPMBase::solve().