ngsolve 6.2.2506.post74.dev0__cp314-cp314-macosx_10_15_universal2.whl

netgen/include/ectypes.hpp

@@ -0,0 +1,121 @@
+#include "expr.hpp"
+
+// arithmetics with error control (better say error monitoring)
+
+
+namespace ngbla
+{
+  
+  
+  class EC_double
+  {
+    double val;
+    double err;
+  public:
+    EC_double () : val(0), err(0) { }
+    EC_double (double a)
+      : val(a), err(fabs(a)*1e-16) { }
+    EC_double (double a, double aerr)
+      : val(a), err(aerr)
+    {
+      Check();
+    }
+    EC_double (EC_double&) = default;    
+
+    EC_double & operator= (EC_double&) = default;
+    EC_double & operator= (double a)
+    {
+      val = a; err = fabs(a)*1e-16;
+      return *this;
+    }
+    
+    double Val() const { return val; }
+    double Err() const { return err; }    
+    operator double() const { return val; }
+
+    EC_double & operator*= (EC_double b)
+    {
+      *this = EC_double(*this)*b;
+      return *this;
+    }
+
+    EC_double & operator+= (EC_double b)
+    {
+      val += b.Val();
+      err += b.Err();
+      Check();
+      return *this;
+    }
+    
+    EC_double & operator-= (EC_double b)
+    {
+      *this = EC_double(*this) + (-b);
+      return *this;
+    }
+
+    EC_double operator- () const
+    {
+      return EC_double(-val, err);
+    }
+
+    void Check() const
+    {
+      if (err > 1e-4 * abs(val))
+        {
+          // throw Exception("roundoff is big, val = "+ToString(val)+", err = ", ToString(err));
+          cout << "roundoff is big, val = "+ToString(val)+", err = " + ToString(err) << endl;
+        }
+    }
+  };
+
+
+  inline EC_double operator* (EC_double a, EC_double b)
+  {
+    return EC_double(a.Val()*b.Val(), abs(a.Val())*b.Err()+abs(b.Val())*a.Err());
+  }
+
+  inline EC_double operator* (EC_double a, double b)
+  {
+    return a * EC_double(b);
+  }
+
+  inline EC_double operator* (double b, EC_double a)
+  {
+    return a * EC_double(b);
+  }
+  
+  inline EC_double operator+ (EC_double a, EC_double b)
+  {
+    return EC_double(a.Val()+b.Val(), a.Err()+b.Err());
+  }
+
+  typedef std::complex<EC_double> EC_Complex;
+
+  inline auto operator* (double a, EC_Complex b)
+  {
+    return EC_double(a) * b;
+  }
+
+  inline auto operator* (EC_Complex a, const Complex b)
+  {
+    return a*EC_Complex(b);
+  }
+
+  inline double L2Norm2 (EC_Complex val)
+  {
+    return sqr(val.real().Val()) + sqr(val.imag().Val());
+  }
+
+
+  // typedef std::complex<__128float> Complex128;
+  // typedef std::complex<double> Complex128;  
+
+
+
+  
+  template<> struct is_scalar_type<ngfem::EC_Complex> { static constexpr bool value = true; };
+  template<> struct is_scalar_type<ngfem::EC_double> { static constexpr bool value = true; };
+  template <> inline constexpr bool IsComplex<ngfem::EC_Complex> () { return true; }
+
+  
+}

netgen/include/eigen.hpp

@@ -0,0 +1,60 @@
+#ifndef FILE_EIGEN
+#define FILE_EIGEN
+
+/**************************************************************************/
+/* File:   eigen.hpp                                                      */
+/* Author: Joachim Schoeberl                                              */
+/* Date:   5. Jul. 96                                                     */
+/**************************************************************************/
+
+
+#include "basematrix.hpp"
+
+namespace ngla
+{
+
+  /**
+     Lanczos Eigen system calculation
+  */ 
+
+  class NGS_DLL_HEADER EigenSystem
+  {
+    ///
+    const BaseMatrix *a, *c;
+    ///
+    Array<double> ai, bi;
+    ///
+    double prec;
+    ///
+    int maxsteps;
+  
+  public:
+    ///
+    EigenSystem (const BaseMatrix & aa);
+    ///
+    EigenSystem (const BaseMatrix & aa, const BaseMatrix & ac);
+    ///
+    void SetMatrix (const BaseMatrix & aa);
+    ///
+    void SetPrecond (const BaseMatrix & ac);
+    ///
+    void SetMaxSteps (int amaxsteps);
+    ///
+    void SetPrecision (double aprec);
+
+    ///
+    int Calc();
+
+    ///
+    double EigenValue (int nr) const;
+    ///
+    double MaxEigenValue () const;
+    ///
+    int NumEigenValues () const;
+    ///
+    void PrintEigenValues (ostream & ost) const;
+  };
+
+}
+
+#endif

netgen/include/eigensystem.hpp

@@ -0,0 +1,18 @@
+#ifndef EIGENSYSTEM_HPP
+#define EIGENSYSTEM_HPP
+
+
+namespace ngbla
+{
+
+  /// Computes eigenvalues and vectors of the symmetric matrix mat.
+  extern NGS_DLL_HEADER void
+  CalcEigenSystem (const FlatMatrix<double> mat,
+                   FlatVector<double> lami,
+                   FlatMatrix<double> eigenvecs);
+
+
+}
+
+
+#endif