--- FGM(0.5) (d=2), grid nx=1024, ny=4000^1
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.30000  ||H-P_N H||~0.04167  bound~1.04167  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=2.19e-01  D1(A,B_N)~0.15338  ||H-P_N H||~0.02083  bound~0.52083  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=4.49e-02  D1(A,B_N)~0.07700  ||H-P_N H||~0.01042  bound~0.26042  OK
N=  8  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=9.64e-03  D1(A,B_N)~0.03853  ||H-P_N H||~0.00521  bound~0.13021  OK
N= 16  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=2.19e-03  D1(A,B_N)~0.01927  ||H-P_N H||~0.00260  bound~0.06510  OK
N= 32  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=5.18e-04  D1(A,B_N)~0.00963  ||H-P_N H||~0.00130  bound~0.03255  OK
N= 64  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.26e-04  D1(A,B_N)~0.00482  ||H-P_N H||~0.00065  bound~0.01628  OK
--- (M+Pi)/2 (d=2), grid nx=1024, ny=4000^1
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.16667  ||H-P_N H||~0.16667  bound~1.16667  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.25e-01  D1(A,B_N)~0.08333  ||H-P_N H||~0.08333  bound~0.58333  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=3.12e-02  D1(A,B_N)~0.04167  ||H-P_N H||~0.04167  bound~0.29167  OK
N=  8  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=7.81e-03  D1(A,B_N)~0.02083  ||H-P_N H||~0.02083  bound~0.14583  OK
N= 16  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.95e-03  D1(A,B_N)~0.01042  ||H-P_N H||~0.01041  bound~0.07291  OK
N= 32  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=4.88e-04  D1(A,B_N)~0.00521  ||H-P_N H||~0.00520  bound~0.03645  OK
N= 64  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.22e-04  D1(A,B_N)~0.00261  ||H-P_N H||~0.00260  bound~0.01823  OK
--- Clayton(2) (d=2), grid nx=1024, ny=4000^1
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.17157  ||H-P_N H||~0.18007  bound~1.18007  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.22e-01  D1(A,B_N)~0.12523  ||H-P_N H||~0.09442  bound~0.59442  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=5.86e-03  D1(A,B_N)~0.06843  ||H-P_N H||~0.04754  bound~0.29754  OK
N=  8  total=1.000000000000  max|margin-1/N|=1.1e-16  minmass=2.98e-04  D1(A,B_N)~0.03522  ||H-P_N H||~0.02368  bound~0.14868  OK
N= 16  total=1.000000000000  max|margin-1/N|=2.2e-16  minmass=1.68e-05  D1(A,B_N)~0.01782  ||H-P_N H||~0.01179  bound~0.07429  OK
N= 32  total=1.000000000000  max|margin-1/N|=1.1e-16  minmass=1.00e-06  D1(A,B_N)~0.00896  ||H-P_N H||~0.00588  bound~0.03713  OK
N= 64  total=1.000000000000  max|margin-1/N|=2.2e-16  minmass=6.10e-08  D1(A,B_N)~0.00449  ||H-P_N H||~0.00293  bound~0.01856  OK
--- FGM3(0.5) (d=3), grid nx=256, ny=256^2
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.31208  ||H-P_N H||~0.00694  bound~2.00694  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.17e-01  D1(A,B_N)~0.14978  ||H-P_N H||~0.00347  bound~1.00347  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.23e-02  D1(A,B_N)~0.07401  ||H-P_N H||~0.00174  bound~0.50174  OK
N=  8  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.30e-03  D1(A,B_N)~0.03689  ||H-P_N H||~0.00087  bound~0.25087  OK
N= 16  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.44e-04  D1(A,B_N)~0.01845  ||H-P_N H||~0.00043  bound~0.12543  OK
N= 32  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.66e-05  D1(A,B_N)~0.00927  ||H-P_N H||~0.00022  bound~0.06272  OK
--- (M3+Pi3)/2 (d=3), grid nx=256, ny=256^2
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.15882  ||H-P_N H||~0.16732  bound~2.16732  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=6.25e-02  D1(A,B_N)~0.07631  ||H-P_N H||~0.08365  bound~1.08365  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=7.81e-03  D1(A,B_N)~0.03770  ||H-P_N H||~0.04182  bound~0.54182  OK
N=  8  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=9.77e-04  D1(A,B_N)~0.01879  ||H-P_N H||~0.02089  bound~0.27089  OK
N= 16  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.22e-04  D1(A,B_N)~0.00940  ||H-P_N H||~0.01042  bound~0.13542  OK
N= 32  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.53e-05  D1(A,B_N)~0.00472  ||H-P_N H||~0.00515  bound~0.06765  OK
--- Pi3 (d=3), grid nx=256, ny=256^2
N=  1  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.00e+00  D1(A,B_N)~0.31765  ||H-P_N H||~0.00000  bound~2.00000  OK
N=  2  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.25e-01  D1(A,B_N)~0.15262  ||H-P_N H||~0.00000  bound~1.00000  OK
N=  4  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.56e-02  D1(A,B_N)~0.07540  ||H-P_N H||~0.00000  bound~0.50000  OK
N=  8  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=1.95e-03  D1(A,B_N)~0.03759  ||H-P_N H||~0.00000  bound~0.25000  OK
N= 16  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=2.44e-04  D1(A,B_N)~0.01880  ||H-P_N H||~0.00000  bound~0.12500  OK
N= 32  total=1.000000000000  max|margin-1/N|=0.0e+00  minmass=3.05e-05  D1(A,B_N)~0.00944  ||H-P_N H||~0.00000  bound~0.06250  OK
