13 August 2014 4:11:51.952 PM ELLIPSOID_MONTE_CARLO_PRB FORTRAN90 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. 1. 2: 1. 4. 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 = 0.500000 Ellipsoid center V: 1: 2.0000000 2: 3.0000000 Ellipsoid matrix A: Col 1 2 Row 1: 9. 1. 2: 1. 4. Ellipsoid volume = 0.132757 N 1 X Y X^2 XY Y^2 X^3 1 0.132757 0.272632 0.392394 0.559883 0.805828 1.15981 1.14979 2 0.132757 0.265034 0.400579 0.529129 0.799709 1.20870 1.05641 4 0.132757 0.269824 0.393395 0.548431 0.799630 1.16611 1.11476 8 0.132757 0.265565 0.394347 0.531333 0.788780 1.17160 1.06327 16 0.132757 0.265597 0.399393 0.531450 0.798998 1.20201 1.06359 32 0.132757 0.264612 0.395392 0.527546 0.788005 1.17813 1.05198 64 0.132757 0.265122 0.398646 0.529571 0.796027 1.19764 1.05801 128 0.132757 0.265579 0.397286 0.531401 0.794681 1.18948 1.06351 256 0.132757 0.265539 0.399178 0.531239 0.798356 1.20079 1.06302 512 0.132757 0.265631 0.398535 0.531609 0.797338 1.19696 1.06413 1024 0.132757 0.265393 0.398042 0.530663 0.795633 1.19403 1.06131 2048 0.132757 0.265436 0.398594 0.530828 0.796874 1.19730 1.06179 4096 0.132757 0.265393 0.398248 0.530661 0.796048 1.19524 1.06130 8192 0.132757 0.265503 0.398403 0.531100 0.796688 1.19616 1.06262 16384 0.132757 0.265500 0.398355 0.531088 0.796580 1.19587 1.06259 32768 0.132757 0.265494 0.398262 0.531066 0.796377 1.19532 1.06252 65536 0.132757 0.265533 0.398247 0.531223 0.796466 1.19522 1.06299 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. 6. 3. 2: 6. 5. 4. 3: 3. 4. 9. 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 4:11:52.301 PM