ossimLagrangeInterpolator.cpp File Reference

#include <iomanip>
#include <string>
#include "ossim/base/ossimLagrangeInterpolator.h"
#include "ossim/base/ossimTrace.h"

Go to the source code of this file.

Functions
ostream &	operator<< (ostream &stream, const ossimLagrangeInterpolator &interpolator)
istream &	operator>> (istream &stream, ossimLagrangeInterpolator &interpolator)

Function Documentation

operator<<()

ostream& operator<<	(	ostream &	stream,
		const ossimLagrangeInterpolator &	interpolator
	)

Definition at line 124 of file ossimLagrangeInterpolator.cpp.

{
   stream << "\n" << HEADER_LABEL
          << "  " << interpolator.theTeeArray.size()
          << "  " << interpolator.theNumElements << endl;

   for (ossim_uint32 i=0; i<interpolator.theTeeArray.size(); ++i)
   {
      stream << setprecision(20);
      stream << interpolator.theTeeArray[i]   << " "
             << interpolator.theNormalizer[i] << endl
             << "( "<<interpolator.theDataArray[i](0)
             << interpolator.theDataArray[i](1)
             << interpolator.theDataArray[2](1) << " )" << endl;
   }
   
   return stream;
}

operator>>()

istream& operator>>	(	istream &	stream,
		ossimLagrangeInterpolator &	interpolator
	)

Definition at line 148 of file ossimLagrangeInterpolator.cpp.

References ossimLagrangeInterpolator::theDataArray, ossimLagrangeInterpolator::theNormalizer, ossimLagrangeInterpolator::theNumElements, and ossimLagrangeInterpolator::theTeeArray.

{
   // Start with ensuring we're at the right spot.  Let it go a max of ten
   // reads.  If we go beyond that there's a problem...
   string s;
   int count = 0;
   stream >> s;
   while (s != "LAGRANGE_INTERP") 
   {
      stream >> s;
      ++count;
      if (count == 10) break; // Avoid infinite loop...
   }

   if (count == 10)
   {
      cerr << "ossimLagrangeInterpolator operator>> ERROR:"
           << "\nStream seems to not contain LAGRANGE_INTERP header!"
           << "\nReturning..." << endl;
      return stream;
   }

   ossim_uint32 numPoints;
   stream >> numPoints >> interpolator.theNumElements;

   interpolator.theTeeArray.clear();
   interpolator.theDataArray.clear();
   interpolator.theNormalizer.clear();

   double tee, normalizer;
   NEWMAT::ColumnVector v(3);
   for (ossim_uint32 i=0; i<numPoints; i++)
   {
      stream >> tee >> v(0) >> v(1) >> v(2) >> normalizer;
      if (!stream.fail())
      {
         interpolator.theTeeArray.push_back(tee);
         interpolator.theDataArray.push_back(v);
         interpolator.theNormalizer.push_back(normalizer);
      }
   }

   return stream;
}