fei_FiniteElementData.hpp

#ifndef _fei_FiniteElementData_hpp_
#define _fei_FiniteElementData_hpp_

#include <fei_Lookup.hpp>

class FiniteElementData {
 public:

  virtual ~FiniteElementData() {};

  virtual int parameters(int numParams, char** params) = 0;


  virtual int setLookup(Lookup& lookup) = 0;


  virtual int describeStructure(int numElemBlocks,
                                const int* numElemsPerBlock,
                                const int* numNodesPerElem,
                                const int* elemMatrixSizePerBlock,
                                int totalNumNodes,
                                int numSharedNodes,
                                int numMultCRs) = 0;

   virtual int setConnectivity(int elemBlockID,
                               int elemID,
                               int numNodes,
                               const int* nodeNumbers,
                               const int* numDofPerNode,
                               const int* dof_ids) = 0;

   virtual int setElemMatrix(int elemBlockID,
                             int elemID,
                             int numNodes,
                             const int* nodeNumbers,
                             const int* numDofPerNode,
                             const int* dof_ids,
                             const double *const * coefs) = 0;

   virtual int setElemVector(int elemBlockID,
                             int elemID,
                             int numNodes,
                             const int* nodeNumbers,
                             const int* numDofPerNode,
                             const int* dof_ids,
                             const double* coefs) = 0;

   virtual int setDirichletBCs(int numBCs,
                               const int* nodeNumbers,
                               const int* dof_ids,
                               const double* values) = 0;

   virtual int sumIntoMatrix(int numRowNodes,
                             const int* rowNodeNumbers,
                             const int* row_dof_ids,
                             const int* numColNodesPerRow,
                             const int* colNodeNumbers,
                             const int* col_dof_ids,
                             const double* coefs) = 0;

   virtual int sumIntoRHSVector(int numNodes,
                             const int* nodeNumbers,
                             const int* dof_ids,
                             const double* coefs) = 0;

   virtual int putIntoRHSVector(int numNodes,
                             const int* nodeNumbers,
                             const int* dof_ids,
                             const double* coefs) = 0;

   virtual int loadComplete() = 0;


   virtual int launchSolver(int& solveStatus, int& iterations) = 0;

   virtual int reset() = 0;

   virtual int deleteConstraints() = 0;

   virtual int getSolnEntry(int nodeNumber,
                            int dof_id,
                            double& value) = 0;

   virtual int getMultiplierSoln(int CRID, double& lagrangeMultiplier) = 0;

   virtual int putNodalFieldData(int fieldID,
                                 int fieldSize,
                                 int numNodes,
                                 const int* nodeNumbers,
                                 const double* coefs) = 0;

   virtual int setMultiplierCR(int CRID,
                               int numNodes,
                               const int* nodeNumbers,
                               const int* dof_ids,
                               const double* coefWeights,
                               double rhsValue) = 0;

   virtual int setPenaltyCR(int CRID,
                            int numNodes,
                            const int* nodeNumbers,
                            const int* dof_ids,
                            const double* coefWeights,
                            double penaltyValue,
                            double rhsValue) = 0;
};

#endif