ngsolve 6.2.2406__cp311-cp311-macosx_10_15_universal2.whl → 6.2.2406.post111.dev1__cp311-cp311-macosx_10_15_universal2.whl

netgen/include/fespace.hpp

@@ -12,7 +12,12 @@
 #include <finiteelement.hpp>
 #include <diffop.hpp>
 #include <symbolicintegrator.hpp>   // for ProxyFunction
-#include <la.hpp>
+
+#include <basevector.hpp>
+#include <basematrix.hpp>
+
+// #include <paralleldofs.hpp>
+
 
 #include "ngsobject.hpp"
 
@@ -283,14 +288,7 @@ ANY                  1 1 1 1 | 15
       order_policy = op;
     }
     
-    virtual void SetOrder (ELEMENT_TYPE et, TORDER order)
-    {
-      if (order_policy == CONSTANT_ORDER || order_policy == OLDSTYLE_ORDER)
-        order_policy = NODE_TYPE_ORDER;
-      et_bonus_order[et] = order - this->order;
-
-      timestamp = 0;  // rerun first_update
-    }
+    virtual void SetOrder (ELEMENT_TYPE et, TORDER order);
 
     virtual void SetOrder (NodeId ni, int order); 
     virtual int GetOrder (NodeId ni) const; 
@@ -924,7 +922,7 @@ ANY                  1 1 1 1 | 15
   class NGS_DLL_HEADER NonconformingFESpace : public FESpace
   {
     ///
-    Array<int> ndlevel;
+    // Array<int> ndlevel;
 
   public:
     NonconformingFESpace (shared_ptr<MeshAccess> ama, const Flags & flags, bool parseflags=false);
@@ -938,7 +936,7 @@ ANY                  1 1 1 1 | 15
 
     virtual FiniteElement & GetFE (ElementId ei, Allocator & lh) const override;
     ///
-    virtual size_t GetNDof () const throw() override;
+    // virtual size_t GetNDof () const throw() override;
     ///
     virtual void GetDofNrs (ElementId ei, Array<DofId> & dnums) const override;
   };
@@ -949,7 +947,7 @@ ANY                  1 1 1 1 | 15
   class NGS_DLL_HEADER NonconformingSurfaceFESpace : public FESpace
   {
     ///
-    Array<int> ndlevel;
+    // Array<int> ndlevel;
 
   public:
     NonconformingSurfaceFESpace (shared_ptr<MeshAccess> ama, const Flags & flags, bool parseflags=false);
@@ -963,7 +961,7 @@ ANY                  1 1 1 1 | 15
 
     virtual FiniteElement & GetFE (ElementId ei, Allocator & lh) const override;
     ///
-    virtual size_t GetNDof () const throw() override;
+    // virtual size_t GetNDof () const throw() override;
     ///
     virtual void GetDofNrs (ElementId ei, Array<DofId> & dnums) const override;
   };
@@ -1433,7 +1431,7 @@ namespace ngcore
   template<>
   struct MPI_typetrait<ngcomp::COUPLING_TYPE>
   {
-    static NG_MPI_Datatype MPIType () 
+    static auto MPIType () 
     { 
       static_assert ( (sizeof(ngcomp::COUPLING_TYPE) == sizeof(char)) ||
                       (sizeof(ngcomp::COUPLING_TYPE) == sizeof(int)) );

netgen/include/finiteelement.hpp

@@ -71,7 +71,7 @@ namespace ngfem
     /// the name of the element family
     virtual string ClassName() const;
 
-    virtual void SetVertexNumbers (FlatArray<int> vnums);
+    virtual FiniteElement * SetVertexNumbers (FlatArray<int> vnums);
     
     virtual IntegrationRule GetIR (int order) const;
 
@@ -103,13 +103,12 @@ namespace ngfem
   {
   protected:
     /// pointers to the components
-    // ArrayMem<const FiniteElement*,10> fea;
-    FlatArray<const FiniteElement*> fea;
+    FlatArray<FiniteElement*> fea;
     bool all_the_same{true};
 
   public:
     /// initialize with pointers to components, copy pointers
-    CompoundFiniteElement (FlatArray<const FiniteElement*> afea);
+    CompoundFiniteElement (FlatArray<FiniteElement*> afea);
 
     HD virtual ELEMENT_TYPE ElementType() const override { return fea[0]->ElementType(); }
     /// number of components
@@ -130,10 +129,11 @@ namespace ngfem
     /// the name of the element family
     virtual string ClassName() const override { return "CompoundFiniteElement"; }
     
-    virtual void SetVertexNumbers (FlatArray<int> vnums) override
+    virtual CompoundFiniteElement * SetVertexNumbers (FlatArray<int> vnums) override
     {
       for (auto pfel : fea)
-        const_cast<FiniteElement*>(pfel) -> SetVertexNumbers(vnums);
+        pfel -> SetVertexNumbers(vnums);
+      return this;
     }
 
     HD virtual bool ComplexShapes() const override;
@@ -174,7 +174,7 @@ namespace ngfem
     /// the name of the element family
     virtual string ClassName() const override { return "VectorFiniteElement"; }
 
-    virtual void SetVertexNumbers (FlatArray<int> vnums) override;
+    virtual VectorFiniteElement * SetVertexNumbers (FlatArray<int> vnums) override;
 
     virtual void Print (ostream & ost) const override;
 

netgen/include/h1amg.hpp

@@ -1,7 +1,9 @@
 #ifndef H1AMG_HPP_
 #define H1AMG_HPP_
 
-#include <la.hpp>
+// #include <la.hpp>
+#include <basematrix.hpp>
+#include <sparsematrix.hpp>
 
 namespace ngcomp
 {

netgen/include/hcurlhdivfes.hpp

@@ -9,6 +9,8 @@
 
 
 #include "fespace.hpp"
+#include <sparsematrix.hpp>
+
 
 namespace ngcomp
 {

netgen/include/hcurlhofespace.hpp

@@ -9,6 +9,7 @@
 
 #include "fespace.hpp"
 #include "h1hofespace.hpp"
+#include <sparsematrix.hpp>
 
 namespace ngcomp
 {
@@ -23,7 +24,7 @@ namespace ngcomp
   protected:
 
     typedef short TORDER;
-
+    
     // Level
     int level;
     Array<DofId> first_edge_dof;
@@ -58,7 +59,7 @@ namespace ngcomp
     bool usegrad;  
     bool var_order; 
   
-    int ndof;
+    // int ndof;
     int nedfine; 
     int uniform_order_inner;
     int uniform_order_face; 
@@ -95,7 +96,7 @@ namespace ngcomp
     ///
     virtual void DoArchive (Archive & archive) override;
     ///
-    virtual size_t GetNDof () const throw() override;
+    // virtual size_t GetNDof () const throw() override;
     virtual void SetOrder (NodeId ni, int order) override;
     virtual int GetOrder (NodeId ni) const override;
     using FESpace::GetOrder;

netgen/include/hdivfes.hpp

@@ -19,7 +19,7 @@ namespace ngcomp
   class RaviartThomasFESpace : public FESpace
   {
     ///
-    Array<int> ndlevel;
+    // Array<int> ndlevel;
   public:
     ///
     /*
@@ -42,9 +42,9 @@ namespace ngcomp
 
     virtual FiniteElement & GetFE (ElementId ei, Allocator & lh) const override;
     ///
-    virtual size_t GetNDof () const throw() override;
+    // virtual size_t GetNDof () const throw() override;
     ///
-    virtual size_t GetNDofLevel (int level) const override;
+    // virtual size_t GetNDofLevel (int level) const override;
 
     ///
     virtual void GetDofNrs (ElementId ei, Array<DofId> & dnums) const override;

netgen/include/hdivhosurfacefespace.hpp

@@ -15,7 +15,7 @@ namespace ngcomp
   {
   protected:
     
-    int ndof;
+    // int ndof;
     
     Array<DofId> first_facet_dof;
     Array<DofId> first_inner_dof;
@@ -28,7 +28,7 @@ namespace ngcomp
     Array<bool> fine_facet; 
     Array<bool> boundary_facet; 
  
-    Array<int> ndlevel;
+    // Array<int> ndlevel;
     int uniform_order_inner; 
     int uniform_order_facet; 
 
@@ -64,7 +64,7 @@ namespace ngcomp
     virtual void UpdateDofTables() override; 
     virtual void UpdateCouplingDofArray() override;   
     
-    virtual size_t GetNDof () const throw() override;
+    // virtual size_t GetNDof () const throw() override;
     
     virtual FiniteElement & GetFE (ElementId ei, Allocator & alloc) const override;            
 

netgen/include/intrule.hpp

@@ -436,7 +436,7 @@ namespace ngfem
     INLINE SCAL operator() (int i) const { return point(i); }
 
     ///
-    INLINE const Vec<R,SCAL> GetNV () const { return normalvec; }
+    INLINE const Vec<R,SCAL> & GetNV () const { return normalvec; }
     /// 
     INLINE void SetNV ( Vec<R,SCAL> vec) { normalvec = vec; }
     ///
@@ -1481,6 +1481,7 @@ namespace ngfem
     int DimSpace() const;
     
     virtual SliceMatrix<> GetPoints() const = 0;
+    virtual SliceMatrix<> GetNormals() const = 0;    
     virtual SliceMatrix<Complex> GetPointsComplex() const
     { throw Exception("don't have complex ir"); }
     virtual void ComputeNormalsAndMeasure (ELEMENT_TYPE et) { throw Exception ("ComputeNormalsAndMeasure(ET) not overloaded"); }
@@ -1551,6 +1552,13 @@ namespace ngfem
                             const_cast<double*> (&mips[0].GetPoint()(0)));
     }
 
+    virtual SliceMatrix<> GetNormals() const
+    {
+      return SliceMatrix<> (mips.Size(), DIM_SPACE*sizeof(SCAL)/sizeof(double),
+                            sizeof(MappedIntegrationPoint<DIM_ELEMENT, DIM_SPACE, SCAL>) / sizeof(double),
+                            const_cast<double*> (&mips[0].GetNV()(0)));
+    }
+
     virtual void ComputeNormalsAndMeasure (ELEMENT_TYPE et);
     virtual void ComputeNormalsAndMeasure (ELEMENT_TYPE et, int facetnr);
     virtual bool IsComplex() const { return false; } 
@@ -1616,6 +1624,12 @@ namespace ngfem
                                    const_cast<Complex*> (&mips[0].GetPointComplex()(0)));
     }
 
+    virtual SliceMatrix<> GetNormals() const
+    {
+      throw Exception("never tested");
+    }
+    
+
     virtual void ComputeNormalsAndMeasure (ELEMENT_TYPE et, int facetnr);
     virtual bool IsComplex() const { return true; }     
   };

netgen/include/jacobi.hpp

@@ -7,6 +7,9 @@
 /* Date:   06. Oct. 96                                                    */
 /* *************************************************************************/
 
+
+#include "sparsematrix.hpp"
+
 namespace ngla
 {
 

netgen/include/la.hpp

@@ -2,15 +2,6 @@
 #define FILE_NGLA
 
 #include <bla.hpp>
-#include <core/mpi_wrapper.hpp>
-
-/*
-namespace ngcomp
-{
-  class Preconditioner;
-  class LocalPreconditioner;
-}
-*/
 
 /** namespace for linear algebra.
  */
@@ -33,13 +24,10 @@ namespace ngla
 #include "order.hpp"
 #include "sparsecholesky.hpp"
 #include "pardisoinverse.hpp"
-// include these only from c++-files
-// #include "umfpackinverse.hpp"
-// #include "superluinverse.hpp"
-// #include "mumpsinverse.hpp"
 #include "jacobi.hpp"
 #include "blockjacobi.hpp"
 #include "commutingAMG.hpp"
+#include "diagonalmatrix.hpp"
 #include "special_matrix.hpp"
 #include "elementbyelement.hpp"
 #include "cg.hpp"

netgen/include/matrix.hpp

@@ -243,11 +243,21 @@ namespace ngbla
       return FlatMatrix (next-first, w, data+first*w);
     }
 
+    INLINE FlatMatrix Rows (size_t next) const
+    {
+      return FlatMatrix (next, w, data);
+    }
+
     INLINE SliceMatrix<T> Cols (size_t first, size_t next) const
     {
       return SliceMatrix<T> (h, next-first, w, data+first);
     }
 
+    INLINE SliceMatrix<T> Cols (size_t next) const
+    {
+      return SliceMatrix<T> (h, next, w, data);
+    }
+
     INLINE FlatMatrix Rows (IntRange range) const
     {
       return FlatMatrix (range.Next()-range.First(), w, data+range.First()*w);
@@ -1356,6 +1366,12 @@ namespace ngbla
         return MatrixView<T,ORD,TH,size_t,TDIST>  (h, next-first, Dist(), Addr(0,first));
     }
 
+    INLINE auto Cols (size_t next) const
+    {
+      return Cols(0, next);
+    }
+
+    
     INLINE auto Rows (size_t first, size_t next) const
     {
       if constexpr (ORD==RowMajor)
@@ -1369,6 +1385,11 @@ namespace ngbla
       return Rows (range.First(), range.Next());
     }
 
+    INLINE auto Rows (size_t next) const
+    {
+      return Rows (0, next);
+    }
+
     template <int R>
     INLINE auto Rows() const
     {

netgen/include/meshaccess.hpp

@@ -9,8 +9,6 @@
 
 
 #include <nginterface_v2.hpp>
-#include <core/ranges.hpp>
-
 #include <elementtopology.hpp>
 
 namespace ngfem
@@ -26,7 +24,6 @@ namespace ngcomp
 {
   class PML_Transformation;
   
-  // using ngcore::INT;
   using netgen::Ng_Node;
   using ngfem::ELEMENT_TYPE;
   
@@ -158,30 +155,6 @@ namespace ngcomp
     INLINE auto end () const { return DimElementIterator<VB,DIM>(ma, r.Next()); }
   };
 
-  /*
-  class NodeIterator
-  {
-    NodeId ni;
-  public:
-    NodeIterator (NodeId ani) : ni(ani) { ; }
-    NodeIterator operator++ () { return NodeIterator(++ni); }
-    INLINE NodeId operator*() const { return ni; }
-    bool operator!=(NodeIterator id2) const { return ni != id2.ni; }
-    bool operator==(NodeIterator id2) const { return ni == id2.ni; }
-  };
-
-  class NodeRange : public IntRange
-  {
-    NODE_TYPE nt;
-  public:
-    NodeRange (NODE_TYPE ant, IntRange ar) 
-      : IntRange(ar), nt(ant) { ; } 
-    NodeIterator begin () const { return NodeIterator(NodeId(nt,IntRange::First())); }
-    NodeIterator end () const { return NodeIterator(NodeId(nt,IntRange::Next())); }
-    NodeId operator[] (size_t nr) { return NodeId(nt, IntRange::First()+nr); }
-  };
-  */
-
   /** 
       Access to mesh topology and geometry.
 
@@ -240,9 +213,12 @@ namespace ngcomp
 
     /// store periodic vertex mapping for each identification number
     // shared ptr because Meshaccess is copy constructible
+    /*
     shared_ptr<Array<Array<IVec<2>>>> periodic_node_pairs[3] = {make_shared<Array<Array<IVec<2>>>>(),
                                                                make_shared<Array<Array<IVec<2>>>>(),
                                                                make_shared<Array<Array<IVec<2>>>>()};
+    */
+    Array<Array<IVec<2>>> periodic_node_pairs[3];
 
     DynamicTable<size_t> neighbours[4][4];
     friend class Region;
@@ -662,7 +638,7 @@ namespace ngcomp
     
     void SetDeformation (shared_ptr<GridFunction> def = nullptr);
 
-    const shared_ptr<GridFunction> & GetDeformation () const
+    shared_ptr<GridFunction> GetDeformation () const
     {
       return deformation;
     }
@@ -817,15 +793,6 @@ namespace ngcomp
     void GetFaceEdges (int fnr, Array<int> & edges) const;
     INLINE auto GetFaceEdges (size_t fnr) const
     { return ArrayObject(mesh.GetFaceEdges(fnr)); }
-    /*
-    {
-      ArrayMem<int,4> f2ed;
-      GetFaceEdges (fnr, f2ed);
-      f2ed.NothingToDelete(); // dynamic allocation never needed
-      return f2ed;
-    }
-    */
-    
     void GetEdgeFaces (int enr, Array<int> & faces) const;
     /// returns elements connected to a face
     void GetFaceElements (int fnr, Array<int> & elnums) const;
@@ -911,7 +878,6 @@ namespace ngcomp
     }
 
 
-    // void GetVertexElements (int vnr, Array<int> & elnrs) const;
     /// element order stored in Netgen
     int GetElOrder (int enr) const
     { return mesh.GetElementOrder (enr+1); } 
@@ -1053,11 +1019,11 @@ namespace ngcomp
 	}
     }
     
-
+    /*
     // (old style optimization)
     [[deprecated("functionality not useful anymore, just remove function call!")]]            
     void SetPointSearchStartElement(const int el) const;
-
+    */
 
     
     ElementId FindElementOfPoint (FlatVector<double> point,
@@ -1083,6 +1049,7 @@ namespace ngcomp
     { return GetElement(ElementId(VOL,elnr)).is_curved; }
       // { return bool (Ng_IsElementCurved (elnr+1)); }
 
+    /*
     [[deprecated("Use GetPeriodicNodes(NT_VERTEX, pairs) instead!")]]
     void GetPeriodicVertices ( Array<IVec<2> > & pairs) const;
     [[deprecated("Use GetNPeriodicNodes(NT_VERTEX) instead!")]]
@@ -1100,10 +1067,11 @@ namespace ngcomp
     void GetPeriodicEdges (int idnr, Array<IVec<2> > & pairs) const;
     [[deprecated("Use GetPeriodicNodes(NT_EDGE, idnr).Size() instead")]]
     int GetNPairsPeriodicEdges (int idnr) const;
-
+    */
+    
     int GetNPeriodicIdentifications() const
     {
-      return periodic_node_pairs[NT_VERTEX]->Size();
+      return periodic_node_pairs[NT_VERTEX].Size();
     }
     // get number of all periodic nodes of nodetype nt
     size_t GetNPeriodicNodes(NODE_TYPE nt) const;
@@ -1147,9 +1115,10 @@ namespace ngcomp
 
     // void PrecomputeGeometryData(int intorder);
 
+    /*
     void InitPointCurve(double red = 1, double green = 0, double blue = 0) const;
     void AddPointCurvePoint(const Vec<3> & point) const;
-
+    */
 
     
     template <int DIMS, int DIMR> friend class Ng_ElementTransformation;
@@ -1169,10 +1138,11 @@ namespace ngcomp
        Returns the global number of the node.
        Currently, this function works only for vertex-nodes.
      */
-
+    /*
     [[deprecated("should not need global numbers")]]                    
     size_t GetGlobalNodeNum (NodeId node) const;
-
+    */
+    
     size_t GetGlobalVertexNum (int locnum) const;
     
     FlatArray<int> GetDistantProcs (NodeId node) const
@@ -1357,21 +1327,21 @@ namespace ngcomp
       for (auto p : GetDistantProcs(Node(nt, i)))
         dist_data[p][cnt[p]++] = data[i];
 
-    Array<NG_MPI_Request> requests;
+    NgMPI_Requests requests;
     for (auto i : cnt.Range())
       if (cnt[i])
 	{
-	  requests.Append (comm.ISend(dist_data[i], i, NG_MPI_TAG_SOLVE));
-	  requests.Append (comm.IRecv(recv_data[i], i, NG_MPI_TAG_SOLVE));
+	  requests += comm.ISend(dist_data[i], i, NG_MPI_TAG_SOLVE);
+	  requests += comm.IRecv(recv_data[i], i, NG_MPI_TAG_SOLVE);
 	}
-    MyMPI_WaitAll (requests);
+    requests.WaitAll();
     
     cnt = 0;
-    NG_MPI_Datatype type = GetMPIType<T>();
+    auto type = GetMPIType<T>();
     for (auto i : Range(GetNNodes(nt)))
       for (auto p : GetDistantProcs(Node(nt, i)))
         NG_MPI_Reduce_local (&recv_data[p][cnt[p]++],
-                          &data[i], 1, type, op);
+                             &data[i], 1, type, op);
   }
 
 

netgen/include/mgpre.hpp

@@ -9,7 +9,7 @@
 /*********************************************************************/
 
 #include <bilinearform.hpp>
-#include <la.hpp>
+// #include <la.hpp>
 
 namespace ngmg
 {

netgen/include/mptools.hpp

@@ -813,7 +813,6 @@ c
       scale = newscale;
     }
 
-    
     template <typename TARGET>
     void Transform (MultiPole<TARGET> & target, Vec<3> dist) const
     {
@@ -2095,12 +2094,12 @@ c
     template <typename TARGET>
     void ShiftZ (double z, MultiPole<TARGET> & target) { mp.ShiftZ(z, target); }
 
+    using CoefficientFunction::Transform;        
     template <typename TARGET>
     void Transform (MultiPoleCF<TARGET> & target)
     {
       mp.Transform (target.MP(), target.Center()-center);
     }
-    
   };
 
   

netgen/include/multivector.hpp

@@ -7,6 +7,8 @@
 /* Date:   June 2020                                                 */
 /*********************************************************************/
 
+#include "basevector.hpp"
+
 namespace ngla {
 
   /* TODO: 

netgen/include/ngsobject.hpp

@@ -2,7 +2,7 @@
 #define FILE_NGSOBJECT
 
 /*********************************************************************/
-/* File:   ngsobject.hh                                              */
+/* File:   ngsobject.hpp                                             */
 /* Author: Joachim Schoeberl                                         */
 /* Date:   2. Aug. 2000                                              */
 /*********************************************************************/

netgen/include/nodalhofe.hpp

@@ -29,8 +29,11 @@ namespace ngfem
     using ET_trait<ET>::PolDimension;
     using ET_trait<ET>::PolBubbleDimension;
 
+    
   public:
     using VertexOrientedFE<ET>::SetVertexNumbers;    
+    NodalHOFE * SetVertexNumbers (FlatArray<int> vnums) override
+    { VertexOrientedFE<ELEMENT_TYPE(ET)>::SetVertexNumbers(vnums); return this; }   // cast for msvc ?
     using ET_trait<ET>::ElementType;
 
     /// builds a functional element of order aorder.

netgen/include/normalfacetfespace.hpp

@@ -25,8 +25,8 @@ namespace ngcomp
     Array<IVec<2> > order_facet;
     Array<bool> fine_facet;
 
-    int ndof;
-    Array<int> ndlevel;
+    // int ndof;
+    // Array<int> ndlevel;
     bool var_order;
     bool print;
 
@@ -51,9 +51,8 @@ namespace ngcomp
     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 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;
 

netgen/include/order.hpp

@@ -2,7 +2,7 @@
 #define FILE_ORDER
 
 /* *************************************************************************/
-/* File:   order.hh                                                       */
+/* File:   order.hpp                                                      */
 /* Author: Joachim Schoeberl                                              */
 /* Date:   18. Jun. 97                                                    */
 /* *************************************************************************/
@@ -10,6 +10,8 @@
 
 namespace ngla
 {
+  using namespace ngcore;
+  using namespace ngstd;  
 
   /*
     reordering for sparse cholesky factoriztion

netgen/include/parallel_matrices.hpp

@@ -7,6 +7,11 @@
 /* Date:   2007,2011                                                      */
 /* ************************************************************************/
 
+
+#include <basematrix.hpp>
+#include <sparsematrix.hpp>
+#include <paralleldofs.hpp>
+
 namespace ngla
 {