ngsolve 6.2.2501.post21.dev1__cp312-cp312-macosx_10_15_universal2.whl → 6.2.2501.post37.dev1__cp312-cp312-macosx_10_15_universal2.whl

ngsolve/webgui.py

@@ -480,7 +480,7 @@ def BuildRenderData(mesh, func, order=2, draw_surf=True, draw_vol=True, intpoint
 
         pmin, pmax = [ngs.Vector(p) for p in zip(*pmima)]
         mesh_center = 0.5*(pmin+pmax)
-        mesh_radius = np.linalg.norm(pmax-pmin)/2
+        mesh_radius = float(np.linalg.norm(pmax-pmin)/2)
         
         pmat = pmat.reshape(-1, n_points_per_trig, 4)
 

{ngsolve-6.2.2501.post21.dev1.dist-info → ngsolve-6.2.2501.post37.dev1.dist-info}/METADATA

@@ -1,9 +1,9 @@
 Metadata-Version: 2.1
 Name: ngsolve
-Version: 6.2.2501.post21.dev1
+Version: 6.2.2501.post37.dev1
 Summary: NGSolve
 Author: The NGSolve team
 License: LGPL2.1
 License-File: LICENSE
-Requires-Dist: netgen-mesher ==6.2.2501.post24.dev0
+Requires-Dist: netgen-mesher ==6.2.2501.post37.dev0
 

ngsolve-6.2.2501.post37.dev1.dist-info/RECORD

@@ -0,0 +1,25 @@
+ngsolve/TensorProductTools.py,sha256=PKwn8c1LLeqFhPyhWtI01JZNztVz5_5sKQUYdnZgo5w,8353
+ngsolve/__console.py,sha256=8vkdKt1TaOxfSk-n8-4xeI4ATkaVIggYMEiZjxOtyWY,2889
+ngsolve/__expr.py,sha256=ZJFhvQ6S14AAwm3dVyEy-ixMwPw_ZOhMnOqyfXNVKnw,4215
+ngsolve/__init__.py,sha256=z2CG5eq5x97RUbozHnklsbG4TBder0zZVf6IpufgBzU,5432
+ngsolve/_scikit_build_core_dependencies.py,sha256=3a7XNqFjWNfA4ADwcNA7NhatXcHY2aHf1RuGtaFmmks,735
+ngsolve/bvp.py,sha256=oRcYhe17uOl3G2zTwf4CYsMDUTcFs0wb4AYu_cKADtk,1987
+ngsolve/directsolvers.py,sha256=ilZmA6_fQulvPZmstKfuyCFI8cj9neEleqKUHxqFpQw,995
+ngsolve/eigenvalues.py,sha256=sQ8cgtK2wRoSVU2OrR3NdJWf1OmLS-KbDnm5K9Av74s,9309
+ngsolve/internal.py,sha256=bXUxNpvi4zaskEt9QJ_IEiv_Hn43kg5JVv1wmMUYBfE,5113
+ngsolve/krylovspace.py,sha256=wBSsE4-lmBnJsiNU1BYBKjfYgFAcqR1Ioe84DMwIjIM,31863
+ngsolve/meshes.py,sha256=fxxtbxkfKuO2nUkXUiYwAwBAVxr9gkdnyIbZX3VpITQ,26786
+ngsolve/ngs2petsc.py,sha256=6978Z7Na6eMOZ5NTEZx0ZVflyX6Kyy4gPQT_9RYuNss,12499
+ngsolve/ngscxx.py,sha256=q3t5DgE3sJazWueas8FTXJ9jZM-LVa9KbGj_RTam6V4,1293
+ngsolve/nonlinearsolvers.py,sha256=8Q_CrD9vjdpZRV9sj2Fth7Y_kJJl7djPgarQmTK7Xfw,8216
+ngsolve/preconditioners.py,sha256=WSr3vnxskLHgnlwFVYpSyWJFpvrDUK6_30tATiBB5V8,208
+ngsolve/solvers.py,sha256=WgZ9X1x9mtQegXDu8VcTjP-wbS2EnrUS96Y8LXqeflw,304
+ngsolve/timestepping.py,sha256=SLbsgE5pY5tXYTLGlVypk4kS2Ui8GxOwWmKVxqneZRU,8350
+ngsolve/timing.py,sha256=Xqkq2byJjNNK9vJjHG2O3cMf6xuKIyySjk8yv6Le2T8,3985
+ngsolve/utils.py,sha256=V2yVF4K3uv1007-DEs_GrzJ3BOOh09c4lDGGR1Xw-PU,4557
+ngsolve/webgui.py,sha256=Oz6RN1_xQ8d6oMIx-KYIRDnBmlJIgwluCtZZjUOHtCs,24233
+ngsolve-6.2.2501.post37.dev1.dist-info/LICENSE,sha256=oZDcnIBDdV2Q-LCnX6ZrnkLUr0yYC_XdxjPwEk2zzuc,26430
+ngsolve-6.2.2501.post37.dev1.dist-info/METADATA,sha256=x-daFhH5e6QL2l7GV7mBGxJ6Ls6voQ2Xeuh-60YsFlo,200
+ngsolve-6.2.2501.post37.dev1.dist-info/WHEEL,sha256=Rs0tttgtJ0omR3rUP4GnW1J4ixuIyT0TIVbeEuGkiG4,110
+ngsolve-6.2.2501.post37.dev1.dist-info/top_level.txt,sha256=1llU44DxNIP85uPT4RChsIuV3h2d0_5j0Ui_Df5vV8E,15
+ngsolve-6.2.2501.post37.dev1.dist-info/RECORD,,

netgen/include/arnoldi.hpp

@@ -1,55 +0,0 @@
-#ifndef FILE_ARNOLDI
-#define FILE_ARNOLDI
-
-
-/**************************************************************************/
-/* File:   arnoldi.hpp                                                    */
-/* Author: Joachim Schoeberl                                              */
-/* Date:   5. Jul. 96                                                     */
-/**************************************************************************/
-
-
-#include "basematrix.hpp"
-
-namespace ngla
-{
-  /**
-     Arnoldi Eigenvalue Solver.
-
-     Solve the generalized evp
-     
-     A x = lam B x
-
-     B must by symmetric and (in theory) positive definite
-     A can be non-symmetric
-
-     It uses a shift-and-invert Arnoldi method 
-   */
-
-  template <typename SCAL>
-  class NGS_DLL_HEADER Arnoldi
-  {
-    shared_ptr<BaseMatrix> a;
-    shared_ptr<BaseMatrix> b;
-    shared_ptr<BitArray> freedofs;
-    SCAL shift;
-    optional<string> inversetype;
-  public:
-    Arnoldi (shared_ptr<BaseMatrix> aa, shared_ptr<BaseMatrix> ab, shared_ptr<BitArray> afreedofs = nullptr)
-      : a(aa), b(ab), freedofs(afreedofs)
-    { 
-      shift = 1.0;
-    }
-
-    void SetShift (SCAL ashift)
-    { shift = ashift; }
-    void SetInverseType (optional<string> ainv)
-    { inversetype = ainv; }
-
-    void Calc (int numval, Array<Complex> & lam, int nev, 
-               Array<shared_ptr<BaseVector>> & evecs, 
-               shared_ptr<BaseMatrix> pre = nullptr) const;
-  };
-}
-
-#endif

netgen/include/bandmatrix.hpp

@@ -1,334 +0,0 @@
-#ifndef FILE_BANDMATRIX
-#define FILE_BANDMATRIX
-
-/****************************************************************************/
-/* File:   bandmatrix.hpp                                                   */
-/* Author: Joachim Schoeberl                                                */
-/* Date:   14. Aug. 2002                                                    */
-/****************************************************************************/
-
-namespace ngbla
-{
-
-  /**
-     A symmetric band-matrix.
-  */
-  template <class T = double>
-  class FlatSymBandMatrix
-  {
-  protected:
-    /// the matrix size
-    int n;
-    /// number of bands in the triangular matrix
-    int bw;
-    /// the matrix entries
-    T *data;
-  public:
-    /// the according vector type
-    typedef typename mat_traits<T>::TV_COL TV;
-
-    /// Construction of FlatSymBandMatirx
-    FlatSymBandMatrix (int an, int abw, T * adata)
-      : n(an), bw(abw), data(adata)
-    { ; }
-
-    /// Matrix vector multiplication
-    void Mult (const FlatVector<TV> & x, FlatVector<TV> & y) const
-    {
-      for (int i = 0; i < n; i++)
-        y(i) = (*this)(i,i) * x(i);
-      for (int i = 0; i < n; i++)
-        for (int j = max2(i-bw+1, 0); j < i; j++)
-          {
-            y(i) += (*this)(i,j) * x(j);
-            y(j) += Trans((*this)(i,j)) * x(i);
-          }
-    }
-
-    /// Print matrix to stream
-    ostream & Print (ostream & ost) const
-    {
-      for (int i = 0; i < n; i++)
-        {
-          for (int j = 0; j < n; j++)
-            if (Used (i, j))
-              ost << setw(8) << (*this)(i,j) << " ";
-            else if (Used (j,i))
-              ost << setw(8) << "sym" << " ";
-            else
-              ost << setw(8) << 0;
-          ost << endl;
-        }
-      return ost;
-    }
-
-    /// matrix height
-    int Height() const { return n; }
-
-    /// matrix width
-    int BandWidth() const { return bw; }
-  
-    /// access operator. Assumes that $i \geq j$ referes to an element in the band
-    const T & operator() (int i, int j) const
-    { return data[i * bw + j - i + bw-1]; }
-
-    /// access operator. Assumes that $i \geq j$ referes to an element in the band
-    T & operator() (int i, int j) 
-    { return data[i * bw + j - i + bw-1]; }
-
-    /// check whether i, j refers to a valid element
-    bool Used (int i, int j) const
-    {
-      return (n > i && i >= j && j >= 0 && i-j < bw);
-    }
-  
-    /// assigns a constant value
-    FlatSymBandMatrix & operator= (const T & val)
-    {
-      for (int i = 0; i < bw * n; i++)
-        data[i] = val;
-      return *this;
-    }
-
-
-    /// computes required memory
-    static int RequiredMem (int n, int bw)
-    { return n*bw; }
-  };
-
-
-
-  /// output operator
-  template<typename T>
-  inline std::ostream & operator<< (std::ostream & s, const FlatSymBandMatrix<T> & m)
-  {
-    m.Print (s);
-    return s;
-  }
-
-
-
-
-  /**
-     A symmetric band-matrix with memory management.
-  */
-  template <class T = double>
-  class SymBandMatrix : public FlatSymBandMatrix<T>
-  {
-  public:
-    typedef typename mat_traits<T>::TV_COL TV;
-
-    /// Generates a symmetric band matrix
-    SymBandMatrix (int an, int abw)
-      : FlatSymBandMatrix<T> (an, abw, new T[an*abw])
-    { ; }
-
-    /// Deletes matrix
-    ~SymBandMatrix ()
-    { delete [] this->data; }
-
-    /// assigns a constant value
-    SymBandMatrix & operator= (const T & val)
-    {
-      for (int i = 0; i < this->bw * this->n; i++)
-        this->data[i] = val;
-      return *this;
-    }
-  };
-
-
-
-
-
-
-
-
-
-
-
-
-
-  /**
-     Cholesky factors of a band matrix.
-     This class does not provide memory management.
-
-
-     storage:
-
-     lfact (bw = 3)
-
-     \begin{verbatim}
-     d0                
-     0   d1            
-     1    2   d2        
-     3    4   d3   
-     \end{verbatim}
-  */
-  template <class T = double>
-  class FlatBandCholeskyFactors
-  {
-  protected:
-    /// matrix dimension
-    int n;
-    /// number of bands in the triangular matrix
-    int bw;
-    /// matrix matrix data, first diags, than lfact
-    T * mem; 
-  public:
-    // typedef typename mat_traits<T>::TV_COL TV;
-    typedef typename mat_traits<T>::TSCAL TSCAL;
-
-    /// assign dimension, bandwidth and memory
-    FlatBandCholeskyFactors (int an, int abw, T * amem)
-    { n = an, bw = abw, mem = amem; }
-
-    /// default constructor
-    FlatBandCholeskyFactors ()
-    { n = bw = 0; mem = 0; }
-
-    /// factor bandmatrix a
-    NGS_DLL_HEADER void Factor (const FlatSymBandMatrix<T> & a);
-
-    /// solve with factored matrices
-    template <class TVX, class TVY>
-    void Mult (const FlatVector<TVX> & x, FlatVector<TVY> & y) const
-    {
-      // const TVX * hx = x.Addr(0);
-      // TVY * hy = y.Addr(0);
-      FlatVector<TVX> hx = x;
-      FlatVector<TVY> hy = y;
-      const T * hm = &mem[0];
-
-      for (int i = 0; i < n; i++)
-        hy[i] = hx[i];
-
-      int i, jj = n;
-      for (i = 0; i < bw-1; i++)
-        {
-          typedef typename mat_traits<TVY>::TSCAL TTSCAL;
-          TVY sum = TTSCAL(0.0);  
-
-          for (int j = 0; j < i; j++, jj++)
-            sum += hm[jj] * hy[j];
-
-          hy[i] -= sum;
-        }
-
-      for (  ; i < n; i++)
-        {
-          typedef typename mat_traits<TVY>::TSCAL TTSCAL;
-          TVY sum = TTSCAL(0.0);
-
-          for (int j = i-bw+1; j < i; j++, jj++)
-            sum += hm[jj] * hy[j];
-
-          hy[i] -= sum;
-        }
-
-      for (int i = 0; i < n; i++)
-        {
-          TVY sum = mem[i] * hy[i];
-          hy[i] = sum;
-        }
-
-      // jj = n + (n-1) * (bw-1) - bw*(bw-1)/2;   
-      for (i = n-1; i >= bw-1; i--)
-        {
-          jj -= bw-1;
-          TVY val = hy[i];
-
-          int firstj = i-bw+1;
-          for (int j = 0; j < bw-1; j++)
-            hy[firstj+j] -= Trans (mem[jj+j]) * val;
-        }
-    
-      for (  ; i >= 0; i--) 
-        {
-          jj -= i;
-          TVY val = hy[i];
-
-          for (int j = 0; j < i; j++)
-            hy[j] -= Trans (mem[jj+j]) * val;
-        }
-    }
-
-
-
-    /// print matrix factors
-    ostream & Print (ostream & ost) const;
-
-    /// compute linear position of matrix element (i,j)
-    int Index (int i, int j) const
-    {
-      if (i < bw)
-        return n + (i * (i-1)) / 2 + j;
-      else
-        return n + i * (bw-2) + j - ((bw-1)*(bw-2))/2;   
-    }
-  
-    /// matrix element (i,j), (i,j) must be a valid position
-    const T & operator() (int i, int j) const
-    { 
-      if (i < bw)
-        return mem[n + (i * (i-1)) / 2 + j];
-      else
-        return mem[n + i * (bw-2) + j - ((bw-1)*(bw-2))/2];   
-    }
-
-    /// matrix element (i,j), (i,j) must be a valid position
-    T & operator() (int i, int j) 
-    {
-      if (i < bw)
-        return mem[n + (i * (i-1)) / 2 + j];
-      else
-        return mem[n + i * (bw-2) + j - ((bw-1)*(bw-2))/2];   
-    }
-
-    /// matrix size
-    int Size() const { return n; }
-    /// band-width of triangular matrix
-    int BandWidth() const { return bw; }
-    /// computes required memory
-    static int RequiredMem (int n, int bw)
-    { return n*bw - (bw * (bw-1)) / 2 + n; }
-  };
-  
-
-  ///  output operator.
-  template<typename T>
-  inline std::ostream & operator<< (std::ostream & s, const FlatBandCholeskyFactors<T> & m)
-  {
-    m.Print (s);
-    return s;
-  }
-
-
-
-
-
-  /**
-     Cholesky factors of a band matrix.
-  */
-  template <class T = double>
-  class BandCholeskyFactors : public FlatBandCholeskyFactors<T>
-  {
-  public:
-    /// allocate memory and factor the matrix a
-    BandCholeskyFactors (const SymBandMatrix<T> & a)
-      : FlatBandCholeskyFactors<T> (a.Height(), 
-                                    a.BandWidth(),
-                                    new T[FlatBandCholeskyFactors<T>::RequiredMem (a.Height(), a.BandWidth())])
-    {
-      this->Factor (a);
-    }
-
-    /// delete memory
-    ~BandCholeskyFactors ()
-    {
-      delete [] this->mem;
-    }
-  };
-
-}
-
-#endif