/* Code for alternance positive algorithm */
////////////////////////////////////////////////////////////////////////////////
// export C interface
extern "C" 
void computeGold( float* reference, float* idata, 
                  const unsigned int size_x, const unsigned int size_y );

//////////////////////////////////////////////////////////////////////////////////
////! Compute reference data set
//////////////////////////////////////////////////////////////////////////////////
void
computeGold( float* reference, float* idata, 
            const unsigned int size_x, const unsigned int size_y ) 
{
  for( unsigned int y = 0; y < size_y; ++y) 
  {
    for( unsigned int x = 0; x < size_x; ++x) 
    {
      if( idata[(y * size_x) + x] < 0.00 )
      {
        reference[(y * size_x) + x] = 0.00;
      }
      else
      {
        reference[(y * size_x) + x] = idata[(y * size_x) + x];
      }
    }
  }
}