13 August 2014 9:53:47.154 PM ELLIPSOID_MONTE_CARLO_PRB FORTRAN77 version Test the ELLIPSOID_MONTE_CARLO library. TEST01 Use ELLIPSOID_SAMPLE to estimate integrals in a 2D ellipse x' * A * x <= r^2. Ellipsoid radius R = 2.00000 Ellipsoid center V: 1: 0.0000000 2: 0.0000000 Ellipsoid matrix A: Col 1 2 Row 1: 9.00000 1.00000 2: 1.00000 4.00000 Ellipsoid volume = 2.12410 N 1 X Y X^2 XY Y^2 X^3 1 2.12410 0.455600 -0.376066 0.977220E-01 -0.806627E-01 0.665814E-01 0.209604E-01 2 2.12410 -0.306358E-01 0.147780 0.449434E-02 -0.307640E-02 0.105018E-01 -0.181719E-03 4 2.12410 0.275870 -0.312000 0.417762E-01 -0.231694E-01 0.140882 0.700437E-02 8 2.12410 0.333166E-02 -0.251031 0.251927E-01 -0.180584E-01 0.825508E-01 -0.312286E-03 16 2.12410 0.534942E-02 0.718986E-01 0.230367E-01 -0.985476E-02 0.117591 -0.846337E-03 32 2.12410 -0.576574E-01 -0.184163 0.317003E-01 -0.195772E-01 0.149220 -0.135138E-02 64 2.12410 -0.250151E-01 0.240724E-01 0.278674E-01 -0.233873E-01 0.145951 -0.251368E-05 128 2.12410 0.421042E-02 -0.629404E-01 0.285123E-01 -0.227768E-01 0.147640 0.484888E-04 256 2.12410 0.164242E-02 0.581420E-01 0.283226E-01 -0.196689E-01 0.135295 -0.829320E-04 512 2.12410 0.756784E-02 0.170093E-01 0.281605E-01 -0.224068E-01 0.141451 0.239302E-03 1024 2.12410 -0.767481E-02 -0.145921E-01 0.298340E-01 -0.230815E-01 0.149816 -0.185755E-03 2048 2.12410 -0.493584E-02 0.207855E-01 0.283452E-01 -0.212838E-01 0.140061 -0.484776E-04 4096 2.12410 -0.768025E-02 -0.135906E-02 0.294489E-01 -0.226947E-01 0.142611 -0.155983E-03 8192 2.12410 -0.663898E-03 0.852250E-02 0.296095E-01 -0.221527E-01 0.141227 -0.336923E-04 16384 2.12410 -0.849378E-03 0.544156E-02 0.295450E-01 -0.227768E-01 0.141865 0.318350E-04 32768 2.12410 -0.121951E-02 -0.488945E-03 0.296223E-01 -0.228369E-01 0.141868 -0.455365E-04 65536 2.12410 0.128864E-02 -0.146042E-02 0.296739E-01 -0.225971E-01 0.140368 0.210398E-04 TEST02 Use ELLIPSOID_SAMPLE to estimate integrals in a 2D ellipse (x-v)' * A * (x-v) <= r^2. Ellipsoid radius R = 2.00000 Ellipsoid center V: 1: 2.0000000 2: 3.0000000 Ellipsoid matrix A: Col 1 2 Row 1: 9.00000 1.00000 2: 1.00000 4.00000 Ellipsoid volume = 2.12410 N 1 X Y X^2 XY Y^2 X^3 1 2.12410 4.70381 5.99625 10.4165 13.2786 16.9271 23.0673 2 2.12410 4.21757 6.52009 8.37837 12.9452 20.0141 16.6520 4 2.12410 4.52408 6.06031 9.64167 12.9251 17.3858 20.5609 8 2.12410 4.25154 6.12128 8.53494 12.2345 17.6933 17.1837 16 2.12410 4.25356 6.44421 8.54085 12.8946 19.6659 17.1944 32 2.12410 4.19055 6.18815 8.29749 12.1838 18.1612 16.4898 64 2.12410 4.22319 6.39639 8.42422 12.6943 19.4073 16.8599 128 2.12410 4.25242 6.30937 8.54177 12.6086 18.8869 17.2145 256 2.12410 4.24985 6.43045 8.53131 12.8462 19.6011 17.1824 512 2.12410 4.25578 6.38932 8.55485 12.7789 19.3604 17.2529 1024 2.12410 4.24053 6.35772 8.49555 12.6693 19.1792 17.0796 2048 2.12410 4.24327 6.39310 8.50502 12.7501 19.3817 17.1036 4096 2.12410 4.24053 6.37095 8.49515 12.6962 19.2514 17.0772 8192 2.12410 4.24754 6.38084 8.52337 12.7375 19.3093 17.1625 16384 2.12410 4.24736 6.37775 8.52256 12.7302 19.2915 17.1599 32768 2.12410 4.24699 6.37182 8.52116 12.7172 19.2559 17.1559 65536 2.12410 4.24950 6.37085 8.53125 12.7230 19.2485 17.1864 TEST03 Use ELLIPSOID_SAMPLE to estimate integrals in a 3D ellipse (x-v)' * A * (x-v) <= r^2. Ellipsoid radius R = 0.500000 Ellipsoid center V: 1: 1.0000000 2: 2.0000000 3: 3.0000000 Ellipsoid matrix A: Col 1 2 3 Row 1: 9.00000 6.00000 3.00000 2: 6.00000 5.00000 4.00000 3: 3.00000 4.00000 9.00000 Ellipsoid volume = 0.872665E-01 N 1 X Y Z X^2 YZ Z^3 1 0.872665E-01 0.118556 0.127999 0.270935 0.161063 1.80968 2.61156 2 0.872665E-01 0.505427E-01 0.232324 0.249758 0.301574E-01 5.06528 2.04716 4 0.872665E-01 0.917626E-01 0.165736 0.265086 0.984727E-01 2.91490 2.45045 8 0.872665E-01 0.885784E-01 0.173441 0.261981 0.972310E-01 3.17841 2.37401 16 0.872665E-01 0.954258E-01 0.160572 0.265292 0.114340 2.81944 2.46932 32 0.872665E-01 0.815371E-01 0.184104 0.258630 0.844829E-01 3.48828 2.28260 64 0.872665E-01 0.873417E-01 0.173498 0.262541 0.947942E-01 3.18914 2.38842 128 0.872665E-01 0.908308E-01 0.169208 0.262950 0.103113 3.05812 2.40171 256 0.872665E-01 0.844516E-01 0.179415 0.260577 0.905178E-01 3.36271 2.33783 512 0.872665E-01 0.863446E-01 0.176134 0.261319 0.954505E-01 3.27468 2.35974 1024 0.872665E-01 0.868068E-01 0.175389 0.261589 0.954095E-01 3.24681 2.36496 2048 0.872665E-01 0.880433E-01 0.173281 0.262051 0.977763E-01 3.18198 2.37747 4096 0.872665E-01 0.870301E-01 0.174981 0.261639 0.957777E-01 3.23227 2.36634 8192 0.872665E-01 0.867427E-01 0.175376 0.261592 0.951869E-01 3.24487 2.36523 16384 0.872665E-01 0.875100E-01 0.174143 0.261904 0.967846E-01 3.20988 2.37371 32768 0.872665E-01 0.870833E-01 0.174799 0.261744 0.959327E-01 3.22908 2.36940 65536 0.872665E-01 0.873715E-01 0.174369 0.261825 0.965099E-01 3.21578 2.37158 ELLIPSOID_MONTE_CARLO_PRB Normal end of execution. 13 August 2014 9:53:47.408 PM