Statistics.h File Reference

#include "StatisticsVector.h"

Go to the source code of this file.

Functions
void	Statistics (int argc, char *argv[])
	This is the function that the user should call for statistics (CG). More...

Function Documentation

Statistics()

void Statistics	(	int	argc,
		char *	argv[]
	)

This is the function that the user should call for statistics (CG).

This function offers an interface to the templated class StatisticsVector. It picks the correct StatType template from the user-specified flag. The other flags are processed by StatisticsVector.

{
    if (argc > 1 && strcmp(argv[1], "-help"))
        logger(INFO, "\nGet statistics for %", argv[1]);
    
    /* Check for the '-stattype' flag */
    StatType T = O; //default value
    for (unsigned int i = 2; i < argc; i++)
    {
        if (!strcmp(argv[i], "-stattype") || !strcmp(argv[i], "-statType"))
        {
            if (!strcmp(argv[i + 1], "XYZ"))
                T = XYZ;
            else if (!strcmp(argv[i + 1], "RAZ"))
                T = RAZ;
            else if (!strcmp(argv[i + 1], "RA"))
                T = RA;
            else if (!strcmp(argv[i + 1], "RZ"))
                T = RZ;
            else if (!strcmp(argv[i + 1], "AZ"))
                T = AZ;
            else if (!strcmp(argv[i + 1], "R"))
                T = R;
            else if (!strcmp(argv[i + 1], "A"))
                T = A;
            else if (!strcmp(argv[i + 1], "XY"))
                T = XY;
            else if (!strcmp(argv[i + 1], "XZ"))
                T = XZ;
            else if (!strcmp(argv[i + 1], "YZ"))
                T = YZ;
            else if (!strcmp(argv[i + 1], "X"))
                T = X;
            else if (!strcmp(argv[i + 1], "Y"))
                T = Y;
            else if (!strcmp(argv[i + 1], "Z"))
                T = Z;
            else if (!strcmp(argv[i + 1], "O"))
                T = O;
            else
            {
                logger(ERROR, "stattype unknown");
            }
        }
    }
    if (T == XY)
    { // averaging in z-direction
        logger(INFO, "averaging in z-direction");
        StatisticsVector<XY> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == XZ)
    { // averaging in y-direction
        logger(INFO, "averaging in y-direction");
        StatisticsVector<XZ> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == YZ)
    { // averaging in x-direction
        logger(INFO, "averaging in x-direction");
        StatisticsVector<YZ> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == X)
    { // averaging in yz-direction
        logger(INFO, "averaging in yz-direction");
        StatisticsVector<X> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == Y)
    { // averaging in yz-direction
        logger(INFO, "averaging in xz-direction");
        StatisticsVector<Y> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == Z)
    { // averaging in yz-direction
        logger(INFO, "averaging in xy-direction");
        StatisticsVector<Z> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == O)
    { // averaging in all directions
        logger(INFO, "averaging in xyz-direction");
        StatisticsVector<O> stats(argc, argv);
        stats.setDoPeriodicWalls(false);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == RAZ)
    { //no averaging
        logger(INFO, "cylindrical, no averaging");
        StatisticsVector<RAZ> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == RA)
    { //no averaging
        logger(INFO, "cylindrical, Z averaging");
        StatisticsVector<RA> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == RZ)
    { //no averaging
        logger(INFO, "cylindrical, A averaging");
        StatisticsVector<RZ> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == AZ)
    { //no averaging
        logger(INFO, "cylindrical, R averaging");
        StatisticsVector<AZ> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == A)
    { //no averaging
        logger(INFO, "cylindrical, RZ averaging");
        StatisticsVector<A> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == R)
    { //no averaging
        logger(INFO, "cylindrical, AZ averaging");
        StatisticsVector<R> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    else if (T == XYZ)
    { //no averaging
        logger(INFO, "no spatial averaging");
        StatisticsVector<XYZ> stats(argc, argv);
        stats.statistics_from_fstat_and_data();
    }
    
}

References A, AZ, ERROR, constants::i, INFO, logger, O, R, RA, RAZ, RZ, StatisticsVector< T >::setDoPeriodicWalls(), StatisticsVector< T >::statistics_from_fstat_and_data(), X, XY, XYZ, XZ, Y, YZ, and Z.

Referenced by main().