AssembleFEBlock_decl.hpp

#ifndef AssembleFEBlock_DECL_hpp
#define AssembleFEBlock_DECL_hpp
 
 
#include "feddlib/core/FEDDCore.hpp"
#include "feddlib/core/LinearAlgebra/Matrix.hpp"
#include "feddlib/core/FE/Helper.hpp"
#include "feddlib/core/AceFemAssembly/AssembleFE.hpp"
 
 
namespace FEDD {
    
    template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
    class AssembleFEBlock: public AssembleFE<SC,LO,GO,NO>  {
    public:
 
 
        typedef SmallMatrix<SC> SmallMatrix_Type;
        typedef Teuchos::RCP<SmallMatrix_Type> SmallMatrixPtr_Type;
 
        // Teuchos:: Array anstatt vec_dbl_Type
 
        typedef AssembleFEBlock<SC,LO,GO,NO> AssembleFEBlock_Type;

        virtual void assembleJacobian() = 0;

        virtual void assembleRHS() = 0;
 
        virtual void assembleJacobianBlock(LO i) =0;
        //virtual void assembleMass(MatrixPtr_Type &A) =0;
 
    protected:

        AssembleFEBlock(int flag,
                   vec2D_dbl_Type nodesRefConfig,
                   ParameterListPtr_Type parameters,
                tuple_disk_vec_ptr_Type tuple);
 
        //void readTuple(); /// TODO: To have tuple information in basis class as well?
        //tuple_disk_vec_ptr_Type getTuple();  
    private:
 
        void assembleMonolithicSystem(SmallMatrixPtr_Type elementMatrix);
        
        int dimSystem_;
 
        std::string FEType_ ; // FEType of Disk
 
        int dofsSolid_ ; // Degrees of freedom per node
        int dofsChem_;
        int numNodesSolid_ ; // Number of nodes of element
        int numNodesChem_ ; // Number of nodes of element
 
 
        int dofsElement_; // "Dimension of return matrix"
        
 
        friend class AssembleFEFactory<SC,LO,GO,NO>;
    };
}
#endif