hack_blitztotable.h

Go to the documentation of this file.

// include guard
#ifndef TF_HACK_BLITZTOTABLE_H_VERSION

#define TF_HACK_BLITZTOTABLE_H_VERSION \
  "TF_HACK_BLITZTOTABLE_H   V1.0   "

#include<fstream>
#include<cmath>
#include<blitz/array.h>

namespace tfxx {
namespace hacks {

template<class T>
void write_blitztable(const blitz::Array<T, 1>& array,
                      const char* filename,
                      const double& dx=1.)
{
  std::ofstream os(filename);
  for (int i=array.lbound(blitz::firstDim);
       i<=array.ubound(blitz::secondDim); i++)
  {
    os << dx*(i-1) << " " 
      << array(i) << " "
      << std::abs(array(i)) << std::endl;
  }
}

} // namespace hacks
} // namespace tfxx

#endif // TF_HACK_BLITZTOTABLE_H_VERSION (includeguard)

/* ----- END OF hack_blitztotable.h ----- */