hack_blitztogrepg.h

Go to the documentation of this file.

// include guard
#ifndef TF_HACK_BLITZTOGREPG_H_VERSION

#define TF_HACK_BLITZTOGREPG_H_VERSION \
  "TF_HACK_BLITZTOGREPG_H   V1.1"

#include<complex>
#include<blitz/array.h>

namespace tfxx {
namespace hacks {

void write_grepg_c(const blitz::Array<std::complex<float>, 2>& array,
                   const char* filename,
                   const double& dx=1., const double& dy=1.);

void write_grepg_f(const blitz::Array<float, 2>& array,
                   const char* filename,
                   const double& dx=1., const double& dy=1.);

template<class T>
void write_grepg(const blitz::Array<T, 2>& array,
                 const char* filename,
                 const double& dx=1., const double& dy=1.)
{
  blitz::Array<std::complex<float>, 2> 
    carray(array.lbound(), array.shape(), blitz::fortranArray);
  carray=blitz::cast<std::complex<float> >(array);
  write_grepg_c(carray, filename, dx, dy);
}

} // namespace hacks
} // namespace tfxx

#endif // TF_HACK_BLITZTOGREPG_H_VERSION (includeguard)

/* ----- END OF hack_blitztogrepg.h ----- */