Merge pull request #13425 from lethe-cfd/add_erf_mu_parser

Add erf function to deal.II muparser functions
dealii · Feb 21, 2022 · e57f860 · e57f860
2 parents 790dc5d + 5a010da
commit e57f860
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 0 deletions.
diff --git a/include/deal.II/base/mu_parser_internal.h b/include/deal.II/base/mu_parser_internal.h
@@ -71,6 +71,9 @@ namespace internal
     double
     mu_pow(double a, double b);
 
+    double
+    mu_erf(double value);
+
     double
     mu_erfc(double value);
 

diff --git a/source/base/function_parser.cc b/source/base/function_parser.cc
@@ -187,6 +187,9 @@ FunctionParser<dim>::init_muparser() const
       fp.get()[component]->DefineFun("pow",
                                      internal::FunctionParser::mu_pow,
                                      true);
+      fp.get()[component]->DefineFun("erf",
+                                     internal::FunctionParser::mu_erf,
+                                     true);
       fp.get()[component]->DefineFun("erfc",
                                      internal::FunctionParser::mu_erfc,
                                      true);

diff --git a/source/base/mu_parser_internal.cc b/source/base/mu_parser_internal.cc
@@ -111,6 +111,12 @@ namespace internal
       return std::pow(a, b);
     }
 
+    double
+    mu_erf(double value)
+    {
+      return std::erf(value);
+    }
+
     double
     mu_erfc(double value)
     {
@@ -186,6 +192,7 @@ namespace internal
       "floor",
       "sec",
       "pow",
+      "erf",
       "erfc",
       "rand",
       "rand_seed"};