00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
00016 
00017 
00018 
00019 
00020 
00021 
00022 
00023 
00024 
00025 
00026 
00027 
00028 
00029 
00030 
00031 #ifndef SUNDANCE_FIELDBASE_H
00032 #define SUNDANCE_FIELDBASE_H
00033 
00034 #include "SundanceDefs.hpp"
00035 #include "PlayaHandleable.hpp"
00036 #include "Teuchos_Array.hpp"
00037 #include "Teuchos_RefCountPtr.hpp"
00038 #include "SundanceCellSet.hpp"
00039 
00040 
00041 
00042 namespace Sundance
00043 {
00044 class CellFilter;
00045 }
00046 
00047 
00048 namespace Sundance
00049 {
00050 class Mesh;
00051 using Sundance::CellFilter;
00052 
00053 using Sundance::Map;
00054 using Teuchos::Array;
00055 using Teuchos::RefCountPtr;
00056 
00057 
00058 
00059 class FieldBase : public Playa::Handleable<FieldBase>
00060 {
00061 public:
00062 
00063   FieldBase(){;}
00064 
00065 
00066   virtual ~FieldBase(){;}
00067 
00068 
00069   virtual int numElems() const {return 1;}
00070 
00071 
00072   virtual double getData(int cellDim, int cellID, int elem) const = 0 ;
00073 
00074 
00075   virtual bool isDefined(int cellDim, int cellID, int elem) const = 0 ;
00076 
00077 
00078   virtual bool isPointData() const = 0 ;
00079 
00080 
00081   virtual bool isCellData() const {return !isPointData();}
00082 
00083 
00084 
00085 
00086 
00087   virtual void getDataBatch(int cellDim, const Array<int>& cellID,
00088     const Array<int>& funcElem, Array<double>& batch) const ;
00089 
00090 
00091 
00092 
00093   virtual const CellFilter& domain() const ;
00094 
00095 };
00096 
00097 using Sundance::CellSet;
00098 using Sundance::CellFilter;
00099 
00100 CellSet connectedNodeSet(const CellFilter& f, const Mesh& mesh);
00101 RCP<Array<int> > cellSetToLIDArray(const CellSet& cs);
00102 }
00103 
00104 #endif