makeflo.cpp File Reference

#include <stdio.h>
#include "makeflo.h"

Include dependency graph for makeflo.cpp:

Go to the source code of this file.

Functions
void	checkError (const char *errCode)
void	printUsage (const char *why)
int	main (int argc, char *argv[])

Variables
makefloParam	parameters

Function Documentation

void checkError

(

const char *

errCode

)

Definition at line 62 of file makeflo.cpp.

Referenced by main().

                                     {
  if (errCode==NULL) return; //No problem
  fprintf(stderr,"FATAL ERROR! %s\n",errCode);
  exit(1);
}

Here is the caller graph for this function:

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 103 of file makeflo.cpp.

References blocks, buildFaces(), checkError(), checkQuality(), endsWith(), i, printUsage(), blockConsumer::read(), readBoundaries(), replaceExtention(), makefloParam::setLevel(), makefloParam::splitAxis, splitBlocks(), makefloParam::splithalf, makefloParam::splitRCB, writeBlocks(), writeFlo(), and writeTop().

{
//Parse the command line arguments
        if (argc<5) printUsage("Not enough command-line arguments");
        int curArg=1;
        while (argv[curArg][0]=='-') 
        {//Parse a flag argument
                const char *flag=argv[curArg++];
                if (0==strcmp(flag,"-h") ||
                    0==strcmp(flag,"-?") ||
                    0==strcmp(flag,"--help")) 
                        printUsage(NULL);
                else if (0==strcmp(flag,"-splitaxis")) {
                        if (1!=sscanf(argv[curArg++],"%d",&parameters.splitAxis))
                                printUsage("Must pass 0,1, or 2 to -splitaxis");
                }
                else if (0==strcmp(flag,"-splitrcb")) {
                        parameters.splitRCB=1;
                }
                else if (0==strcmp(flag,"-splithalf")) {
                        parameters.splithalf=1;
                }
                else if (0==strcmp(flag,"-multilevel")) {
                        int level=0;
                        if ((1!=sscanf(argv[curArg++],"%d",&level))||(level<=0))
                                printUsage("Must pass a positive integer to -multilevel");
                        parameters.setLevel(level);
                }               
                else printUsage("Unrecognized flag");
        }
        if (curArg+4!=argc) printUsage("Too many arguments");
        const char *inMesh=argv[curArg+0];
        string inpFile=replaceExtention(argv[curArg+0],".inp");
        int nPieces=atoi(argv[curArg+1]);
        const char *floFile=argv[curArg+2];
        const char *outMesh=argv[curArg+3];

        bool forFlo;
        const char *bcExt="?";
        if (endsWith(floFile,".flo")) {
                forFlo=true;
                bcExt=".bc";
        }
        else if (endsWith(floFile,".top")) {
                forFlo=false;
                bcExt=".bcmp";
        }
        else {
                printf("UNRECOGNIZED TOPOLOGY NAME: '%s'\n",floFile);
                exit(1);
        }
        string bcFile=replaceExtention(argv[curArg+0],bcExt);
        
        vector<block *> blocks;
        
//Read the blocks of the mesh
        { //<- scoping block so r's destructor gets called
                blockReader r(blocks);
                checkError(r.read(inMesh));
                printf("Mesh file read successfully.\n");
        }

// Check the quality of the mesh
        for (unsigned int i=0;i<blocks.size();i++)
                checkQuality(*blocks[i]);
        
//Read the boundary conditions
        checkError(readBoundaries(blocks,inpFile.c_str()));
        printf("Boundary conditions read successfully.\n");
        
        if (parameters.splithalf) {
                // split in half for every block
                nPieces = blocks.size()*2;
                printf("Split every block in half for %d pieces.\n", nPieces);
        }

//Split the blocks (for greater parallelism)
        checkError(splitBlocks(blocks,nPieces));
        printf("Split blocks for %d PEs (%d blocks generated).\n",
                nPieces,(int)blocks.size());    

//Build the block faces & associate all shared nodes
        buildFaces(blocks,forFlo);
        printf("Nodes matched successfully.\n");
        
//Write out the block patches to the topology file
        if (forFlo) { 
                checkError(writeFlo(blocks,nPieces,bcFile.c_str(),floFile));
                printf(".flo file written successfully\n");
        } else /*forTop*/ {
                checkError(writeTop(blocks,bcFile.c_str(),floFile));
                printf(".top file written successfully\n");
        }

//Write out the blocks themselves
        checkError(writeBlocks(blocks,outMesh));
        printf("Program finished successfully\n");
        return 0;
}

Here is the call graph for this function:

void printUsage

(

const char *

why

)

Definition at line 68 of file makeflo.cpp.

References MAKEFLO_VERSION.

Referenced by ComSwitch::initializeFileStreams(), main(), and ComSwitch::parseOptions().

{
        printf("Usage: makeflo [ -splitaxis <-1, 0, 1, or 2)> ]\n"
          "      [ -splitrcb ]      [ -splithalf ]\n"
          "      [ -multilevel <maximum level> ]\n"
          "     <.grd, .grdb, .msh, or .hdf; and .inp and .bc file> \n"
          "     <# PEs> <.flo or .top output file> <.msh or .hdf output file>\n"
          "\n"
          "  Makeflo reads a mesh (first parameter) and boundary condition\n"
          "(.inp) list and description (.bc), partitions the mesh into the\n"
          "requested number of processors, and writes out a rocflo control\n"
          "file (.flo or .top) and the partitioned mesh (last parameter).\n"
          "  The mesh formats supported are:\n"
          "<*.grd>  A Gridgen ASCII double-precision mesh description\n"
          "<*.msh>  A double-precision binary Rocflo mesh file\n"
          "<*.hdf>  A double-precision 3D HDF mesh file\n"
          "  When the mesh consists of multiple files, give only the name of\n"
          "the first file (e.g., 'tstflo_001.hdf').  The number of blocks and\n"
          "numeric format will be automatically determined.\n"
          "\n"
          "Options:\n"
          " -splitaxis <-1, 0, 1, or 2)>: axis to split along. \n"
          " -splitrcb:   split using recursive bisection.\n"
          " -multilevel <maximum level>: number of multigrid levels to create \n"
          " -splithalf:  split every block in half, <# PEs> is ignored.\n"
          "\n"
          "Part of the CSAR Tools.  Version " MAKEFLO_VERSION "\n"
        ); 
        if (why!=NULL)
                printf("Exiting> %s\n",why);
        exit(1);
}

Here is the caller graph for this function:

Variable Documentation

makefloParam parameters

Definition at line 101 of file makeflo.cpp.

Referenced by checkSpan(), blockReader::consume(), face::face(), block::split(), splitBlocks(), writeBlock(), and writeMblock().