test nz

Author: tk-yoshimura

.gitignore

@@ -1,13 +1,54 @@
-/.vs
+# User-specific files
+*.rsuser
+*.suo
+*.user
+*.userosscache
+*.sln.docstates
+*.userprefs
+
+# Images
+*.png
+*.jpg
+*.jpeg
+*.bmp
+*.webp
+
+# WinMarge backups
+*.bak
+
+# Cache/Options directory
+/.vs
 /dll
-/packages
-/**/bin
-/**/obj
-/**/*.png
-/**/*.jpg
-/**/*.bak
-/**/*.user
-/results_disused
-/TestResults
+
+# NuGet Packages
+*.nupkg
+*.snupkg
+/**/[Pp]ackages/*
+
+# Build results
+/**/[Bb]in/
+/**/[Oo]bj/
+/**/[Ll]og/
+/**/[Ll]ogs/
+/**/bld/
+
+# MSTest test Results
+[Tt]est[Rr]esult*/
+[Bb]uild[Ll]og.*
+
+# Publish
+/**/publish/
+/**/PublishProfiles
+*.[Pp]ublish.xml
+*.azurePubxml
+*.pubxml
+*.publishproj
+
+# Python cache
 /**/__pycache__
-/**/PublishProfiles
+*.pyc
+
+# Other
+/results_disused
+
+# Safelist

ErrorFunctionApproximation/Program.cs

@@ -63,6 +63,21 @@ static void Main(string[] args) {
                     sw.WriteLine($"{z},{y_true},{y_approx},{mp_err},{double_err}");
                 }
             }
+
+            using (StreamWriter sw = new StreamWriter("../../../../results/inverf_approx_nz.csv")) {
+                sw.WriteLine("z,y_true,y_approx,mp_error,double_error");
+
+                for (int i = -100; i <= +100; i++) {
+                    double z = 1.220703125e-4 + i / 1000000.0;
+
+                    MultiPrecision<Pow2.N4> y_true = MultiPrecision<Pow2.N4>.InverseErf(z);
+                    double y_approx = ErrorFunction.InverseErf(z);
+                    MultiPrecision<Pow2.N4> mp_err = y_true - y_approx;
+                    double double_err = (double)y_true - y_approx;
+
+                    sw.WriteLine($"{z},{y_true},{y_approx},{mp_err},{double_err}");
+                }
+            }
         }
     }
 }

results/inverf_approx.csv

@@ -1,5 +1,5 @@
 z,y_true,y_approx,mp_error,double_error
--1,-inf,-inf,NaN,NaN
+-1,-∞,-∞,NaN,NaN
 -0.9999,-2.7510639057120796917425935925819583,-2.75106390571208,1.2391489564199963139825013762384876e-16,0
 -0.9998,-2.6297417762102925939779247515988718,-2.6297417762102926,-1.9803387646251539898754354370091143e-18,0
 -0.9997,-2.5564031509202212422794006041040893,-2.556403150920221,-5.2244946919756407637267832582595266e-17,0
@@ -19999,4 +19999,4 @@ z,y_true,y_approx,mp_error,double_error
 0.9997,2.5564031509202212422794006041040893,2.556403150920221,5.2244946919756407637267832582595266e-17,0
 0.9998,2.6297417762102925939779247515988718,2.6297417762102926,1.9803387646251539898754354370091143e-18,0
 0.9999,2.7510639057120796917425935925819583,2.75106390571208,-1.2391489564199963139825013762384876e-16,0
-1,inf,inf,NaN,NaN
+1,∞,∞,NaN,NaN