ASMmxBase.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ASM_MX_BASE_H
#define _ASM_MX_BASE_H
 
#include "MatVec.h"
#include <memory>
 
namespace Go {
  class SplineSurface;
  class SplineVolume;
}
 
namespace SplineUtils {
  enum class AdjustOp;
}
 
 
class ASMmxBase
{
protected:
  explicit ASMmxBase(const std::vector<unsigned char>& n_f) : nfx(n_f) {}
 
  void initMx(const std::vector<int>& MLGN, const int* sysMadof);
 
  void extractNodeVecMx(const RealArray& glbVec, RealArray& nodVec,
                        int basis) const;
 
  void injectNodeVecMx(RealArray& glbVec, const RealArray& nodVec,
                       int basis) const;
 
  bool getSolutionMx(Matrix& sField, const Vector& locSol,
                     const std::vector<int>& nodes) const;
 
public:
  enum MixedType {
    NONE = 0,                  
    FULL_CONT_RAISE_BASIS1,    
    REDUCED_CONT_RAISE_BASIS1, 
    FULL_CONT_RAISE_BASIS2,    
    REDUCED_CONT_RAISE_BASIS2, 
    DIV_COMPATIBLE,            
    SUBGRID,                   
  };
 
  static MixedType Type; 
  static char  itgBasis; 
  bool piola = false;    
  static bool includeExtra; 
 
protected:
  typedef std::vector<std::shared_ptr<Go::SplineSurface>> SurfaceVec; 
  typedef std::vector<std::shared_ptr<Go::SplineVolume>> VolumeVec; 
 
  static SurfaceVec establishBases(Go::SplineSurface* surf, MixedType type);
 
  static VolumeVec establishBases(Go::SplineVolume* svol, MixedType type);
 
  static Go::SplineSurface* adjustBasis(const Go::SplineSurface& surf,
                                        const std::array<SplineUtils::AdjustOp,2>& ops);
  static Go::SplineVolume* adjustBasis(const Go::SplineVolume& svol,
                                       const std::array<SplineUtils::AdjustOp,3>& ops);
 
private:
  std::vector<int> MADOF; 
 
protected:
  std::vector<size_t> elem_size;
  std::vector<size_t> nb;
 
  std::vector<unsigned char> nfx; 
};
 
#endif