write_grepg_c()

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

create grepg file from complex array

Parameters

array	2D complex float blitz array
filename	name of file to write to
dx	x-sampling to simulate
dy	y-sampling to simulate

Definition at line 48 of file hack_blitztogrepg.cc.

References tfxx::fortranio::FortranBinOutput::end_block(), and tfxx::ioswap::file_magic_write().

Referenced by write_grepg(), and write_grepg_f().

{
  std::ofstream os(filename);
  tfxx::ioswap::file_magic_write(os, "1234", true);
  tfxx::fortranio::FortranBinOutput fos(os);
  int nom=array.extent(blitz::firstDim);
  int fom=array.lbound(blitz::firstDim);
  int nslo=array.extent(blitz::secondDim);
  int fslo=array.lbound(blitz::secondDim);
  fos << nom << nslo;
  fos.end_block();
  for (int i=0; i<nom; i++) 
  { fos << float(3.14159265358979311*2.*(i+fom-1)*dx); }
  for (int i=0; i<nslo; i++) { fos << float(1.e-3*(i+fslo-1)*dy); }
  fos.end_block();
  for (int i=0; i<nslo; i++)
  {
    for (int j=0; j<nom; j++)
    {
      fos << array(j+fom,i+fslo);
    }
  }
  fos.end_block();
}

Here is the call graph for this function:

Here is the caller graph for this function: