JEigenValues3D.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <vector>

#include "JGeometry3D/JVector3D.hh"
#include "JGeometry3D/JEigenValues3D.hh"

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"


/**
 * \file
 *
 * Example program to print eigen values in 3D.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;

  string  inputFile;
  int     debug;

  try { 

    JParser<> zap("Example program to print eigen values in 3D.");
    
    zap['f'] = make_field(inputFile);
    zap['d'] = make_field(debug)      = 1;
    
    zap(argc, argv);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }


  using namespace JPP;

  cout.tie(&cerr);


  vector<JVector3D> buffer;

  if (inputFile != "") {

    ifstream in(inputFile.c_str());

    for (double x, y, z; in >> x >> y >> z; ) {
      
      buffer.push_back(JVector3D(x,y,z));
    }
    
    in.close();
  }

  JEigenValues3D eigen(buffer.begin(), buffer.end());

  cout << "Number of eigen values "  << eigen.getN() << endl;
  cout << FIXED(7,3) << eigen.getX() << endl;
  cout << FIXED(7,3) << eigen.getY() << endl;
  cout << FIXED(7,3) << eigen.getZ() << endl;
}