SLIP  1.4
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
regularvector3dfield3d_op.cpp
#include <iostream>
#include <vector>
#include <numeric>
#include <algorithm>
#include <slip/arithmetic_op.hpp>
#include <slip/histo.hpp>
#include <slip/Volume.hpp>
//#include <slip/border_treatment.hpp>
#include <slip/RegularVector3dField3d.hpp>
#include <slip/dynamic.hpp>
#include <slip/statistics.hpp>
#include <cstdlib>
#include <ctime>
static void usage(const char * );
int main(int argc, char **argv)
{
//-----------------------------------------------------------------
//get the command line argmuments
//-----------------------------------------------------------------
std::vector<std::string> args;
for(int i = 0; i < argc; i++)
{
args.push_back(argv[i]);
}
if((args.size() == 1) || (args[1] == "-h") || (args[1] == "--help"))
{
usage(args[0].c_str());
exit(0);
}
//order of finite difference derivative schemes
const std::size_t order = static_cast<std::size_t>(atoi(argv[4]));
//------------------------------------------------------------------
if(args.size() == 13)
{
//-----------------------------------------------------------------
//read the vector field
//-----------------------------------------------------------------
slip::RegularVector3dField3d<double> VF;
VF.read_tecplot(args[2]);
//store the number of rows and columns of the vector field
const std::size_t slices = VF.slices();
const std::size_t rows = VF.rows();
const std::size_t cols = VF.cols();
//-----------------------------------------------------------------
//computes the (Euclidean) norm map of the vector field
//and get some statistics from it
//-----------------------------------------------------------------
slip::Volume<double> Norm(slices,rows,cols);
VF.norm(Norm);
double sum_n = Norm.sum();
double energy = sum_n * sum_n;
double mean_n = slip::mean<double>(Norm.begin(),Norm.end());
double std_n = slip::std_dev<double>(Norm.begin(),Norm.end(),mean_n);
double min_n = *std::min_element(Norm.begin(),Norm.end());
double median_n = slip::median_from_data_n(Norm.begin(),Norm.end(),Norm.size());
double max_n = *std::max_element(Norm.begin(),Norm.end());
double rms_n = slip::rms<double>(Norm.begin(),Norm.end());
//-----------------------------------------------------------------
//computes a real histogram of norm
//-----------------------------------------------------------------
//create the steps array [min_n, min_n + step_n,...,max_n]
slip::Array<double> histo_step_n(9);
double step_n = (max_n-min_n)/double(histo_step_n.size()-1);
slip::iota(histo_step_n.begin(),histo_step_n.end(),min_n,step_n);
//computes the histogram
slip::Array<int> histo_n(histo_step_n.size()-1);
slip::real_histogram(Norm.begin(),Norm.end(),
histo_step_n.begin(),
histo_step_n.end(),histo_n.begin());
std::cout<<"Vector field energy = "<<energy<<std::endl;
std::cout<<"mean_n = "<<mean_n<<std::endl;
std::cout<<"std dev n = "<<std_n<<std::endl;
std::cout<<"[min_n,median_n,max_n] = "<<"["<<min_n<<","<<median_n<<","<<max_n<<"]"<<std::endl;
std::cout<<"rms_n = "<<rms_n<<std::endl;
std::cout<<"histo_n ="<<std::endl;
std::cout<<histo_n<<std::endl;
//-----------------------------------------------------------------
//computes some statistical informations on the first component u
//of the vector field
//-----------------------------------------------------------------
double mean_u = slip::mean<double>(VF.begin(0),VF.end(0));
double std_u = slip::std_dev<double>(VF.begin(0),VF.end(0),mean_u);
double min_u = *std::min_element(VF.begin(0),VF.end(0));
double median_u = slip::median_from_data_n(VF.begin(0),VF.end(0),VF.size());
double max_u = *std::max_element(VF.begin(0),VF.end(0));
double rms_u = slip::rms<double>(VF.begin(0),VF.end(0));
//create the steps array [min_u, min_u + step_u,...,max_u]
slip::Array<double> histo_step_u(9);
double step_u = (max_u-min_u)/double(histo_step_u.size()-1);
slip::iota(histo_step_u.begin(),histo_step_u.end(),min_u,step_u);
//computes the histogram
slip::Array<int> histo_u(histo_step_u.size()-1);
slip::real_histogram(VF.begin(0),VF.end(0),
histo_step_u.begin(),
histo_step_u.end(),histo_u.begin());
std::cout<<"mean_u = "<<mean_u<<std::endl;
//std::cout<<"rms_u = "<<rms_u<<std::endl;
std::cout<<"std dev u = "<<std_u<<std::endl;
std::cout<<"[min_u,median_u,max_u] = "<<"["<<min_u<<","<<median_u<<","<<max_u<<"]"<<std::endl;
std::cout<<"rms_u = "<<rms_u<<std::endl;
std::cout<<"histo_u ="<<std::endl;
std::cout<<histo_u<<std::endl;
//-----------------------------------------------------------------
//computes some statistical informations on the second component v
//of the vector field
//-----------------------------------------------------------------
double mean_v = slip::mean<double>(VF.begin(1),VF.end(1));
double std_v = slip::std_dev<double>(VF.begin(1),VF.end(1),mean_v);
double min_v = *std::min_element(VF.begin(1),VF.end(1));
double median_v = slip::median_from_data_n(VF.begin(1),VF.end(1),VF.size());
double max_v = *std::max_element(VF.begin(1),VF.end(1));
double rms_v = slip::rms<double>(VF.begin(1),VF.end(1));
//create the steps array [min_v, min_v + step_v,...,max_v]
slip::Array<double> histo_step_v(9);
double step_v = (max_v-min_v)/double(histo_step_v.size()-1);
slip::iota(histo_step_v.begin(),histo_step_v.end(),min_v,step_v);
//computes the histogram
slip::Array<int> histo_v(histo_step_v.size()-1);
slip::real_histogram(VF.begin(1),VF.end(1),
histo_step_v.begin(),
histo_step_v.end(),histo_v.begin());
std::cout<<"mean_v = "<<mean_v<<std::endl;
std::cout<<"std dev u = "<<std_v<<std::endl;
std::cout<<"[min_v,median_v,max_v] = "<<"["<<min_v<<","<<median_v<<","<<max_v<<"]"<<std::endl;
std::cout<<"rms_v = "<<rms_v<<std::endl;
std::cout<<"histo_v ="<<std::endl;
std::cout<<histo_v<<std::endl;
//-----------------------------------------------------------------
//computes some statistical informations on the third component w
//of the vector field
//-----------------------------------------------------------------
double mean_w = slip::mean<double>(VF.begin(2),VF.end(2));
double std_w = slip::std_dev<double>(VF.begin(2),VF.end(2),mean_w);
double min_w = *std::min_element(VF.begin(2),VF.end(2));
double median_w = slip::median_from_data_n(VF.begin(2),VF.end(2),VF.size());
double max_w = *std::max_element(VF.begin(2),VF.end(2));
double rms_w = slip::rms<double>(VF.begin(2),VF.end(2));
//create the steps array [min_w, min_w + step_w,...,max_w]
slip::Array<double> histo_step_w(9);
double step_w = (max_w-min_w)/double(histo_step_w.size()-1);
slip::iota(histo_step_w.begin(),histo_step_w.end(),min_w,step_w);
//computes the histogram
slip::Array<int> histo_w(histo_step_w.size()-1);
slip::real_histogram(Norm.begin(),Norm.end(),
histo_step_w.begin(),
histo_step_w.end(),histo_w.begin());
std::cout<<"mean_w = "<<mean_w<<std::endl;
//std::cout<<"rms_w = "<<rms_w<<std::endl;
std::cout<<"std dev u = "<<std_w<<std::endl;
std::cout<<"[min_w,median_w,max_w] = "<<"["<<min_w<<","<<median_w<<","<<max_w<<"]"<<std::endl;
std::cout<<"rms_w = "<<rms_w<<std::endl;
std::cout<<"histo_w ="<<std::endl;
std::cout<<histo_w<<std::endl;
//-----------------------------------------------------------------
//compute the divergence with a scheme of order order
//-----------------------------------------------------------------
slip::Volume<double> Div(slices,rows,cols);
VF.divergence(order,Div);
//-----------------------------------------------------------------
//compute the vorticity with a scheme of order order
//-----------------------------------------------------------------
slip::RegularVector3dField3d<double> Vort(slices,rows,cols);
VF.vorticity(order,Vort);
//-----------------------------------------------------------------
//compute the vorticity with a scheme of order order
//-----------------------------------------------------------------
slip::RegularVector3dField3d<double> Lambda(slices,rows,cols);
VF.lambda(order,Lambda);
//-----------------------------------------------------------------
//Put the images between [0,1]
//-----------------------------------------------------------------
slip::change_dynamic_01(Norm.begin(),Norm.end(),Norm.begin(),slip::AFFINE_FUNCTION);
slip::change_dynamic_01(Div.begin(),Div.end(),Div.begin(),slip::SIGMOID_FUNCTION);
//-----------------------------------------------------------------
//Write the images
//-----------------------------------------------------------------
Norm.write_to_images(args[6],0,(slices-1));
Div.write_to_images(args[8],0,(slices-1));
Vort.write_tecplot(args[10],"title","zone");
Lambda.write_tecplot(args[12],"title","zone");
}
else
{
usage(args[0].c_str());
exit(1);
}
return 0;
}
static void usage(const char* cmd)
{
std::cout<<std::endl;
std::cout<<"NAME:"<<std::endl;
std::cout<<"\t"<<cmd<<" - Computes some statistics and differential operators on a regular vector field"<<std::endl;
std::cout<<std::endl;
std::cout<<"SYNOPSIS:"<<std::endl;
std::cout<<"\t"<<cmd<<" -i input_tecplot_file -n derivative_scheme_order -o norm_image_file -o divergence_image_file -o vorticity_tecplot_file -o lambda2_tecplot_file"<<std::endl;
std::cout<<std::endl;
std::cout<<"DESCRIPTION: Computes some statistics and differential operators on a regular vector field"<<std::endl;
std::cout<<"\tThe finite difference derivatives scheme order is specify by the -n order option"<<std::endl;
std::cout<<"\t "<<std::endl;
std::cout<<std::endl;
std::cout<<"\t -i input gnuplot_file"<<std::endl;
std::cout<<"\t -n finite difference derivative scheme order of the divergence and the vorticity"<<std::endl;
std::cout<<"\t -o norm output image file"<<std::endl;
std::cout<<"\t -o divergence output image file"<<std::endl;
std::cout<<"\t -o vorticity output image file"<<std::endl;
std::cout<<"\t -o lambda2 output image file"<<std::endl;
std::cout<<"\t -h, --help\t display this help and exit"<<std::endl;
std::cout<<std::endl;
std::cout<<"AUTHOR:"<<std::endl;
std::cout<<"\t Benoit Tremblais <tremblais@sic.univ-poitiers.fr>"<<std::endl;
std::cout<<std::endl;
std::cout<<"REPORTING BUGS:"<<std::endl;
std::cout<<"\t Report bugs to slip trac serveur"<<std::endl;
std::cout<<std::endl;
std::cout<<"COPYRIGHT:"<<std::endl;
std::cout<<"\t This is free software; see the source for copying conditions."<<std::endl;
std::cout<<std::endl;
std::cout<<"EXAMPLE: "<<std::endl;
std::cout<<"./regularvector3dfield3d_op -i 3D3C.dat -n 2 -o norm.png -o div.png -o vort.dat -o lambda2.dat"<<std::endl;
std::cout<<"\t "<<std::endl;
std::cout<<std::endl;
std::cout<<"SEE ALSO: "<<std::endl;
std::cout<<std::endl;
}
Generated on Tue Sep 30 2014 10:05:55 for SLIP by   doxygen 1.8.6