main.h File Reference

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <ostream>
#include <fstream>
#include <string>
#include <time.h>
#include <getopt.h>
#include <vector>
#include <limits>
#include "file.h"
#include "fileselect.h"
#include "metric.h"
#include "interpolation.h"
#include "datatypedef.h"
#include "util.h"
#include "BoundingBox.hpp"

Include dependency graph for main.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	ComSwitch

Functions
string	gUsage ("Usage: [--outfile | -o] <output file prefix>\n\t""[--infile1 | -1] <intput file 1 name>:<dimensions>\n\t""[--infile2 | -2] <intput file 2 name>:<dimensions>\n\t""[--fieldmap | -f] <File (1,2)>:<File Dimension>/<Mapped Dimension or \"null\">\n\t""[--index-order | -i] <File (1,2)>:<Partition>:<Dimension 1>,<Dimension 2>,...\n\t""[--conv | -c] <File (1,2)>:<Variable>:<Conversion Factor> \n\t""[--norm | -n] <File (1,2)>:<Variable>:<Normalization Factor> \n\t""[--range | -r] <xmin>,<ymin>,<zmin>,<xmax>,<ymax>,<zmax> \n\t""[--metric | -m] <metric number>=<file 1 partition 1>,<file 1 partition 2>,...:<file 1 variable>/\n\t"" <file 2 partition>:<file 2 variable>\n\t""[--verbose | -v]\n")
void	printUsage (string rErrorMessage)
	Print the rc_cmp usage message. More...
vector< datafile * >	openAndParseFiles ()
	Open the files specified at the command line and parse them. More...
void	compareFiles (vector< datafile * > pFiles)
	Compare the two files, interpolating the points from the second onto the first mesh. More...
std::vector< points * >	retrievePoints (datafile *pFile, vector< int > rPartitions)
	Construct a vector of pointers to datasets (point objects) from a list of partition numbers. More...
void	printPointValues (points *pPoints)
	Write point values to standard out. More...
long double	calcAverageValue (points *pPoints, int rVar)
	Calculates the average value of the points. More...
long double	minPoint (points *pPoints, int rVar)
	Find the minimum value. More...
long double	maxPoint (points *pPoints, int rVar)
	Find the maximum value. More...

Function Documentation

long double calcAverageValue	(	points *	pPoints,
		int	rVar
	)

Calculates the average value of the points.

Parameters

pPoints	Datapoints
rVar	Variable number to average

Returns

Returns the average

Definition at line 276 of file main.cpp.

References points::get_num_points(), points::get_point(), and pnt::vals.

                                                       {

    long double sum = 0.0;

    //Sum the points
    for(int point_num = 0; point_num < pPoints->get_num_points(); point_num++){
        pnt curr = pPoints->get_point(point_num);
        sum += curr.vals[rVar];
    }

    return (sum / (long double)pPoints->get_num_points());
}

Here is the call graph for this function:

void compareFiles

(

vector< datafile * >

pFiles

)

Compare the two files, interpolating the points from the second onto the first mesh.

Then, ouput results.

Definition at line 95 of file main.cpp.

References Interpolate::bilinear(), points::contains_null(), points::delete_all(), cmp_map::file1_partitions, cmp_map::file2_partition, metric::get_metric_name(), datafile::get_points(), metric::get_res(), ComSwitch::getComparisonList(), ComSwitch::getLogFile(), ComSwitch::getOutPrefix(), ComSwitch::getRange(), ComSwitch::getVerbosity(), i, ComSwitch::isRangeSpecified(), itoa(), cmp_map::metric, metric_select(), retrievePoints(), same_elements(), metric::setRange(), cmp_map::var1, and cmp_map::var2.

Referenced by main().

                                           {

    datafile *file1 = pFiles[0];
    datafile *file2 = pFiles[1];

    if(gCmdSwitches.getVerbosity())
        gCmdSwitches.getLogFile() << "\n";


    //Get comparison list
    std::vector<cmp_map> comparison_list = gCmdSwitches.getComparisonList();

    std::vector<points*> interp_pts;

    for(int i=0; i<comparison_list.size(); i++){

        cmp_map curr_comp = comparison_list[i];

        //Check if previous interp_pts satisfies current interp_pts
        if( i > 0 ){
            cmp_map prev_comp = comparison_list[(i-1)];

            if( !same_elements(curr_comp.file1_partitions, curr_comp.file1_partitions) ||
                curr_comp.file2_partition != prev_comp.file2_partition                 ||
                curr_comp.var1 != prev_comp.var1                                       ||
                curr_comp.var2 != prev_comp.var2                                       ){

                points::delete_all(interp_pts);

                for(int partition_1_num = 0; partition_1_num < curr_comp.file1_partitions.size(); partition_1_num++){
                    points *temp_points;
                    temp_points = Interpolate::bilinear( file1->get_points(curr_comp.file1_partitions[partition_1_num]),
                                                         file2->get_points(curr_comp.file2_partition),
                                                         curr_comp.var1,
                                                         curr_comp.var2 );
                    interp_pts.push_back(temp_points);
                }

            }
        }
        else{

            for(int partition_1_num = 0; partition_1_num < curr_comp.file1_partitions.size(); partition_1_num++){
                points *temp_points;
                temp_points = Interpolate::bilinear( file1->get_points(curr_comp.file1_partitions[partition_1_num]),
                                                     file2->get_points(curr_comp.file2_partition),
                                                     curr_comp.var1,
                                                     curr_comp.var2 );
                interp_pts.push_back(temp_points);
            }
        }

        if(points::contains_null(interp_pts)){

            string out_msg;
            out_msg  = "ERROR   |  Interpolation Failed, ";
            out_msg += "Possible Invalid Partition or Variable Argument: \"";
            out_msg += itoa(curr_comp.metric) + "=";
            
            for(int partition_1_num = 0; partition_1_num < curr_comp.file1_partitions.size(); partition_1_num++){
                
                if(partition_1_num != 0)
                    out_msg += ",";
                
                out_msg += itoa(curr_comp.file1_partitions[partition_1_num]);
            }
            
            out_msg += ":";
            out_msg += itoa(curr_comp.var1) + "/" + itoa(curr_comp.file2_partition) + ":";
            out_msg += itoa(curr_comp.var2) + "\"";

            gCmdSwitches.getLogFile() << out_msg << endl;
            cout << out_msg << endl;

            continue;
        }

        //Construct the outprefix string for the metric
        //Used by the metric to prefix output files
        string metric_outprefix = gCmdSwitches.getOutPrefix();
        metric_outprefix += "_part1_";

        for(int partition_1_num = 0; partition_1_num < curr_comp.file1_partitions.size(); partition_1_num++){
                
            if(partition_1_num != 0)
                metric_outprefix += "_";
                
            metric_outprefix += itoa(curr_comp.file1_partitions[partition_1_num]);
        }

        metric_outprefix += "_part2_" + itoa(curr_comp.file2_partition);

        //Create metric object
        metric *comp = metric_select(curr_comp.metric, retrievePoints(file1, curr_comp.file1_partitions),
                                     interp_pts, curr_comp.var1, 2, //Interpolation only returns the desired point
                                     gCmdSwitches.getLogFile(), metric_outprefix);

        if( gCmdSwitches.isRangeSpecified( ) )
        {
                BoundingBox range = gCmdSwitches.getRange( );
                comp ->setRange( &range );
        }

        //Check metric object
        if(comp == NULL){
            gCmdSwitches.getLogFile() << "ERROR   |  " << "Invalid Metric Number: \"";
            gCmdSwitches.getLogFile() <<  curr_comp.metric << "\"" << endl;

            cout << "ERROR   |  " << "Invalid Metric Number: \"";
            cout <<  curr_comp.metric << "\"" << endl;

            continue;
        }

        //Calculate metric
        string res = comp->get_res();

        //Print result
        gCmdSwitches.getLogFile().precision(10);
        gCmdSwitches.getLogFile().width(65);
        string msg_out = "RESULT  |  (" + itoa(curr_comp.file1_partitions[0]) + ":"
                                         + itoa(curr_comp.var1) + " "
                                         + itoa(curr_comp.file2_partition) + ":"
                                         + itoa(curr_comp.var2) + ") "
                                         + comp->get_metric_name() + ": ";

        gCmdSwitches.getLogFile() << std::left << msg_out;
        gCmdSwitches.getLogFile() << std::right << res << "\n";

        cout.precision(10);
        cout.width(65);
        cout << std::left << msg_out;
        cout << std::right << res << "\n";

        delete comp;
    }

    //Cleanup Memory
    points::delete_all(interp_pts);
}

Here is the call graph for this function:

Here is the caller graph for this function:

string gUsage	(	"Usage: <output file prefix>\n\t"" <intput file 1 name>:<dimensions>\n\t"" <intput file 2 name>:<dimensions>\n\t"" <File (1,2)>:<File Dimension>/<Mapped Dimension or \"null\">\n\t"" <File (1,2)>:<Partition>:<Dimension 1>	[--outfile|-o][--infile1|-1][--infile2|-2][--fieldmap|-f][--index-order|-i],
		< Dimension 2 >	,
		...\n\t""< File(1, 2)>:< Variable >:< Conversion Factor >\n\t""< File(1, 2)>:< Variable >:< Normalization Factor >\n\t""< xmin >	[--conv|-c][--norm|-n][--range|-r],
		< ymin >	,
		< zmin >	,
		< xmax >	,
		< ymax >	,
		< zmax >\n\t""< metric number >	[--metric|-m] = <file 1 partition 1>,
		< file 1 partition 2 >	,
		...:< file 1 variable >/\n\t""< file 2 partition >:< file 2 variable >\n\t""\n"	[--verbose|-v]
	)

Referenced by printUsage().

Here is the caller graph for this function:

long double maxPoint	(	points *	pPoints,
		int	rVar
	)

Find the maximum value.

Parameters

pPoints	Points
rVar	Variable number

Definition at line 309 of file main.cpp.

References points::get_num_points(), points::get_point(), max(), and pnt::vals.

                                               {
    long double max = -std::numeric_limits<long double>::infinity();

    for(int point_num = 0; point_num < pPoints->get_num_points(); point_num++){
        pnt curr = pPoints->get_point(point_num);
        if(curr.vals[rVar] > max)
            max = curr.vals[rVar];
    }

    return max;
}

Here is the call graph for this function:

long double minPoint	(	points *	pPoints,
		int	rVar
	)

Find the minimum value.

Parameters

pPoints	Points
rVar	Variable number

Definition at line 293 of file main.cpp.

References points::get_num_points(), points::get_point(), min(), and pnt::vals.

                                               {
    long double min = std::numeric_limits<long double>::infinity();

    for(int point_num = 0; point_num < pPoints->get_num_points(); point_num++){
        pnt curr = pPoints->get_point(point_num);
        if(curr.vals[rVar] < min)
            min = curr.vals[rVar];
    }

    return min;
}

Here is the call graph for this function:

vector<datafile*> openAndParseFiles

(

)

Open the files specified at the command line and parse them.

Returns

Pointers to the files

Definition at line 40 of file main.cpp.

References Fileselect::detectFiletype(), ComSwitch::getConvFactor(), ComSwitch::getFieldMappings(), ComSwitch::getFileDimensions(), ComSwitch::getFileName(), ComSwitch::getIndexOrder(), ComSwitch::getInputFile(), ComSwitch::getLogFile(), ComSwitch::getNormVal(), ComSwitch::getVerbosity(), and datafile::parse().

Referenced by main().

                                     {

    vector<datafile*> file_pointers;

    //Create file objects
    datafile *file1 = Fileselect::detectFiletype(gCmdSwitches.getInputFile(1),
                                                 gCmdSwitches.getLogFile(),
                                                 gCmdSwitches.getVerbosity(),
                                                 gCmdSwitches.getFileName(1),
                                                 gCmdSwitches.getFileDimensions(1),
                                                 gCmdSwitches.getFieldMappings(1),
                                                 gCmdSwitches.getIndexOrder(1),
                                                 gCmdSwitches.getConvFactor(1),
                                                 gCmdSwitches.getNormVal(1));

    datafile *file2 = Fileselect::detectFiletype(gCmdSwitches.getInputFile(2),
                                                 gCmdSwitches.getLogFile(),
                                                 gCmdSwitches.getVerbosity(),
                                                 gCmdSwitches.getFileName(2),
                                                 gCmdSwitches.getFileDimensions(2),
                                                 gCmdSwitches.getFieldMappings(2),
                                                 gCmdSwitches.getIndexOrder(2),
                                                 gCmdSwitches.getConvFactor(2),
                                                 gCmdSwitches.getNormVal(2));

    //Check to make sure files have been properly opened
    if(file1 == NULL || file2 == NULL){

        if(file1 != NULL)
            delete file1;

        if(file2 != NULL)
            delete file2;

        file1 = NULL;
        file2 = NULL;

        return file_pointers;

    }

    //Parse file content
    file1->parse();
    file2->parse();

    file_pointers.push_back(file1);
    file_pointers.push_back(file2);

    return file_pointers;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void printPointValues

(

points *

pPoints

)

Write point values to standard out.

Definition at line 259 of file main.cpp.

References points::get_num_points(), points::get_point(), pnt::size, and pnt::vals.

                                      {

    //Print points
    for(int point_num = 0; point_num < pPoints->get_num_points(); point_num++){
        pnt curr = pPoints->get_point(point_num);

        for(int dimension = 0; dimension < curr.size; dimension++)
            cout << curr.vals[dimension] << ", ";
        cout << endl;
    }

    cout << endl << "Printed " << pPoints->get_num_points() << " points" << endl;
}

Here is the call graph for this function:

void printUsage

(

string

rErrorMessage

)

Print the rc_cmp usage message.

Parameters

rErrorString

The error message

Definition at line 326 of file main.cpp.

References gUsage(), print_metrics(), and Fileselect::printFiletypes().

                                     {

    if(rErrorMessage.compare("") != 0)
        cout << "\n" << rErrorMessage << endl;
    cout << "\n" << gUsage << endl;

    cout << "Comparisons are made at the point locations of input file 1" << endl << endl;

    cout << "Valid Filetypes:\n";
    Fileselect::printFiletypes();
    cout << endl;

    cout << "Metrics:\n";
    print_metrics();
    cout << endl;

}

Here is the call graph for this function:

std::vector<points*> retrievePoints	(	datafile *	pFile,
		vector< int >	rPartitions
	)

Construct a vector of pointers to datasets (point objects) from a list of partition numbers.

Definition at line 240 of file main.cpp.

References datafile::get_num_partitions(), and datafile::get_points().

Referenced by compareFiles().

                                                                           {
    
    std::vector<points*> datasets;
    
    for(int partition_record = 0; partition_record < rPartitions.size(); partition_record++){
        if(rPartitions[partition_record] < pFile->get_num_partitions()){
            datasets.push_back( pFile->get_points(rPartitions[partition_record]) );
        }
        else
            datasets.push_back( NULL );
    }

    return datasets;
}

Here is the call graph for this function:

Here is the caller graph for this function: