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

netgen/include/symmetricmatrix.hpp

@@ -1,144 +0,0 @@
-#ifndef FILE_SYMMETRICMATRIX
-#define FILE_SYMMETRICMATRIX
-
-/****************************************************************************/
-/* File:   symmetricmatrix.hpp                                              */
-/* Author: Joachim Schoeberl                                                */
-/* Date:   12. Dec. 2004                                                    */
-/****************************************************************************/
-
-
-namespace ngbla
-{
-
-
-
-  /// A symmetric band-matrix
-  template <class T>
-  class FlatSymmetricMatrix
-  {
-  protected:
-    ///
-    int n;
-    ///
-    T *data;
-  public:
-    typedef typename mat_traits<T>::TV_COL TV;
-
-    ///
-    FlatSymmetricMatrix () { ; }
-
-    ///
-    FlatSymmetricMatrix (int an, T * adata)
-      : n(an), data(adata)
-    { ; }
-
-    /// set size, and assign mem
-    void AssignMemory (int an, T * mem) throw()
-    {
-      n = an;
-      data = mem;
-    }
-  
-
-
-    ///
-    void Mult (const FlatVector<TV> & x, FlatVector<TV> & y) const
-    {
-      for (int i = 0; i < n; i++)
-        {
-          TV xi = x(i);
-          TV sum = (*this)(i,i) * xi;
-          for (int j = 0; j < i; j++)
-            {
-              sum += (*this)(i,j) * x(j);
-              y(j) += Trans((*this)(i,j)) * xi;
-            }
-          y(i) = sum;
-        }
-    }
-
-    ///
-    ostream & Print (ostream & ost) const
-    {
-      for (int i = 0; i < n; i++)
-        {
-          for (int j = 0; j < n; j++)
-            if (j <= i)
-              ost << setw(8) << (*this)(i,j) << " ";
-            else 
-              ost << setw(8) << "sym" << " ";
-          ost << endl;
-        }
-      return *this;
-    }
-
-    ///
-    int Height() const { return n; }
-
-  
-    ///
-    const T & operator() (int i, int j) const
-    { return data[ i * (i+1) / 2 + j ]; }
-
-    ///
-    T & operator() (int i, int j) 
-    { return data[ i * (i+1) / 2 + j ]; }
-
-  
-    ///
-    FlatSymmetricMatrix & operator= (const T & val)
-    {
-      int nel = n * (n+1) / 2;
-      for (int i = 0; i < nel; i++)
-        data[i] = val;
-      return *this;
-    }
-  };
-
-
-  template<typename T>
-  inline std::ostream & operator<< (std::ostream & s, const FlatSymmetricMatrix<T> & m)
-  {
-    m.Print (s);
-    return s;
-  }
-
-
-
-  ///  A symmetric band-matrix with memory management
-  template <class T>
-  class SymmetricMatrix : public FlatSymmetricMatrix<T>
-  {
-  public:
-    typedef typename mat_traits<T>::TV_COL TV;
-
-    ///
-    SymmetricMatrix (int an)
-      : FlatSymmetricMatrix<T> (an, new T[an*(an+1)/2])
-    { ; }
-
-    ///
-    ~SymmetricMatrix ()
-    { delete [] this->data; }
-
-    ///
-    SymmetricMatrix & operator= (const T & val)
-    {
-      int nel = this->n * (this->n+1) / 2;
-      for (int i = 0; i < nel; i++)
-        this->data[i] = val;
-      return *this;
-    }
-  };
-
-
-
-
-
-
-
-}
-
-
-#endif

netgen/include/tangentialfacetfe.hpp

@@ -1,224 +0,0 @@
-#ifndef VECTOR_FACET_FE_HPP__
-#define VECTOR_FACET_FE_HPP__
-
-/*********************************************************************/
-/* File:   vectorfacetfe.hpp                                         */
-/* Author: A. Sinwel, (J. Schoeberl)                                 */
-/* Date:   2008                                                      */
-/*********************************************************************/
-
-
-// #include <fem.hpp>
-#include "hcurlfe.hpp"
-#include "recursive_pol.hpp"
-#include <cassert>
-
-namespace ngfem 
-{
-  /*
-    facet element with tangential facet components.
-    to be changed similar to scalar facetfe
-  */
-  
-  template <ELEMENT_TYPE ET>
-  class TangentialFacetFacetFE :
-    public HCurlFiniteElement<ET_trait<ET>::DIM>,
-    public VertexOrientedFE<ET>,
-    public ET_trait<ET>
-  {
-  protected:
-    IVec<2> order_inner;
-    using VertexOrientedFE<ET>::vnums;
-    using ET_trait<ET>::DIM;
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::order;
- 
-  public:
-    using VertexOrientedFE<ET>::SetVertexNumber;
-    using VertexOrientedFE<ET>::SetVertexNumbers;
-
-    TangentialFacetFacetFE (int aorder)
-    {
-      order = aorder;
-      order_inner = IVec<2>(aorder,aorder);
-      ComputeNDof();
-    }
-
-    TangentialFacetFacetFE () { ; }
-
-    HD virtual ELEMENT_TYPE ElementType() const override { return ELEMENT_TYPE(ET); }
-
-    INLINE void SetOrder (int aorder)
-    {
-      order = aorder;
-      order_inner = IVec<2>(aorder,aorder);
-      ComputeNDof();
-    }
-  
-    INLINE void SetOrder (IVec<2> oi)
-    {
-      order = max2 (oi[0], oi[1]);
-      order_inner = oi;
-      ComputeNDof();
-    }
-
-    virtual void ComputeNDof ();
-
-    virtual void CalcShape(const IntegrationPoint & ip,
-                           BareSliceMatrix<> shape) const override;
-
-    template<typename Tx, typename TFA>  
-    void T_CalcShape (TIP<DIM,Tx> tip, TFA & shape) const;
-  };
-
-
-  template <ELEMENT_TYPE ET>
-  class TangentialFacetVolumeFE : public HCurlFiniteElement<ET_trait<ET>::DIM>, public VertexOrientedFE<ET>
-  {
-  protected:
-    using ET_T = ET_trait<ET>;
-    IVec<2> facet_order[ET_T::N_FACET];
-    int first_facet_dof[ET_T::N_FACET+1];
-    bool highest_order_dc;
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::order;
-    using VertexOrientedFE<ET>::vnums;
-    enum { DIM = ET_trait<ET>::DIM };
-    
-  public:
-    using VertexOrientedFE<ET>::SetVertexNumbers;
-    
-    TangentialFacetVolumeFE () { highest_order_dc=false; }
-    
-    HD virtual ELEMENT_TYPE ElementType() const override { return ELEMENT_TYPE(ET); }
-
-    void SetHighestOrderDC(bool set) { highest_order_dc=set; }
-
-    void SetOrder(int ao)
-    {
-      order = ao;
-      for ( int i = 0; i < ET_T::N_FACET; i++ )
-        facet_order[i] = IVec<2> (ao, ao);
-      ComputeNDof();
-    }
-
-    void SetOrder(FlatArray<int> & ao)
-    {
-      order = 0;
-      assert(ao.Size()==ET_T::N_FACET);
-      for ( int i = 0; i < ET_T::N_FACET; i++ )
-        {
-          order = max2 ( order, ao[i] );
-          facet_order[i] = IVec<2> (ao[i], ao[i]);
-        }
-      ComputeNDof();
-    }
-
-    void SetOrder(FlatArray<IVec<2> > & ao)
-    {
-      order = 0;
-      assert(ao.Size()==ET_T::N_FACET);
-      for ( int i = 0; i < ET_T::N_FACET; i++ )
-        {
-          order = max3 ( order, ao[i][0], ao[i][1] );
-          facet_order[i] = ao[i];
-        }
-      ComputeNDof();
-    }
-
-    IVec<2> GetFacetOrder(int j) const { return facet_order[j]; }
-    int GetVertexNumber(int j) const { return vnums[j]; }
-    
-    virtual void CalcShape (const IntegrationPoint & ip, BareSliceMatrix<> shape) const override
-    {
-      int fnr = ip.FacetNr();
-      if (fnr >= 0)
-        {
-          CalcShape (ip, fnr, shape);
-          return;
-        }
-      throw Exception("TangentialFacetVolumeFiniteElement<D>::CalcShape in global coordinates disabled");
-    }
-
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::CalcMappedShape;
-    virtual void CalcMappedShape (const SIMD_BaseMappedIntegrationRule & mir, 
-				  BareSliceMatrix<SIMD<double>> shapes) const override;
-
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::Evaluate;
-    virtual void Evaluate (const SIMD_BaseMappedIntegrationRule & ir, BareSliceVector<> coefs,
-                           BareSliceMatrix<SIMD<double>> values) const override;
-    
-    virtual void AddTrans (const SIMD_BaseMappedIntegrationRule & ir, BareSliceMatrix<SIMD<double>> values,
-                           BareSliceVector<> coefs) const override;
-
-    template<typename Tx, typename TFA>  
-    void T_CalcShape (Tx hx[DIM], int fnr, TFA & shape) const;
-
-    template<typename MIP, typename TFA>  
-    void CalcDualShape2 (const MIP & mip, int fnr, TFA & shape) const;
-    
-    virtual void CalcShape (const IntegrationPoint & ip, int facet, BareSliceMatrix<> shape) const;
-
-    virtual int GetNExtraShapes( int facet) const {return 0;}
-    virtual void CalcExtraShape (const IntegrationPoint & ip, int facet, FlatMatrixFixWidth<ET_T::DIM> xshape) const {xshape = 0.0;}
-
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::ndof;
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::CalcDualShape;
-    virtual void CalcDualShape (const BaseMappedIntegrationPoint & bmip, BareSliceMatrix<> shape) const override;
-    virtual void CalcDualShape (const SIMD_BaseMappedIntegrationRule & bmir, BareSliceMatrix<SIMD<double>> shape) const override;
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::EvaluateDual;
-    virtual void EvaluateDual (const SIMD_BaseMappedIntegrationRule & bmir, BareSliceVector<> coefs, BareSliceMatrix<SIMD<double>> values) const override;
-    using HCurlFiniteElement<ET_trait<ET>::DIM>::AddDualTrans;
-    virtual void AddDualTrans (const SIMD_BaseMappedIntegrationRule & bmir, BareSliceMatrix<SIMD<double>> values,
-                               BareSliceVector<double> coefs) const override;
-
-    virtual void GetFacetDofNrs(int afnr, Array<int>& fdnums) const
-    {
-      int first = first_facet_dof[afnr];
-      int next = first_facet_dof[afnr+1];
-      assert(next > first);
-      fdnums.SetSize(next-first);
-      for (auto i : Range(next-first))
-        fdnums[i] = first+i;
-    }
-
-    virtual int GetFacetNDof(int afnr) const
-    { return first_facet_dof[afnr+1] - first_facet_dof[afnr]; };
-    virtual int GetFirstFacetDof(int afnr) const {
-      return first_facet_dof[afnr];};
-
-    virtual void ComputeNDof ();
-
-    /// degrees of freedom sitting inside the element, used for static condensation
-    virtual void GetInternalDofs (Array<int> & idofs) const
-    {
-      idofs.SetSize(0);
-      if (highest_order_dc)
-        {
-          if (ET_T::DIM==2)
-            {
-              for (int i=0; i < ET_T::N_FACET; i++)
-                idofs.Append(first_facet_dof[i+1]-1);
-            }
-          else
-            {
-              for (int i=0; i < ET_T::N_FACET; i++)
-                {
-                  int pos = first_facet_dof[i]-2;
-                  int fac = 4 - ElementTopology::GetNVertices(ElementTopology::GetFaceType(ET,i));
-                  for (int k = 0; k <= facet_order[i][0]; k++){
-                    pos += 2*(facet_order[i][0]+1-fac*k);
-                    idofs.Append(pos);
-                    idofs.Append(pos+1);
-                  }
-                }
-            }//end if Dimension
-        }
-    }
-  };
-
-}
-
-
-
-#endif
-
-

netgen/include/tangentialfacetfespace.hpp

@@ -1,106 +0,0 @@
-#ifndef TANGENTIAL_FACET_FESPACE_HPP
-#define TANGENTIAL_FACET_FESPACE_HPP
-
-/*********************************************************************/
-/* File:   tangentialfacetfespace.hpp                                */
-/* Author: A. Sinwel, Joachim Schoeberl                              */
-/* Date:   2008                                                      */
-/*********************************************************************/
-
-
-#include "fespace.hpp"
-
-namespace ngcomp
-{
-
-  class NGS_DLL_HEADER TangentialFacetFESpace : public FESpace
-  {
-  protected:
-    /// Level
-    // int level;
-    /// Number of Facets
-    // int nfacets;
-    /// 
-    // int ncfacets;
-    ///
-    // int nel;
-
-    Array<int> first_facet_dof;
-    Array<int> first_inner_dof;  // for highest_order_dc
-    // int ndof_lo;
-
-    int rel_order;
-
-    Array<IVec<2> > order_facet;
-    Array<bool> fine_facet;
-
-    // int ndof;
-    // Array<int> ndlevel;
-    bool var_order;
-    bool print;
-
-    bool highest_order_dc;
-    bool hide_highest_order_dc;
-    
-  public:
-    ///
-    TangentialFacetFESpace (shared_ptr<MeshAccess> ama, const Flags & flags, 
-                            bool parseflags = false );
-
-    virtual ~TangentialFacetFESpace () { ; }
-    static DocInfo GetDocu ();
-
-    virtual string GetClassName () const override
-    {
-      return "TangentialFacetFESpace";
-    }
-
-    void Update() override;
-    virtual void UpdateCouplingDofArray() override;
-
-    virtual void SetOrder (NodeId ni, int order) override;
-    virtual int GetOrder (NodeId ni) const override;
-
-    
-    // virtual size_t GetNDof() const throw() override { return ndof; }
-    // virtual size_t GetNDofLevel ( int i ) const override { return ndlevel[i]; }
-
-    virtual FlatArray<VorB> GetDualShapeNodes (VorB vb) const override;
-
-    // virtual int GetNDofLowOrder () const
-    // { return ndof_lo; }
-
-    virtual FiniteElement & GetFE(ElementId ei, Allocator & lh) const override;
-    
-    // virtual const FiniteElement & GetFE ( int elnr, LocalHeap & lh ) const;
-    // virtual const FiniteElement & GetSFE ( int selnr, LocalHeap & lh ) const;
-
-    virtual void GetFacetDofNrs (int felnr, Array<DofId> & dnums) const;
-
-    virtual int GetNFacetDofs (int felnr) const;
-
-    virtual void GetDofNrs (ElementId ei, Array<DofId> & dnums) const override;
-    ///
-    virtual shared_ptr<Array<int>> CreateDirectSolverClusters (const Flags & precflags) const override;
-  
-    // some utility functions for convenience
-    ///
-    // virtual void GetVertexNumbers(int elnr, Array<int>& vnums) const;
-    ///
-    virtual IVec<2> GetFacetOrder(int fnr) const;
-
-    virtual int GetFirstFacetDof(int fanr) const;
-
-    virtual bool UsesHighestOrderDiscontinuous() const {return highest_order_dc;};
-
-    virtual void GetVertexDofNrs (int elnum, Array<DofId> & dnums) const override;
-    virtual void GetEdgeDofNrs (int elnum, Array<DofId> & dnums) const override;
-    virtual void GetFaceDofNrs (int felnr, Array<DofId> & dnums) const override;
-    virtual void GetInnerDofNrs (int felnr, Array<DofId> & dnums) const override;
-
-  };
-
-}
-
-#endif
-