Line data    Source code

       1             : #include <synthesis/MeasurementComponents/GlobalFTMachineCallbacks.h>
       2             : #include <synthesis/TransformMachines/IlluminationConvFunc.h>
       3             : 
       4             : using namespace casacore;
       5             : using namespace casa;
       6             : 
       7             : extern "C" 
       8             : {
       9             :   //
      10             :   // The Gridding Convolution Function (GCF) used by the underlying
      11             :   // gridder written in FORTRAN.
      12             :   //
      13             :   // The arguments must all be pointers and the value of the GCF at
      14             :   // the given (u,v) point is returned in the weight variable.  Making
      15             :   // this a function which returns a complex value (namely the weight)
      16             :   // has problems when called in FORTRAN - I (SB) don't understand
      17             :   // why.
      18             :   //
      19             :   //---------------------------------------------------------------
      20             :   //
      21             :   IlluminationConvFunc casa::gwEij;
      22           0 :   void gcppeij(Double *griduvw, Double *area,
      23             :                Double *raoff1, Double *decoff1,
      24             :                Double *raoff2, Double *decoff2, 
      25             :                Int *doGrad,
      26             :                Complex *weight,
      27             :                Complex *dweight1,
      28             :                Complex *dweight2,
      29             :                Double *currentCFPA)
      30             :   {
      31           0 :     Complex w,d1,d2;
      32           0 :     gwEij.getValue(griduvw, raoff1, raoff2, decoff1, decoff2,
      33             :                    area,doGrad,w,d1,d2,*currentCFPA);
      34           0 :     *weight   = w;
      35           0 :     *dweight1 = d1;
      36           0 :     *dweight2 = d2;
      37           0 :   }
      38             : };