---

ExactRaytracer.hpp

/*
 *  $Id: ExactRaytracer_8hpp-source.html,v 1.1.1.1 2002/06/12 22:58:44 cgagne Exp $
 *
 *  LOSR -- Library for Optical Systems Raytracing
 *  Copyright (C) 2001 by Christian Gagné and Julie Beaulieu
 *
 *  This library is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU Lesser General Public
 *  License as published by the Free Software Foundation; either
 *  version 2.1 of the License, or (at your option) any later version.
 *
 *  This library is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *  Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public
 *  License along with this library; if not, write to the Free Software
 *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

#ifndef LOSR_ExactRaytracer_hpp
#define LOSR_ExactRaytracer_hpp

#include <vector>
#include <iostream>

#include "losr/Surface.hpp"
#include "losr/OpticalSystem.hpp"

namespace LOSR {

struct QU {
  QU();
  QU(const LOSR::Surface& inSurface);

  QU& operator=(const LOSR::Surface& inSurface);

  long double mT;      
  long double mN;      
  long double mC;      
  long double mR;      
  long double mY;      
  long double mU;      
  long double mQ;      
  long double mSinU;   
  long double mCosU;   
  long double mSinI;   
  long double mCosI;   
  long double mSinI_U; 
  long double mCosI_U; 
  long double mSinIp;  
  long double mCosIp;  
  long double mSinUp;  
  long double mCosUp;  
  long double mQp;     
};


class ExactRaytracer : public std::vector<QU> {

public: 
  ExactRaytracer(unsigned int inSize=0);
  ExactRaytracer(const LOSR::OpticalSystem& inOpticalSystem);
  ~ExactRaytracer() { }   

  ExactRaytracer& operator=(const LOSR::OpticalSystem& inOpticalSystem);

  void initialize(const LOSR::OpticalSystem& inOpticalSystem);
  void read(std::istream& ioIs=std::cin);
  bool trace(ExactRaytracer::iterator inStartingPoint);
  bool trace(ExactRaytracer::iterator inStartingPoint,long double inY,long double inU);
  bool traceBackward(ExactRaytracer::iterator inActualEntry,ExactRaytracer::iterator inStopEntry);
  bool traceForward(ExactRaytracer::iterator inActualEntry,ExactRaytracer::iterator inStopEntry);
  void write(std::ostream& ioOs=std::cout) const;

};

}


std::istream& operator>>(std::istream& ioIs, LOSR::QU& outQU);
std::ostream& operator<<(std::ostream& ioOs, const LOSR::QU& inQU);

inline std::istream& operator>>(
  std::istream& ioIs,                     
  LOSR::ExactRaytracer& outExactRaytracer 
)
{
  outExactRaytracer.read(ioIs);
  return ioIs;
}


inline std::ostream& operator<<(
  std::ostream& ioOs,                     
  const LOSR::ExactRaytracer& inExactRaytracer  
)
{
  inExactRaytracer.write(ioOs);
  return ioOs;
}


#endif // LOSR_ExactRaytracer_hpp

---