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

ddfem/data/2.11.dev20250709/{localfunction_d75993d7ae5919d23117b153e900840d_a524c1196983e65de1c06d7d6afdeb44.cc → localfunction_3d957c5b2628e5120636b8e3f461bc0d_a524c1196983e65de1c06d7d6afdeb44.cc}

@@ -1,5 +1,5 @@
-#ifndef GUARD_d75993d7ae5919d23117b153e900840d
-#define GUARD_d75993d7ae5919d23117b153e900840d
+#ifndef GUARD_3d957c5b2628e5120636b8e3f461bc0d
+#define GUARD_3d957c5b2628e5120636b8e3f461bc0d
 
 #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_d75993d7ae5919d23117b153e900840d
+namespace UFLLocalFunctions_3d957c5b2628e5120636b8e3f461bc0d
 {
 
   // UFLLocalFunction
@@ -39,10 +39,7 @@ struct UFLLocalFunction
   typedef typename EntityType::Geometry Geometry;
   typedef typename Geometry::GlobalCoordinate GlobalCoordinateType;
   typedef Dune::Fem::IntersectionSide Side;
-  typedef double Conc3;
-  typedef std::tuple< std::shared_ptr< Conc3 > > 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< Conc3 >( (Conc3(0)) );
-  }
+  {}
 
   void bind ( const IntersectionType &intersection, Side side )
   {
@@ -76,20 +71,71 @@ 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 = -0.4 + tmp0[ 0 ];
+    const auto tmp2 = tmp0[ 1 ] * tmp0[ 1 ];
+    const auto tmp3 = tmp0[ 0 ] * tmp0[ 0 ];
+    const auto tmp4 = tmp3 + tmp2;
+    const auto tmp5 = 1e-10 + tmp4;
+    const auto tmp6 = std::sqrt( tmp5 );
+    const auto tmp7 = -1.4 + tmp6;
+    const auto tmp8 = std::min( tmp7, tmp1 );
+    const auto tmp9 = -1 * tmp0[ 0 ];
+    const auto tmp10 = -0.4 + tmp9;
+    const auto tmp11 = std::min( tmp7, tmp10 );
+    const auto tmp12 = std::max( tmp11, tmp8 );
+    const auto tmp13 = 1 + tmp0[ 0 ];
+    const auto tmp14 = tmp13 * tmp13;
+    const auto tmp15 = tmp2 + tmp14;
+    const auto tmp16 = 1e-10 + tmp15;
+    const auto tmp17 = std::sqrt( tmp16 );
+    const auto tmp18 = -0.5 + tmp17;
+    const auto tmp19 = -1 + tmp0[ 0 ];
+    const auto tmp20 = tmp19 * tmp19;
+    const auto tmp21 = tmp2 + tmp20;
+    const auto tmp22 = 1e-10 + tmp21;
+    const auto tmp23 = std::sqrt( tmp22 );
+    const auto tmp24 = -0.5 + tmp23;
+    const auto tmp25 = 0.8 + tmp0[ 1 ];
+    const auto tmp26 = tmp25 * tmp25;
+    const auto tmp27 = tmp3 + tmp26;
+    const auto tmp28 = 1e-10 + tmp27;
+    const auto tmp29 = std::sqrt( tmp28 );
+    const auto tmp30 = -0.5 + tmp29;
+    const auto tmp31 = -1 * tmp30;
+    const auto tmp32 = -0.8 + tmp0[ 1 ];
+    const auto tmp33 = tmp32 * tmp32;
+    const auto tmp34 = tmp3 + tmp33;
+    const auto tmp35 = 1e-10 + tmp34;
+    const auto tmp36 = std::sqrt( tmp35 );
+    const auto tmp37 = -0.5 + tmp36;
+    const auto tmp38 = -1 * tmp37;
+    const auto tmp39 = -1 + tmp6;
+    const auto tmp40 = std::max( tmp39, tmp38 );
+    const auto tmp41 = std::max( tmp40, tmp31 );
+    const auto tmp42 = std::min( tmp41, tmp24 );
+    const auto tmp43 = std::min( tmp42, tmp18 );
+    const auto tmp44 = std::max( tmp43, tmp12 );
+    result[ 0 ] = (tmp44 <= 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 +149,20 @@ struct UFLLocalFunction
   {
     return *std::get< i >( constants_ );
   }
-
-  const Conc3 &conc3 () const
-  {
-    return *std::get< 0 >( constants_ );
-  }
-
-  Conc3 &conc3 ()
-  {
-    return *std::get< 0 >( constants_ );
-  }
   ConstantTupleType constants_;
   std::tuple<  > coefficients_;
 };
 
-} // namespace UFLLocalFunctions_d75993d7ae5919d23117b153e900840d
+} // namespace UFLLocalFunctions_3d957c5b2628e5120636b8e3f461bc0d
 
-PYBIND11_MODULE( localfunction_d75993d7ae5919d23117b153e900840d_a524c1196983e65de1c06d7d6afdeb44, module )
+PYBIND11_MODULE( localfunction_3d957c5b2628e5120636b8e3f461bc0d_a524c1196983e65de1c06d7d6afdeb44, module )
 {
-  typedef UFLLocalFunctions_d75993d7ae5919d23117b153e900840d::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_3d957c5b2628e5120636b8e3f461bc0d::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_d75993d7ae5919d23117b153e900840d::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_d75993d7ae5919d23117b153e900840d_a524c1196983e65de1c06d7d6afdeb44.cc"})).first;
+      auto cls = Dune::Python::insertClass<LocalFunctionType>(module,"UFLLocalFunction",Dune::Python::GenerateTypeName("UFLLocalFunctions_3d957c5b2628e5120636b8e3f461bc0d::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_3d957c5b2628e5120636b8e3f461bc0d_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( "c3", [] ( LocalFunctionType &self ) -> UFLLocalFunctions_d75993d7ae5919d23117b153e900840d::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 > > >::Conc3 { return self.conc3(); }, [] ( LocalFunctionType &self, const UFLLocalFunctions_d75993d7ae5919d23117b153e900840d::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 > > >::Conc3 &v ) { self.conc3() = v; } );
       cls.def_property_readonly( "virtualized", [] ( LocalFunctionType& ) -> bool { return true;});
   }
 }