ngsolve 6.2.2501.post16.dev1__cp313-cp313-macosx_10_15_universal2.whl → 6.2.2501.post34.dev1__cp313-cp313-macosx_10_15_universal2.whl

netgen/include/hdivdivsurfacespace.hpp

@@ -1,76 +0,0 @@
-#ifndef FILE_HDIVSYMSURFACESPACE
-#define FILE_HDIVSYMSURFACESPACE
-
-
-/*********************************************************************/
-/* File:   hdivdivsurfacespace.hpp                                   */
-/* Author: Michael Neunteufel, Joachim Schöberl, Astrid Pechstein    */
-/* Date:   April 2016                                                */
-/*********************************************************************/
-
-
-#include "fespace.hpp"
-
-namespace ngcomp
-{
-  
-  class HDivDivSurfaceSpace : public FESpace
-  {
-  private:
-    int nfa;
-    int nel;
-
-    int order;
-    Array<int> first_face_dof;
-    Array<int> first_element_dof;
-    int ndof;
-
-    Array<bool> fine_face; 
-
-    int noncontinuous;
-  public:
-
-    HDivDivSurfaceSpace(shared_ptr<MeshAccess> ama, const Flags & aflags, bool parseflags = false);
-
-    virtual ~HDivDivSurfaceSpace ();
-    static DocInfo GetDocu ();
-    
-
-    string GetClassName () const override
-    {
-      return "HDivSymSurfaceSpace";
-    }
-
-    void Update() override;
-
-    void UpdateCouplingDofArray() override;
-  
-    size_t GetNDof () const override
-    {
-      return ndof;
-    }
-
-    int GetFirstFaceDof (const int f) const { return first_face_dof[f]; }
-    int GetFirstCellDof (const int c) const { return first_element_dof[c]; }
-
-    FiniteElement & GetFE (ElementId ei, Allocator & lh) const override;
-    void GetDofNrs(ElementId id, Array<int> & dnums) const override;
-
-    virtual void GetWireBasketDofNrs (int elnr, Array<int> & dnums) const
-    { dnums.SetSize(0); }
-    void GetVertexDofNrs (int vnr, Array<int> & dnums) const override
-    { dnums.SetSize(0); }
-    void GetEdgeDofNrs (int ednr, Array<int> & dnums) const override
-    { dnums = IntRange(first_face_dof[ednr],first_face_dof[ednr+1]); }
-    void GetFaceDofNrs (int fanr, Array<int> & dnums) const override
-    { dnums.SetSize(0); }
-    void GetInnerDofNrs (int elnr, Array<int> & dnums) const override
-    { GetDofNrs (elnr, dnums); }
-    std::shared_ptr<ngmg::Prolongation> GetProlongation() const override
-    { return 0; }
-
-  };
-
-}
-
-#endif

netgen/include/hdivfe.hpp

@@ -1,206 +0,0 @@
-#ifndef FILE_HDIVFE
-#define FILE_HDIVFE
-
-/*********************************************************************/
-/* File:   hdivfe.hpp                                                */
-/* Author: Joachim Schoeberl                                         */
-/* Date:   5. Jul. 2001                                              */
-/*********************************************************************/
-
-
-#include "finiteelement.hpp"
-
-namespace ngfem
-{
-
-  /**
-     Finite Elements for H(div)
-     Raviart-Thomas, BDM, BDFM
-  */
-
-
-  class NGS_DLL_HEADER  BaseHDivFiniteElement : public FiniteElement
-  {
-  public:
-    INLINE BaseHDivFiniteElement () { ; } 
-    INLINE BaseHDivFiniteElement (int andof, int aorder)
-      : FiniteElement (andof, aorder) { ; }
-
-    /// compute shape
-    virtual void CalcShape (const IntegrationPoint & ip, 
-			    BareSliceMatrix<> shape) const = 0;
-
-    /// compute div of shape
-    virtual void CalcDivShape (const IntegrationPoint & ip,
-			       BareSliceVector<> divshape) const = 0;
-
-    /// calc normal components of facet shapes, ip has facet-nr
-    virtual void CalcNormalShape (const IntegrationPoint & ip, 
-                                  BareSliceVector<> nshape) const = 0;
-  };
-
-  template <int D>
-  class NGS_DLL_HEADER HDivFiniteElement : public BaseHDivFiniteElement
-  {
-  public:
-    enum { DIM = D };
-
-  public:
-    using BaseHDivFiniteElement::BaseHDivFiniteElement;
-    
-    ///
-    HD virtual ~HDivFiniteElement () { ; }
-
-    /// 
-    virtual string ClassName() const;
-    
-    virtual void CalcDivShape (const IntegrationPoint & ip,
-			       BareSliceVector<> divshape) const;
-    
-    /// calc normal components of facet shapes, ip has facet-nr
-    virtual void CalcNormalShape (const IntegrationPoint & ip, 
-                                  BareSliceVector<> nshape) const;
-
-    /// compute shape
-    virtual void CalcMappedShape (const BaseMappedIntegrationPoint & bmip,
-                                  BareSliceMatrix<> shape) const;
-
-    virtual void CalcMappedShape (const SIMD<MappedIntegrationPoint<DIM,DIM>> & mip,
-				  BareSliceMatrix<SIMD<double>> shape) const;
-
-    virtual void CalcMappedShape (const BaseMappedIntegrationRule & bmir, BareSliceMatrix<> shapes) const;
-
-    virtual void CalcMappedShape (const SIMD_BaseMappedIntegrationRule & mir, 
-                                  BareSliceMatrix<SIMD<double>> shapes) const;
-
-    virtual void CalcMappedNormalShape (const SIMD_BaseMappedIntegrationRule & mir, 
-                                        BareSliceMatrix<SIMD<double>> shapes) const;
-
-    /// compute div of shape
-    virtual void CalcMappedDivShape (const MappedIntegrationPoint<DIM,DIM> & sip,
-				     BareSliceVector<> divshape) const;
-
-    virtual void CalcMappedDivShape (const SIMD_BaseMappedIntegrationRule & mir, 
-                                     BareSliceMatrix<SIMD<double>> divshapes) const;
-
-
-    INLINE const FlatMatrixFixWidth<DIM> GetShape (const IntegrationPoint & ip,
-                                                   LocalHeap & lh) const
-    {
-      FlatMatrixFixWidth<DIM> shape(ndof, lh);
-      CalcShape (ip, shape);
-      return shape;
-    }
-
-    INLINE const FlatVector<> GetDivShape (const IntegrationPoint & ip,
-                                           LocalHeap & lh) const
-    {
-      FlatVector<> divshape(ndof, lh);
-      CalcDivShape (ip, divshape);
-      return divshape;
-    }
-
-    
-    virtual void Evaluate (const IntegrationRule & ir, FlatVector<double> coefs, 
-			   BareSliceMatrix<> vals) const;
-
-    virtual void EvaluateTrans (const IntegrationRule & ir, 
-                                BareSliceMatrix<> vals,
-                                FlatVector<double> coefs) const;
-
-    virtual void Evaluate (const SIMD_BaseMappedIntegrationRule & ir, BareSliceVector<> coefs,
-                           BareSliceMatrix<SIMD<double>> values) const;
-    virtual void AddTrans (const SIMD_BaseMappedIntegrationRule & ir, BareSliceMatrix<SIMD<double>> values,
-                           BareSliceVector<> coefs) const;
-    virtual void EvaluateDiv (const SIMD_BaseMappedIntegrationRule & ir, BareSliceVector<> coefs,
-                              BareVector<SIMD<double>> values) const;
-    virtual void AddDivTrans (const SIMD_BaseMappedIntegrationRule & ir, BareVector<SIMD<double>> values,
-                              BareSliceVector<> coefs) const;
-    
-    virtual void GetFacetDofs(int i, Array<int> & dnums) const;
-    virtual void CalcDualShape (const BaseMappedIntegrationPoint & bmip, BareSliceMatrix<> shape) const;
-
-  protected:
-
-    /// compute basis, will be orthogonalized
-    virtual void CalcShape1 (const IntegrationPoint & ip,
-			     FlatMatrixFixWidth<DIM> shape) const { ; }
-
-    ///
-    virtual void CalcShape2 (const IntegrationPoint & ip,
-			     FlatMatrixFixWidth<DIM> shape) const { ; }
-
-    ///
-    void ComputeFaceMoments (int fnr, ScalarFiniteElement<DIM-1> & testfe,
-			     FlatMatrix<> & moments,
-			     int order, int shape = 1) const;
-
-
-    virtual std::list<std::tuple<std::string,double>> Timing () const;
-  };
-
-
-
-
-  /**
-    HDivNormalFiniteElement
-  */
-
-  template <int D>
-  class NGS_DLL_HEADER HDivNormalFiniteElement : public FiniteElement
-  {
-  public:
-    enum { DIM = D };
-
-  public:
-    ///
-    HDivNormalFiniteElement (int andof, int aorder)
-      : FiniteElement (andof, aorder){;}
-
-    ///
-    HD virtual ~HDivNormalFiniteElement () { ; }
-
-    /// compute shape
-    virtual void CalcShape (const IntegrationPoint & ip,
-			    FlatVector<> shape) const = 0;
-
-    ///
-    const FlatVector<> GetShape (const IntegrationPoint & ip,
-				 LocalHeap & lh) const
-    {
-      FlatVector<> shape(ndof, lh);
-      CalcShape (ip, shape);
-      return shape;
-    }
-
-    virtual void CalcMappedShape (const BaseMappedIntegrationPoint & bmip,
-                                  SliceMatrix<> shape) const
-    {
-      throw Exception(string("CalcMappedShape not implemented for H(div) normal element ")+typeid(*this).name());
-    }
-
-    virtual void Evaluate (const SIMD_BaseMappedIntegrationRule & ir,
-                           BareSliceVector<> coefs,
-                           BareSliceMatrix<SIMD<double>> values) const
-    {
-      throw ExceptionNOSIMD ("HDivNormalFE::Evaluate (simd) not overloaded");
-    }
-    
-    // values * normal
-    virtual void AddTrans (const SIMD_BaseMappedIntegrationRule & ir,
-                           BareSliceMatrix<SIMD<double>> values,
-                           BareSliceVector<> coefs) const
-    {
-      throw ExceptionNOSIMD ("HDivNormalFE::AddTrans (simd) not overloaded");
-    }
-    
-  };
-
-
-
-}
-
-#endif
-
-
-