IntxRllCssphere.hpp

Go to the documentation of this file.

/*
 * IntxRllCssphere.hpp
 *
 */
 
#ifndef INTXRLLCSSPHERE_HPP_
#define INTXRLLCSSPHERE_HPP_
 
#include "Intx2Mesh.hpp"
 
namespace moab
{
 
class IntxRllCssphere : public moab::Intx2Mesh
{
  public:
    IntxRllCssphere( Interface* mbimpl );
 
    virtual ~IntxRllCssphere();
 
    void set_radius( double radius )
    {
        R = radius;
    }
 
    double setup_tgt_cell( EntityHandle tgt, int& nsTgt );
 
    // src cell will be always lat lon cell, so it will be a rectangle in lat-lon coors
    // it will be used for "interior" determinations of other points
    // double setup_src_cell(EntityHandle src, int & nsSrc);
 
    // main method to intersect meshes on a sphere
 
    ErrorCode computeIntersectionBetweenTgtAndSrc( EntityHandle tgt,
                                                   EntityHandle src,
                                                   double* P,
                                                   int& nP,
                                                   double& area,
                                                   int markb[MAXEDGES],
                                                   int markr[MAXEDGES],
                                                   int& nsSrc,
                                                   int& nsTgt,
                                                   bool check_boxes_first = false );
 
    ErrorCode findNodes( EntityHandle tgt, int nsTgt, EntityHandle src, int nsSrc, double* iP, int nP );
 
  private:
    double R;            // radius of the sphere
    int plane;           // current gnomonic plane, will still be used for projection
    int srcEdgeType[4];  // at most 4
    // these could be from [-PI/2, +PI/2] and [0 to 2*PI]
};
 
} /* namespace moab */
#endif /* INTXRLLCSSPHERE_HPP_ */