JGeometry.cc

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
#include <fstream>

#include "JAcoustics/JGeometry.hh"
#include "JAcoustics/JModel.hh"

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


/**
 * Example application for acoustics geometry.
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;

  JMODEL::JString parameters;
  JMechanics      mechanics;
  string          outputFile;
  int             debug;

  try {

    JParser<> zap("Example application for acoustics geometry.");

    zap['S'] = make_field(parameters)     = JPARSER::initialised();
    zap['M'] = make_field(mechanics)      = JPARSER::initialised();
    zap['o'] = make_field(outputFile)     = "";
    zap['d'] = make_field(debug)          = 3;

    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  // build string

  const double x0 =   0.0;
  const double y0 =   0.0;
  const double z0 =   0.0;
  const double z1 =  20.0;
  const double dz =  10.0;

  JGEOMETRY::JString string(JPosition3D(x0, y0, z0));

  string[0] = JGEOMETRY::JFloor(z0);

  for (int i = 1; i <= 18; ++i) {
    string[i] = JGEOMETRY::JFloor(z1 + (i - 1)*dz);   
  }

  string.mechanics = mechanics;

  for (int i = 1; i <= 18; ++i) {

    const JPosition3D pos[] = {
      string.getPosition(parameters, i - 1),
      string.getPosition(parameters, i)
    };

    cout << setw(2)    << i             << ' '
         << FIXED(7,3) << pos[1].getX() << ' ' 
         << FIXED(7,3) << pos[1].getY() << ' ' 
         << FIXED(7,3) << pos[1].getZ() << ' '
         << FIXED(7,3) << pos[1].getX() - pos[0].getX() << ' ' 
         << FIXED(7,3) << pos[1].getY() - pos[0].getY() << endl;
  }

  if (outputFile != "") {

    ofstream out(outputFile.c_str());

    for (int i = 1; i <= 18; ++i) {
      
      const JPosition3D pos = string.getPosition(parameters, i);

      out << FIXED(7,3) << pos.getX() << ' ' << FIXED(7,3) << pos.getZ() << endl;
    }

    out.close();
  }
}