CGSelectRegionSelfTest.cpp File Reference

#include <Walls/InfiniteWall.h>
#include "Mercury3D.h"
#include "Species/LinearViscoelasticSpecies.h"

Classes
class	Chain

Functions
int	main ()

Function Documentation

main()

int main

(

)

A simple DPM problem is coarse-grained in a select region to show selective CG works

Returns

{
    logger(INFO,"\nRun a simple DPM problem\n");
    Chain chain;
    chain.solve();
 
    logger(INFO,"\nCoarse-grain a select region (1<z<2)\n");
    std::string cmd = "./MercuryCG Chain -tMin 1 -timeAverage -z 1 2 -averageBeyondDomain 0";
    if (system(cmd.c_str())==-1) {
        logger(WARN,"system call failed");
    }
 
    logger(INFO,"\nRead in CG\n");
    std::ifstream statFile("Chain.stat");
    //ignore 2 header lines and coordinate
    statFile.ignore(2000,'\n');
    statFile.ignore(2000,'\n');
    statFile.ignore(2000,' ');
    statFile.ignore(2000,' ');
    //read in volume fraction
    Mdouble density = 0;
    statFile >> density;
    logger(INFO,"Density % (should be 1)",density);
}

References UniformPSDSelfTest::density, INFO, logger, DPMBase::solve(), oomph::Global_string_for_annotation::string(), and WARN.