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/special_matrix.hpp

@@ -1,461 +0,0 @@
-#ifndef FILE_NGS_SPECIALMATRIX
-#define FILE_NGS_SPECIALMATRIX
-
-/* ************************************************************************/
-/* File:   special_matrix.hpp                                             */
-/* Author: Joachim Schoeberl                                              */
-/* Date:   14 Mar. 02                                                     */
-/* ************************************************************************/
-
-#include "basematrix.hpp"
-#include "vvector.hpp"
-
-namespace ngla
-{
-
-  // Convert RowMajor to ColMajor matrix (stored as vector)
-  class TransposeVector : public BaseMatrix
-  {
-    int h, w; // result matrix
-  public:
-    typedef double TSCAL;
-    
-    TransposeVector (int ah, int aw);
-    bool IsComplex() const override { return false; } 
-
-    int VHeight() const override { return h*w; }
-    int VWidth() const override { return h*w; }
-
-    ostream & Print (ostream & ost) const override;
-    
-    AutoVector CreateRowVector () const override;
-    AutoVector CreateColVector () const override;
-
-    void Mult (const BaseVector & x, BaseVector & y) const override;    
-    void MultTrans (const BaseVector & x, BaseVector & y) const override;    
-  };
-
-  
-  
-  class PermutationMatrix : public BaseMatrix
-  {
-    size_t width;
-    Array<size_t> ind;
-  public:
-    PermutationMatrix (size_t awidth, Array<size_t> aind)
-      : width(awidth), ind(aind) { ; } 
-
-    virtual bool IsComplex() const override { return false; } 
-
-    virtual int VHeight() const override { return ind.Size(); }
-    virtual int VWidth() const override { return width; }
-
-    virtual AutoVector CreateRowVector () const override
-    {
-      return CreateBaseVector(width, false, 1);
-    }
-    
-    virtual AutoVector CreateColVector () const override
-    {
-      return CreateBaseVector(ind.Size(), false, 1);
-    }
-
-    virtual void Mult (const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTrans (const BaseVector & x, BaseVector & y) const override;
-
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-  };
-
-
-  class Embedding : public BaseMatrix
-  {
-    size_t height;
-    IntRange range;
-    bool is_complex;
-  public:
-    Embedding (size_t aheight, IntRange arange, bool ais_complex = false)
-      : height(aheight), range(arange), is_complex(ais_complex) { ; }
-
-    virtual bool IsComplex() const override { return is_complex; } 
-
-    virtual int VHeight() const override { return height; }
-    virtual int VWidth() const override { return range.Size(); }
-
-    virtual AutoVector CreateRowVector () const override
-    {
-      return CreateBaseVector(range.Size(), is_complex, 1);
-    }
-    
-    virtual AutoVector CreateColVector () const override
-    {
-      return CreateBaseVector(height, is_complex, 1);
-    }
-
-    auto GetRange() const { return range; }
-    
-    virtual void Mult (const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTrans (const BaseVector & x, BaseVector & y) const override;
-
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual shared_ptr<BaseSparseMatrix> CreateSparseMatrix() const override;    
-  };
-
-
-  
-  class EmbeddedMatrix : public BaseMatrix
-  {
-    size_t height;
-    IntRange range;
-    shared_ptr<BaseMatrix> mat;
-  public:
-    EmbeddedMatrix (size_t aheight, IntRange arange, shared_ptr<BaseMatrix> amat)
-      : height(aheight), range(arange), mat(amat) { ; }
-
-    virtual bool IsComplex() const override { return mat->IsComplex(); } 
-
-    virtual int VHeight() const override { return height; }
-    virtual int VWidth() const override { return mat->VWidth(); }
-
-    virtual BaseMatrix::OperatorInfo GetOperatorInfo () const override;
-
-    virtual AutoVector CreateRowVector () const override
-    {
-      return mat->CreateRowVector();      
-    }
-    
-    virtual AutoVector CreateColVector () const override
-    {
-      return CreateBaseVector(height, IsComplex(), 1);      
-    }
-
-    virtual void Mult (const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTrans (const BaseVector & x, BaseVector & y) const override;
-
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    shared_ptr<BaseMatrix> GetMatrix() const { return mat; }
-    IntRange GetRange() const { return range; }
-
-    virtual shared_ptr<BaseSparseMatrix> CreateSparseMatrix() const override;
-    
-    
-    /*
-    virtual shared_ptr<BaseMatrix> CreateDeviceMatrix() const override
-    {
-      return make_shared<EmbeddedMatrix>(height, range, mat->CreateDeviceMatrix());
-    }
-    */
-  };
-
-
-  class EmbeddingTranspose : public BaseMatrix
-  {
-    size_t width;
-    IntRange range;
-    bool is_complex;
-  public:
-    EmbeddingTranspose (size_t awidth, IntRange arange, bool ais_complex = false)
-      : width(awidth), range(arange), is_complex(ais_complex) { ; }
-
-    virtual bool IsComplex() const override { return is_complex; } 
-    
-    virtual int VHeight() const override { return range.Size(); }
-    virtual int VWidth() const override { return width; }
-
-    virtual AutoVector CreateRowVector () const override
-    {
-      return CreateBaseVector(width, is_complex, 1);
-    }
-    
-    virtual AutoVector CreateColVector () const override
-    {
-      return CreateBaseVector(range.Size(), is_complex, 1);
-    }
-
-    auto GetRange() const { return range; }
-    
-    virtual void Mult (const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTrans (const BaseVector & x, BaseVector & y) const override;
-
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-  };
-
-
-  class EmbeddedTransposeMatrix : public BaseMatrix
-  {
-    size_t width;
-    IntRange range;
-    shared_ptr<BaseMatrix> mat;
-  public:
-    EmbeddedTransposeMatrix (size_t awidth, IntRange arange, shared_ptr<BaseMatrix> amat)
-      : width(awidth), range(arange), mat(amat) { ; }
-
-    virtual bool IsComplex() const override { return mat->IsComplex(); } 
-
-    virtual BaseMatrix::OperatorInfo GetOperatorInfo () const override;
-    
-    virtual int VHeight() const override { return mat->Height(); }
-    virtual int VWidth() const override { return width; }
-
-    virtual AutoVector CreateRowVector () const override
-    {
-      return CreateBaseVector(width, mat->IsComplex(), 1);
-    }
-    
-    virtual AutoVector CreateColVector () const override
-    {
-      return mat->CreateColVector();
-    }
-
-    virtual void Mult (const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTrans (const BaseVector & x, BaseVector & y) const override;
-
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    shared_ptr<BaseMatrix> GetMatrix() const { return mat; }
-    IntRange GetRange() const { return range; }
-
-    /*
-    virtual shared_ptr<BaseMatrix> CreateDeviceMatrix() const override
-    {
-      return make_shared<EmbeddedTransposeMatrix>(width, range, mat->CreateDeviceMatrix());
-    }
-    */
-  };
-
-
-
-  
-  
-  template <class TVR, class TVC>
-  class Real2ComplexMatrix : public BaseMatrix
-  {
-    shared_ptr<BaseMatrix> realmatrix;
-    VVector<TVR> hx, hy;
-  public:
-    NGS_DLL_HEADER Real2ComplexMatrix (shared_ptr<BaseMatrix> arealmatrix = nullptr);
-    bool IsComplex() const override { return true; }     
-    void SetMatrix (shared_ptr<BaseMatrix> arealmatrix);
-    const BaseMatrix & GetMatrix () const { return *realmatrix; }
-    void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    void MultAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-
-    AutoVector CreateRowVector() const override;
-    AutoVector CreateColVector() const override;
-  };
-
-
-
-
-  ////////////////////////////////////////////////////////////////////////////////
-  // added 08/19/2003
-  template <class TVR>
-  class NGS_DLL_HEADER Sym2NonSymMatrix : public BaseMatrix
-  {
-    const BaseMatrix * base;
-    VVector<TVR> hx, hy;
-  public:
-    Sym2NonSymMatrix (const BaseMatrix * abasematrix = 0);
-    void SetMatrix (const BaseMatrix * abasematrix);
-    const BaseMatrix & GetMatrix () const { return *base; }
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const;
-    //  virtual void MultAdd (Complex s, const BaseVector & x, BaseVector & y) const;
-  };
-
-
-
-  ////////////////////////////////////////////////////////////////////////////////
-  // added 09/02/2003
-  template <class TVSMALL, class TVBIG>
-  class NGS_DLL_HEADER Small2BigNonSymMatrix : public BaseMatrix
-  {
-    const BaseMatrix * base;
-    VVector<TVSMALL> hx1, hx2, hy1, hy2;
-  public:
-    Small2BigNonSymMatrix (const BaseMatrix * abasematrix = 0);
-    void SetMatrix (const BaseMatrix * abasematrix);
-    bool IsComplex() const override { return base->IsComplex(); }
-    const BaseMatrix & GetMatrix () const { return *base; }
-    void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    //  virtual void MultAdd (Complex s, const BaseVector & x, BaseVector & y) const;
-    void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    AutoVector CreateRowVector() const override
-    { throw Exception("CreateRowVector not implemented for Small2BigNonSymMatrix!"); }
-    AutoVector CreateColVector() const override
-    { throw Exception("CreateColVector not implemented for Small2BigNonSymMatrix!"); }
-  };
-
-
-  class BlockMatrix : public BaseMatrix
-  {
-    Array<Array<shared_ptr<BaseMatrix>>> mats;
-    size_t h, w;
-
-    // one matrix per row/col that can be used to create vectors etc.
-    Array<shared_ptr<BaseMatrix>> row_reps;
-    Array<shared_ptr<BaseMatrix>> col_reps;
-
-  public:
-    BlockMatrix (const Array<Array<shared_ptr<BaseMatrix>>> & amats);
-
-    bool IsComplex() const override { return row_reps[0]->IsComplex(); }
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    const shared_ptr<BaseMatrix> & operator()(size_t i, size_t j)
-    {
-      if (i >= h) throw Exception("Tried to access BlockMatrix row that is out of range");
-      if (j >= w) throw Exception("Tried to access BlockMatrix col that is out of range");
-      return mats[i][j];
-    }
-
-    size_t BlockRows() const { return h; }
-    size_t BlockCols() const { return w; }
-
-    virtual int VHeight() const override { throw Exception("VHeight does not make sense for BlockMatrix");}
-    virtual int VWidth() const override { throw Exception("VWidth does not make sense for BlockMatrix");}
-
-    virtual AutoVector CreateRowVector () const override;
-    virtual AutoVector CreateColVector () const override;
-  };
-
-
-  
-  class BaseMatrixFromVector : public BaseMatrix
-  {
-    shared_ptr<BaseVector> vec;
-
-  public:
-    BaseMatrixFromVector (shared_ptr<BaseVector> avec);
-
-    bool IsComplex() const override { return vec->IsComplex(); }
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    virtual int VHeight() const override { return vec->Size(); }
-    virtual int VWidth() const override { return 1; }
-
-    virtual AutoVector CreateRowVector () const override;
-    virtual AutoVector CreateColVector () const override;
-  };
-
-
-  class BaseMatrixFromMultiVector : public BaseMatrix
-  {
-    shared_ptr<MultiVector> vec;
-
-  public:
-    BaseMatrixFromMultiVector (shared_ptr<MultiVector> avec);
-
-    bool IsComplex() const override { return vec->IsComplex(); }
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-
-    virtual int VHeight() const override { return vec->RefVec()->Size(); }
-    virtual int VWidth() const override { return vec->Size(); }
-
-    virtual AutoVector CreateRowVector () const override;
-    virtual AutoVector CreateColVector () const override;
-  };
-
-
-  template <typename T = double>
-  class BaseMatrixFromMatrix : public BaseMatrix
-  {
-    Matrix<T> mat;
-
-  public:
-    NGS_DLL_HEADER BaseMatrixFromMatrix (Matrix<T> amat);
-
-    bool IsComplex() const override { return typeid(T) == typeid(Complex); }
-    virtual void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-    virtual void MultTransAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-
-    virtual int VHeight() const override { return mat.Height(); }
-    virtual int VWidth() const override { return mat.Width(); }
-    virtual size_t NZE() const override { return mat.Height()*mat.Width(); }
-    virtual AutoVector CreateRowVector () const override;
-    virtual AutoVector CreateColVector () const override;
-  };
-
-
-
-
-  class LoggingMatrix : public BaseMatrix
-  {
-    shared_ptr<BaseMatrix> mat;
-    string label;
-    unique_ptr<ostream> out;
-    optional<NgMPI_Comm> comm;
-  public:
-    LoggingMatrix (shared_ptr<BaseMatrix> amat, string alabel, string filename,
-                   optional<NgMPI_Comm> acomm);
-    ~LoggingMatrix ();
-    int VHeight() const override { return mat->VHeight(); }
-    int VWidth() const override { return mat->VWidth(); }
-    bool IsComplex() const override { return mat->IsComplex(); }
-    
-    BaseVector & AsVector() override;
-    const BaseVector & AsVector() const override;
-    void SetZero() override;
-
-    ostream & Print (ostream & ost) const override { return mat->Print(ost); }
-    Array<MemoryUsage> GetMemoryUsage () const override { return mat->GetMemoryUsage(); }
-    size_t NZE () const override { return mat->NZE(); }
-
-    void Update() override { mat->Update(); }
-    shared_ptr<BaseMatrix> CreateMatrix () const override { return mat->CreateMatrix(); }
-    AutoVector CreateRowVector () const override;
-    AutoVector CreateColVector () const override;
-
-    void Mult (const BaseVector & x, BaseVector & y) const override;
-    void MultTrans (const BaseVector & x, BaseVector & y) const override;
-    void MultAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    void MultAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-    void MultTransAdd (double s, const BaseVector & x, BaseVector & y) const override;
-    void MultTransAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-    void MultConjTransAdd (Complex s, const BaseVector & x, BaseVector & y) const override;
-    void MultAdd (FlatVector<double> alpha, const MultiVector & x, MultiVector & y) const override;
-    
-    void MultAdd1 (double s, const BaseVector & x, BaseVector & y,
-                   const BitArray * ainner = NULL,
-                   const Array<int> * acluster = NULL) const override
-    { mat->MultAdd1 (s, x, y, ainner, acluster); }
-    
-    void MultAdd2 (double s, const BaseVector & x, BaseVector & y,
-                   const BitArray * ainner = NULL,
-                   const Array<int> * acluster = NULL) const override
-    { mat->MultAdd2 (s, x, y, ainner, acluster); }
-
-    shared_ptr<BaseMatrix> InverseMatrix (shared_ptr<BitArray> subset = nullptr) const override
-    { return mat->InverseMatrix(subset); }
-    
-    shared_ptr<BaseMatrix> InverseMatrix (shared_ptr<const Array<int>> clusters) const override
-    { return mat->InverseMatrix(clusters); }
-    
-    INVERSETYPE SetInverseType ( INVERSETYPE ainversetype ) const override
-    { return mat->SetInverseType(ainversetype); }
-    INVERSETYPE SetInverseType ( string ainversetype ) const override
-    { return mat->SetInverseType(ainversetype); }
-    INVERSETYPE  GetInverseType () const override
-    { return mat->GetInverseType(); }
-    void DoArchive (Archive & ar) override
-    { mat->DoArchive(ar); }
-  };
-  
-
-
-
-  
-}
-
-
-#endif

netgen/include/specialelement.hpp

@@ -1,125 +0,0 @@
-#ifndef FILE_SPECIALELEMENT
-#define FILE_SPECIALELEMENT
-
-/*********************************************************************/
-/* File:   specialelement.hpp                                        */
-/* Author: Joachim Schoeberl                                         */
-/* Date:   28. Mar. 2002                                             */
-/*********************************************************************/
-
-#include "integrator.hpp"
-
-
-namespace ngfem
-{
-
-  /*
-    Something special ...
-
-    Finite-Element + Integrator
-
-    E.g., a Contact-Element  
-  */
-
-
-
-  class NGS_DLL_HEADER SpecialElement
-  {
-  public:
-    SpecialElement ()
-    {
-      ;
-    }
-    virtual ~SpecialElement() { ; }
-  
-
-
-    virtual void GetDofNrs (Array<int> & dnums) const = 0;
-    
-    // dofs for test-space (which might be different from trial-space)
-    virtual void GetDofNrs2 (Array<int> & dnums) const
-    {
-      GetDofNrs(dnums);
-    }
-    
-    virtual double Energy (FlatVector<double> elx,
-			   LocalHeap & lh) const
-    {
-      return 0;
-    }
-    virtual double Energy (FlatVector<Complex> elx,
-			   LocalHeap & lh) const 
-    {
-      cerr << "SpecialElement::Energy (complex) called" << endl;
-      return 0;
-    }
-
-    template<int S, class T>
-    void Apply (FlatVector< Vec<S,T> > elx,
-		FlatVector< Vec<S,T> > ely,
-		LocalHeap & lh) const
-    {
-      cerr << "SpecialElement::Apply (Vec) called" << endl;
-    }
-
-    virtual void Apply (FlatVector<double> elx, FlatVector<double> ely, 
-			LocalHeap & lh) const;
-
-    virtual void Apply (FlatVector<Complex> elx,
-			FlatVector<Complex> ely,
-			LocalHeap & lh) const
-    {
-      cerr << "SpecialElement::Apply (complex) called" << endl;
-    }
-
-
-    virtual void CalcElementMatrix(FlatMatrix<double> elmat,
-			   LocalHeap& lh) const;
-
-    virtual void CalcElementMatrix(FlatMatrix<Complex> elmat,
-                                   LocalHeap & lh) const;
-    /*
-    {
-      cerr << "SpecialElement::CalcElementMatrix(complex) called" << endl;
-      exit(10);
-      FlatMatrix<double> relmat;
-      CalcElementMatrix(relmat, lh);
-      elmat.AssignMemory (relmat.Height(), relmat.Width(), lh);
-      elmat = relmat;
-    }
-    */
-    
-    virtual void CalcElementVector(FlatVector<double> elvec,
-			   LocalHeap & lh) const;
-
-    virtual void CalcElementVector(FlatVector<Complex> elvec,
-                                   LocalHeap & lh) const;
-    /*
-    {
-      cerr << "SpecialElement::CalcElementMatrix(complex) called" << endl;
-      exit(10);
-      FlatVector<double> relvec;
-      CalcElementVector(relvec, lh);
-      elvec.AssignMemory (relvec.Height(),lh);
-      elvec = relvec;
-    }
-    */
-    
-    virtual void CalcLinearizedElementMatrix(FlatVector<double> elx,
-                                             FlatMatrix<double> elmat,
-                                             LocalHeap& lh) const
-    {
-      CalcElementMatrix(elmat, lh);
-    }
-    virtual void CalcLinearizedElementMatrix(FlatVector<Complex> elx,
-                                             FlatMatrix<Complex> elmat,
-                                             LocalHeap& lh) const
-    {
-      throw Exception("Complex CalcLinearizedElementMatrix not implemented!");
-    }
-
-  };
-
-}
-
-#endif

netgen/include/statushandler.hpp

@@ -1,33 +0,0 @@
-#ifndef _STATUSHANDLER_HPP
-#define _STATUSHANDLER_HPP
-
-namespace ngstd
-{
-
-  /** Access to statusbar. (and more)
-   */
-
-  class NGS_DLL_HEADER BaseStatusHandler
-  {
-  public:
-    static void PushStatus (const char * str);
-    static void PopStatus ();
-    static void SetThreadPercentage (double percent);
-
-    static void GetStatus (string & str, double & percent);
-
-    static void SetTerminate(void);
-    static void UnSetTerminate(void);
-    static bool ShouldTerminate(void);
-
-    class Region
-    {
-    public:
-      Region(const char * str) { PushStatus(str); }
-      ~Region() { PopStatus(); }
-    };
-  };
-
-}
-
-#endif // _STATUSHANDLER_HPP

netgen/include/stringops.hpp

@@ -1,12 +0,0 @@
-#ifndef _STRINGOPS_HPP
-#define _STRINGOPS_HPP
-
-namespace ngstd
-{
-
-  // NGS_DLL_HEADER bool StringFitsPattern(const string & str, const string & pattern);
-  NGS_DLL_HEADER bool StringFitsPattern(string_view str, string_view pattern);
-
-}
-
-#endif // _STRINGOPS_HPP