GridWithCartesianVertices.hpp

#ifndef __FIBER_GRID_TYPES_GRIDWITHCARTESIANVERTICES_HPP
#define __FIBER_GRID_TYPES_GRIDWITHCARTESIANVERTICES_HPP
 
#include "TypedGrid.hpp"
 
 
namespace Fiber
{
 
struct  gridtypes_API   GridWithCartesianVertices : virtual public TypedGrid
{
        RefPtr<Skeleton>
                Vertices;
 
        RefPtr<Representation>
                CartesianVertices; 
 
        RefPtr<Field>
                CartesianPositions; 
 
        GridWithCartesianVertices(const RefPtr<Grid>&theGrid = nullptr)
        {
                assign( theGrid );
        }
 
        bool    assign(const RefPtr<Grid>&theGrid)
        {
                Vertices           = nullptr; 
                CartesianVertices  = nullptr; 
                CartesianPositions = nullptr;
 
                if (!TypedGrid::assign(theGrid) ) return false;
 
                Vertices = theGrid->findVertices();
                if (!Vertices)
                        return false; 
 
                CartesianVertices = theGrid->getCartesianRepresentation(*Vertices ); 
                if (!CartesianVertices)
                        return false;
 
                CartesianPositions = theGrid->CartesianPositions(); 
                if (!CartesianPositions)
                        return false; 
 
                return true;
        } 
 
        ~GridWithCartesianVertices();
 
 
        bool    isValid() const
        {
                return CartesianPositions.valid();
        }
 
static  GridWithCartesianVertices
        create( Grid&myGrid, int Dims, const RefPtr<Field>&Positions)
        {
        Representation&myCartesianRepresentation = myGrid.makeCartesianRepresentation(Dims);
                myCartesianRepresentation[ FIBER_POSITIONS ] = Positions; 
 
                return GridWithCartesianVertices( myGrid.self() );
        }
};
 
 
} // namespace Fiber
 
#endif // __FIBER_GRID_TYPES_GRIDWITHCARTESIANVERTICES_HPP