ddfem 1.0.10__py3-none-any.whl → 1.0.11__py3-none-any.whl

ddfem/data/2.11.dev20250709/{localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44.cc → localfunction_4adfaf8a08fac23151af0b1e790db1f7_a524c1196983e65de1c06d7d6afdeb44.cc}

@@ -1,5 +1,5 @@
-#ifndef GUARD_a8c726f728cd35d137188b33301aeef1
-#define GUARD_a8c726f728cd35d137188b33301aeef1
+#ifndef GUARD_4adfaf8a08fac23151af0b1e790db1f7
+#define GUARD_4adfaf8a08fac23151af0b1e790db1f7
 
 #define USING_DUNE_PYTHON 1
 #include <config.h>
@@ -20,13 +20,13 @@
 #include <dune/fempy/function/virtualizedgridfunction.hh>
 #include <dune/fempy/py/ufllocalfunction.hh>
 
-namespace UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1
+namespace UFLLocalFunctions_4adfaf8a08fac23151af0b1e790db1f7
 {
 
   // UFLLocalFunction
 // ----------------
 
-template< class GridPart, class Coeffddm >
+template< class GridPart, class Coeffuh >
 struct UFLLocalFunction
   : public Dune::Fem::BindableGridFunctionWithSpace<GridPart,Dune::Dim<1>>
 {
@@ -42,14 +42,14 @@ struct UFLLocalFunction
   typedef Dune::Fem::IntersectionSide Side;
   typedef std::tuple<> ConstantTupleType;
   typedef std::tuple< Dune::Fem::GridFunctionSpace< GridPartType, Dune::FieldVector< double, 1 > > > CoefficientFunctionSpaceTupleType;
-  typedef std::tuple< Coeffddm > CoefficientTupleType;
+  typedef std::tuple< Coeffuh > CoefficientTupleType;
   template< std::size_t i >
   using CoefficientFunctionSpaceType = std::tuple_element_t< i, CoefficientFunctionSpaceTupleType >;
   template< std::size_t i >
   using CoefficientRangeType = typename CoefficientFunctionSpaceType< i >::RangeType;
   template< std::size_t i >
   using CoefficientJacobianRangeType = typename CoefficientFunctionSpaceType< i >::JacobianRangeType;
-  static constexpr bool gridPartValid = Dune::Fem::checkGridPartValid<GridPartType,Dune::Fem::ConstLocalFunction<Coeffddm>>();
+  static constexpr bool gridPartValid = Dune::Fem::checkGridPartValid<GridPartType,Dune::Fem::ConstLocalFunction<Coeffuh>>();
   template< std::size_t i >
   using CoefficientType = std::tuple_element_t< i, CoefficientTupleType >;
   template< std::size_t i >
@@ -57,9 +57,9 @@ struct UFLLocalFunction
   using BaseType::entity;
   using BaseType::geometry;
 
-  UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Coeffddm &coeffddm, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
+  UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Coeffuh &coeffuh, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
     : BaseType(gridPart,name,order),
-      coefficients_( Dune::Fem::ConstLocalFunction< Coeffddm >( coeffddm ) )
+      coefficients_( Dune::Fem::ConstLocalFunction< Coeffuh >( coeffuh ) )
   {}
 
   void bind ( const IntersectionType &intersection, Side side )
@@ -156,19 +156,19 @@ struct UFLLocalFunction
     return result;
   }
   ConstantTupleType constants_;
-  std::tuple< Dune::Fem::ConstLocalFunction< Coeffddm > > coefficients_;
+  std::tuple< Dune::Fem::ConstLocalFunction< Coeffuh > > coefficients_;
 };
 
-} // namespace UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1
+} // namespace UFLLocalFunctions_4adfaf8a08fac23151af0b1e790db1f7
 
-PYBIND11_MODULE( localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44, module )
+PYBIND11_MODULE( localfunction_4adfaf8a08fac23151af0b1e790db1f7_a524c1196983e65de1c06d7d6afdeb44, module )
 {
-  typedef UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > > LocalFunctionType;
+  typedef UFLLocalFunctions_4adfaf8a08fac23151af0b1e790db1f7::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > > LocalFunctionType;
   if constexpr( LocalFunctionType::gridPartValid )
   {
-      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_a8c726f728cd35d137188b33301aeef1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_a8c726f728cd35d137188b33301aeef1_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_4adfaf8a08fac23151af0b1e790db1f7::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > >, Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_4adfaf8a08fac23151af0b1e790db1f7_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
       Dune::FemPy::registerUFLLocalFunction( module, cls );
-      cls.def( pybind11::init( [] ( pybind11::object gridView, const std::string &name, int order, const Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > &coeffddm ) {return new LocalFunctionType( Dune::FemPy::gridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>(gridView),name,order, coeffddm); } ), pybind11::keep_alive< 1, 2 >(), pybind11::keep_alive< 1, 3 >() );
+      cls.def( pybind11::init( [] ( pybind11::object gridView, const std::string &name, int order, const Dune::FemPy::VirtualizedGridFunction< Dune::FemPy::GridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>, Dune::FieldVector< double, 1 > > &coeffuh ) {return new LocalFunctionType( Dune::FemPy::gridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>(gridView),name,order, coeffuh); } ), pybind11::keep_alive< 1, 2 >(), pybind11::keep_alive< 1, 3 >() );
       cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
   }
 }

ddfem/data/2.11.dev20250709/{localfunction_00f2e5593eeb23f01906255d67244f6e_a524c1196983e65de1c06d7d6afdeb44.cc → localfunction_4f6a526abfe63959f81b09a5796c5eb1_a524c1196983e65de1c06d7d6afdeb44.cc}

@@ -1,5 +1,5 @@
-#ifndef GUARD_00f2e5593eeb23f01906255d67244f6e
-#define GUARD_00f2e5593eeb23f01906255d67244f6e
+#ifndef GUARD_4f6a526abfe63959f81b09a5796c5eb1
+#define GUARD_4f6a526abfe63959f81b09a5796c5eb1
 
 #define USING_DUNE_PYTHON 1
 #include <config.h>
@@ -19,7 +19,7 @@
 #include <dune/common/exceptions.hh>
 #include <dune/fempy/py/ufllocalfunction.hh>
 
-namespace UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e
+namespace UFLLocalFunctions_4f6a526abfe63959f81b09a5796c5eb1
 {
 
   // UFLLocalFunction
@@ -39,10 +39,7 @@ struct UFLLocalFunction
   typedef typename EntityType::Geometry Geometry;
   typedef typename Geometry::GlobalCoordinate GlobalCoordinateType;
   typedef Dune::Fem::IntersectionSide Side;
-  typedef double Conc2;
-  typedef std::tuple< std::shared_ptr< Conc2 > > ConstantTupleType;
-  template< std::size_t i >
-  using ConstantsRangeType = typename std::tuple_element_t< i, ConstantTupleType >::element_type;
+  typedef std::tuple<> ConstantTupleType;
   typedef std::tuple<> CoefficientTupleType;
   static constexpr bool gridPartValid = true;
   template< std::size_t i >
@@ -54,9 +51,7 @@ struct UFLLocalFunction
 
   UFLLocalFunction ( const GridPartType &gridPart, const std::string &name, int order, const Dune::Fem::ParameterReader &parameter = Dune::Fem::Parameter::container() )
     : BaseType(gridPart,name,order)
-  {
-    std::get< 0 >( constants_ ) = std::make_shared< Conc2 >( (Conc2(0)) );
-  }
+  {}
 
   void bind ( const IntersectionType &intersection, Side side )
   {
@@ -76,20 +71,51 @@ struct UFLLocalFunction
   template< class Point >
   void evaluate ( const Point &x, typename FunctionSpaceType::RangeType &result ) const
   {
-    double tmp0 = constant< 0 >();
-    result[ 0 ] = tmp0;
+    using std::max;
+    using std::min;
+    using std::sqrt;
+    GlobalCoordinateType tmp0 = geometry().global( Dune::Fem::coordinate( x ) );
+    const auto tmp1 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp2 = tmp0[ 0 ] * tmp0[ 0 ];
+    const auto tmp3 = tmp2 + tmp1;
+    const auto tmp4 = 1e-10 + tmp3;
+    const auto tmp5 = std::sqrt( tmp4 );
+    const auto tmp6 = -0.4 + tmp5;
+    const auto tmp7 = 0.15 + tmp0[ 1 ];
+    const auto tmp8 = tmp7 * tmp7;
+    const auto tmp9 = 0.15 + tmp0[ 0 ];
+    const auto tmp10 = tmp9 * tmp9;
+    const auto tmp11 = tmp10 + tmp8;
+    const auto tmp12 = 1e-10 + tmp11;
+    const auto tmp13 = std::sqrt( tmp12 );
+    const auto tmp14 = -0.3 + tmp13;
+    const auto tmp15 = -0.15 + tmp0[ 1 ];
+    const auto tmp16 = tmp15 * tmp15;
+    const auto tmp17 = -0.15 + tmp0[ 0 ];
+    const auto tmp18 = tmp17 * tmp17;
+    const auto tmp19 = tmp18 + tmp16;
+    const auto tmp20 = 1e-10 + tmp19;
+    const auto tmp21 = std::sqrt( tmp20 );
+    const auto tmp22 = -0.3 + tmp21;
+    const auto tmp23 = std::min( tmp22, tmp14 );
+    const auto tmp24 = std::max( tmp23, tmp6 );
+    result[ 0 ] = (tmp24 <= 0.0 ? 1 : 0.0);
   }
 
   template< class Point >
   void jacobian ( const Point &x, typename FunctionSpaceType::JacobianRangeType &result ) const
   {
-    result=typename FunctionSpaceType::JacobianRangeType(0);
+    (result[ 0 ])[ 0 ] = 0.0;
+    (result[ 0 ])[ 1 ] = 0.0;
   }
 
   template< class Point >
   void hessian ( const Point &x, typename FunctionSpaceType::HessianRangeType &result ) const
   {
-    result=typename FunctionSpaceType::HessianRangeType(0);
+    ((result[ 0 ])[ 0 ])[ 0 ] = 0.0;
+    ((result[ 0 ])[ 0 ])[ 1 ] = 0.0;
+    ((result[ 0 ])[ 1 ])[ 0 ] = 0.0;
+    ((result[ 0 ])[ 1 ])[ 1 ] = 0.0;
   }
 
   template< std::size_t i >
@@ -103,31 +129,20 @@ struct UFLLocalFunction
   {
     return *std::get< i >( constants_ );
   }
-
-  const Conc2 &conc2 () const
-  {
-    return *std::get< 0 >( constants_ );
-  }
-
-  Conc2 &conc2 ()
-  {
-    return *std::get< 0 >( constants_ );
-  }
   ConstantTupleType constants_;
   std::tuple<  > coefficients_;
 };
 
-} // namespace UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e
+} // namespace UFLLocalFunctions_4f6a526abfe63959f81b09a5796c5eb1
 
-PYBIND11_MODULE( localfunction_00f2e5593eeb23f01906255d67244f6e_a524c1196983e65de1c06d7d6afdeb44, module )
+PYBIND11_MODULE( localfunction_4f6a526abfe63959f81b09a5796c5eb1_a524c1196983e65de1c06d7d6afdeb44, module )
 {
-  typedef UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > > LocalFunctionType;
+  typedef UFLLocalFunctions_4f6a526abfe63959f81b09a5796c5eb1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > > LocalFunctionType;
   if constexpr( LocalFunctionType::gridPartValid )
   {
-      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_00f2e5593eeb23f01906255d67244f6e_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_4f6a526abfe63959f81b09a5796c5eb1::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > >"), Dune::Python::IncludeFiles({"python/dune/generated/localfunction_4f6a526abfe63959f81b09a5796c5eb1_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
       Dune::FemPy::registerUFLLocalFunction( module, cls );
       cls.def( pybind11::init( [] ( pybind11::object gridView, const std::string &name, int order ) {return new LocalFunctionType( Dune::FemPy::gridPart<Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true >>(gridView),name,order); } ), pybind11::keep_alive< 1, 2 >() );
-      cls.def_property( "c2", [] ( LocalFunctionType &self ) -> UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > >::Conc2 { return self.conc2(); }, [] ( LocalFunctionType &self, const UFLLocalFunctions_00f2e5593eeb23f01906255d67244f6e::UFLLocalFunction< typename Dune::FemPy::GridPart< Dune::Fem::FilteredGridPart< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > >, Dune::Fem::SimpleFilter< Dune::FemPy::GridPart< Dune::Fem::AdaptiveLeafGridPart< Dune::ALUGrid< 2, 2, Dune::simplex > > > >, true > > >::Conc2 &v ) { self.conc2() = v; } );
       cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
   }
 }