Collaboration diagram for ToolContext:

Public Member Functions
	ToolContext (moab::Interface *icore)
	~ToolContext ()
void	timer_push (std::string operation)
void	timer_pop ()
void	ParseCLOptions (int argc, char *argv[])
Public Attributes
moab::Interface *	mbcore
const int	proc_id
const int	n_procs
moab::DebugOutput	outputFormatter
int	blockSize
std::vector< std::string >	inFilenames
std::vector< Mesh * >	meshes
std::vector< moab::EntityHandle >	meshsets
std::vector< int >	disc_orders
std::vector< std::string >	disc_methods
std::vector< std::string >	doftag_names
std::string	outFilename
std::string	intxFilename
std::string	baselineFile
moab::TempestRemapper::TempestMeshType	meshType
bool	computeDual
bool	computeWeights
bool	verifyWeights
bool	enforceConvexity
int	ensureMonotonicity
bool	fNoConservation
bool	fVolumetric
bool	rrmGrids
bool	kdtreeSearch
bool	fNoBubble
bool	fInputConcave
bool	fOutputConcave
bool	fCheck
Private Attributes
moab::CpuTimer *	timer
double	timer_ops
std::string	opName

Detailed Description

Definition at line 42 of file mbtempest.cpp.

Constructor & Destructor Documentation

ToolContext::ToolContext ( moab::Interface * icore ) [inline]

Definition at line 77 of file mbtempest.cpp.

References doftag_names, inFilenames, outputFormatter, moab::DebugOutput::set_prefix(), and timer.

        : mbcore( icore ), proc_id( 0 ), n_procs( 1 ), outputFormatter( std::cout, 0, 0 ),
#endif
          blockSize( 5 ), outFilename( "output.exo" ), intxFilename( "" ), baselineFile( "" ),
          meshType( moab::TempestRemapper::DEFAULT ), computeDual( false ), computeWeights( false ),
          verifyWeights( false ), enforceConvexity( false ), ensureMonotonicity( 0 ), fNoConservation( false ),
          fVolumetric( false ), rrmGrids( false ), kdtreeSearch( true ), fNoBubble( true ), fInputConcave( false ),
          fOutputConcave( false ), fCheck( n_procs > 1 ? false : true )
    {
        inFilenames.resize( 2 );
        doftag_names.resize( 2 );
        timer = new moab::CpuTimer();

        outputFormatter.set_prefix( "[MBTempest]: " );
    }

ToolContext::~ToolContext ( ) [inline]

Definition at line 93 of file mbtempest.cpp.

References baselineFile, disc_methods, disc_orders, doftag_names, inFilenames, intxFilename, meshes, meshsets, outFilename, and timer.

    {
        // for (unsigned i=0; i < meshes.size(); ++i) delete meshes[i];
        meshes.clear();
        inFilenames.clear();
        disc_orders.clear();
        disc_methods.clear();
        doftag_names.clear();
        outFilename.clear();
        intxFilename.clear();
        baselineFile.clear();
        meshsets.clear();
        delete timer;
    }

Member Function Documentation

void ToolContext::ParseCLOptions	(	int	argc,
		char *	argv[]
	)		`[inline]`

Definition at line 135 of file mbtempest.cpp.

References ProgOptions::addOpt(), baselineFile, blockSize, computeDual, computeWeights, moab::TempestRemapper::CS, moab::TempestRemapper::DEFAULT, disc_methods, disc_orders, doftag_names, enforceConvexity, ensureMonotonicity, fCheck, fNoBubble, fNoConservation, fVolumetric, ProgOptions::getOptAllArgs(), moab::TempestRemapper::ICO, inFilenames, intxFilename, kdtreeSearch, meshType, ProgOptions::numOptSet(), outFilename, moab::TempestRemapper::OVERLAP_FILES, moab::TempestRemapper::OVERLAP_MEMORY, moab::TempestRemapper::OVERLAP_MOAB, ProgOptions::parseCommandLine(), proc_id, moab::TempestRemapper::RLL, rrmGrids, and verifyWeights.

Referenced by main().

    {
        ProgOptions opts;
        int imeshType                  = 0;
        std::string expectedFName      = "output.exo";
        std::string expectedMethod     = "fv";
        std::string expectedDofTagName = "GLOBAL_ID";
        int expectedOrder              = 1;
        bool fBubble                   = false;

        if( !proc_id )
        {
            std::cout << "Command line options provided to mbtempest:\n  ";
            for( int iarg = 0; iarg < argc; ++iarg )
                std::cout << argv[iarg] << " ";
            std::cout << std::endl << std::endl;
        }

        opts.addOpt< int >( "type,t",
                            "Type of mesh (default=CS; Choose from [CS=0, RLL=1, ICO=2, OVERLAP_FILES=3, "
                            "OVERLAP_MEMORY=4, OVERLAP_MOAB=5])",
                            &imeshType );
        opts.addOpt< int >( "res,r", "Resolution of the mesh (default=5)", &blockSize );
        opts.addOpt< void >( "dual,d", "Output the dual of the mesh (relevant only for ICO mesh type)", &computeDual );
        opts.addOpt< std::string >( "file,f", "Output computed mesh or remapping weights to specified filename",
                                    &outFilename );
        opts.addOpt< std::string >(
            "load,l", "Input mesh filenames for source and target meshes. (relevant only when computing weights)",
            &expectedFName );
        opts.addOpt< void >( "noconserve,c",
                             "Do not apply conservation to the resultant weights (relevant only "
                             "when computing weights)",
                             &fNoConservation );
        opts.addOpt< void >( "volumetric,v",
                             "Apply a volumetric projection to compute the weights (relevant only "
                             "when computing weights)",
                             &fVolumetric );
        opts.addOpt< void >( "advfront,a",
                             "Use the advancing front intersection instead of the Kd-tree based algorithm "
                             "to compute mesh intersections. (relevant only when computing weights)" );
        opts.addOpt< int >( "monotonic,n", "Ensure monotonicity in the weight generation", &ensureMonotonicity );
        opts.addOpt< int >( "order,o", "Discretization orders for the source and target solution fields",
                            &expectedOrder );
        opts.addOpt< std::string >( "method,m", "Discretization method for the source and target solution fields",
                                    &expectedMethod );
        opts.addOpt< std::string >( "global_id,g",
                                    "Tag name that contains the global DoF IDs for source and target solution fields",
                                    &expectedDofTagName );

        opts.addOpt< void >( "weights,w",
                             "Compute and output the weights using the overlap mesh (generally "
                             "relevant only for OVERLAP mesh)",
                             &computeWeights );
        opts.addOpt< std::string >( "intx,i", "Output TempestRemap intersection mesh filename", &intxFilename );

        opts.addOpt< std::string >( "baseline", "Output baseline file", &baselineFile );

        opts.addOpt< void >( "rrmgrids",
                             "At least one of the meshes is a regionally refined grid (relevant to "
                             "accelerate intersection computation)",
                             &rrmGrids );
        opts.addOpt< void >( "checkmap", "Check the generated map for conservation and consistency", &fCheck );
        opts.addOpt< void >( "verify",
                             "Verify the accuracy of the maps by projecting analytical functions "
                             "from source to target "
                             "grid by applying the maps",
                             &verifyWeights );
        opts.addOpt< void >( "enforce_convexity", "check convexity of input meshes to compute mesh intersections",
                             &enforceConvexity );
        opts.addOpt< void >( "bubble", "use bubble on interior of spectral element nodes", &fBubble );

        opts.parseCommandLine( argc, argv );

        // By default - use Kd-tree based search; if user asks for advancing front, disable Kd-tree
        // algorithm
        kdtreeSearch = opts.numOptSet( "advfront,a" ) == 0;

        // negate the option
        fNoBubble = !fBubble;

        switch( imeshType )
        {
            case 0:
                meshType = moab::TempestRemapper::CS;
                break;

            case 1:
                meshType = moab::TempestRemapper::RLL;
                break;

            case 2:
                meshType = moab::TempestRemapper::ICO;
                break;

            case 3:
                meshType = moab::TempestRemapper::OVERLAP_FILES;
                break;

            case 4:
                meshType = moab::TempestRemapper::OVERLAP_MEMORY;
                break;

            case 5:
                meshType = moab::TempestRemapper::OVERLAP_MOAB;
                break;

            default:
                meshType = moab::TempestRemapper::DEFAULT;
                break;
        }

        if( meshType > moab::TempestRemapper::ICO )
        {
            opts.getOptAllArgs( "load,l", inFilenames );
            opts.getOptAllArgs( "order,o", disc_orders );
            opts.getOptAllArgs( "method,m", disc_methods );
            opts.getOptAllArgs( "global_id,i", doftag_names );

            if( disc_orders.size() == 0 )
            {
                disc_orders.resize( 2, 1 );
            }

            if( disc_orders.size() == 1 )
            {
                disc_orders.push_back( 1 );
            }

            if( disc_methods.size() == 0 )
            {
                disc_methods.resize( 2, "fv" );
            }

            if( disc_methods.size() == 1 )
            {
                disc_methods.push_back( "fv" );
            }

            if( doftag_names.size() == 0 )
            {
                doftag_names.resize( 2, "GLOBAL_ID" );
            }

            if( doftag_names.size() == 1 )
            {
                doftag_names.push_back( "GLOBAL_ID" );
            }

            assert( inFilenames.size() == 2 );
            assert( disc_orders.size() == 2 );
            assert( disc_methods.size() == 2 );
        }

        expectedFName.clear();
    }

void ToolContext::timer_pop ( ) [inline]

Definition at line 114 of file mbtempest.cpp.

References n_procs, opName, proc_id, moab::CpuTimer::time_since_birth(), timer, and timer_ops.

Referenced by main().

    {
        double locElapsed = timer->time_since_birth() - timer_ops, avgElapsed = 0, maxElapsed = 0;
#ifdef MOAB_HAVE_MPI
        MPI_Reduce( &locElapsed, &maxElapsed, 1, MPI_DOUBLE, MPI_MAX, 0, pcomm->comm() );
        MPI_Reduce( &locElapsed, &avgElapsed, 1, MPI_DOUBLE, MPI_SUM, 0, pcomm->comm() );
#else
        maxElapsed = locElapsed;
        avgElapsed = locElapsed;
#endif
        if( !proc_id )
        {
            avgElapsed /= n_procs;
            std::cout << "[LOG] Time taken to " << opName.c_str() << ": max = " << maxElapsed
                      << ", avg = " << avgElapsed << "\n";
        }
        // std::cout << "\n[LOG" << proc_id << "] Time taken to " << opName << " = " <<
        // timer->time_since_birth() - timer_ops << std::endl;
        opName.clear();
    }

void ToolContext::timer_push ( std::string operation ) [inline]

Definition at line 108 of file mbtempest.cpp.

References opName, moab::CpuTimer::time_since_birth(), timer, and timer_ops.

Referenced by main().

    {
        timer_ops = timer->time_since_birth();
        opName    = operation;
    }

Member Data Documentation

std::string ToolContext::baselineFile

Definition at line 59 of file mbtempest.cpp.

Referenced by main(), ParseCLOptions(), and ~ToolContext().

int ToolContext::blockSize

Definition at line 50 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), and ParseCLOptions().

bool ToolContext::computeDual

Definition at line 61 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), and ParseCLOptions().

bool ToolContext::computeWeights

Definition at line 62 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

std::vector< std::string > ToolContext::disc_methods

Definition at line 55 of file mbtempest.cpp.

Referenced by main(), ParseCLOptions(), and ~ToolContext().

std::vector< int > ToolContext::disc_orders

Definition at line 54 of file mbtempest.cpp.

Referenced by main(), ParseCLOptions(), and ~ToolContext().

std::vector< std::string > ToolContext::doftag_names

Definition at line 56 of file mbtempest.cpp.

Referenced by main(), ParseCLOptions(), ToolContext(), and ~ToolContext().

bool ToolContext::enforceConvexity

Definition at line 64 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

int ToolContext::ensureMonotonicity

Definition at line 65 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

bool ToolContext::fCheck

Definition at line 70 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

bool ToolContext::fInputConcave

Definition at line 70 of file mbtempest.cpp.

Referenced by main().

bool ToolContext::fNoBubble

Definition at line 70 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

bool ToolContext::fNoConservation

Definition at line 66 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

bool ToolContext::fOutputConcave

Definition at line 70 of file mbtempest.cpp.

Referenced by main().

bool ToolContext::fVolumetric

Definition at line 67 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

std::vector< std::string > ToolContext::inFilenames

Definition at line 51 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), ParseCLOptions(), ToolContext(), and ~ToolContext().

std::string ToolContext::intxFilename

Definition at line 58 of file mbtempest.cpp.

Referenced by main(), ParseCLOptions(), and ~ToolContext().

bool ToolContext::kdtreeSearch

Definition at line 69 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

moab::Interface* ToolContext::mbcore

Definition at line 44 of file mbtempest.cpp.

Referenced by CreateTempestMesh().

std::vector< Mesh* > ToolContext::meshes

Definition at line 52 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), and ~ToolContext().

std::vector< moab::EntityHandle > ToolContext::meshsets

Definition at line 53 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), and ~ToolContext().

moab::TempestRemapper::TempestMeshType ToolContext::meshType

Definition at line 60 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), and ParseCLOptions().

const int ToolContext::n_procs

Definition at line 48 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), and timer_pop().

std::string ToolContext::opName [private]

Definition at line 294 of file mbtempest.cpp.

Referenced by timer_pop(), and timer_push().

std::string ToolContext::outFilename

Definition at line 57 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), ParseCLOptions(), and ~ToolContext().

moab::DebugOutput ToolContext::outputFormatter

Definition at line 49 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), and ToolContext().

const int ToolContext::proc_id

Definition at line 48 of file mbtempest.cpp.

Referenced by CreateTempestMesh(), main(), ParseCLOptions(), and timer_pop().

bool ToolContext::rrmGrids

Definition at line 68 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

moab::CpuTimer* ToolContext::timer [private]

Definition at line 292 of file mbtempest.cpp.

Referenced by timer_pop(), timer_push(), ToolContext(), and ~ToolContext().

double ToolContext::timer_ops [private]

Definition at line 293 of file mbtempest.cpp.

Referenced by timer_pop(), and timer_push().

bool ToolContext::verifyWeights

Definition at line 63 of file mbtempest.cpp.

Referenced by main(), and ParseCLOptions().

List of all members.

The documentation for this struct was generated from the following file:

mbtempest.cpp

Public Member Functions

Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation