May 17 2007 11:02:35.851 AM TEST_TRIANGULATION_PRB FORTRAN90 version Test the routines in the TEST_TRIANGULATION library. TEST01 P00_TEST_NUM reports the number of problems. P00_TITLE returns a title for each problem. P00_HEADER prints some information about each problem. The number of tests available = 15 Test number: 1 Title: "#1: The unit circle." P01: Strang and Persson example #1 The unit circle. Radius = 1.00000 Center = .000000 .000000 A uniform mesh density is requested. Element sizes tried were 0.4, 0.2, 0.1. Number of boundary segments = 1 Number of fixed points = 0 Number of holes = 0 Test number: 2 Title: "#2: The circle with concentric circular hole." P02: Strang and Persson example #2 The unit circle, with a concentric hole. Radius1 = 1.00000 Radius2 = .400000 Center = .000000 .000000 A uniform mesh density is requested. Element sizes tried were 0.4, 0.2, 0.1. Number of boundary segments = 2 Number of fixed points = 0 Number of holes = 1 Test number: 3 Title: "#3: The unit square with circular hole." P03: Strang and Persson example #3 The unit square, with a hole. The hole is a concentric circle of radius 0.4. A uniform mesh density is requested. Element sizes tried were 0.4, 0.2, 0.1. Number of boundary segments = 2 Number of fixed points = 4 Number of holes = 1 Test number: 4 Title: "#4: The unit hexagon with hexagonal hole." P04: Strang and Persson example #4 The hexagon with hexagonal hole. Radius of outer hexagon R1 = 1.00000 Radius of outer hexagon R2 = .500000 A uniform mesh density is requested. Element sizes tried were ? Number of boundary segments = 2 Number of fixed points = 12 Number of holes = 1 Test number: 5 Title: "#5: The horn." P05: Strang and Persson example #5 The horn. Circle C1 has center = (0,0) Radius R1 = 1.00000 Circle C2 has center = (-0.4,0) Radius R2 = .550000 Points in the region are: in C1 and not in C2 and have 0 <= Y. A uniform mesh density is requested. Element sizes tried were 0.4, 0.2, 0.1. Number of boundary segments = 1 Number of fixed points = 4 Number of holes = 0 Test number: 6 Title: "#6: The superellipse with superelliptical hole." P06: Strang and Persson example #6 Superellipse with superelliptical hole. Radius R1 = 1.00000 Radius R2 = .500000 A uniform mesh density is requested. Element sizes tried were 0.4, 0.2, 0.1. Number of boundary segments = 2 Number of fixed points = 0 Number of holes = 1 Test number: 7 Title: "#7: Bicycle seat (implicit)." P07: Strang and Persson example #7 Bicycle seat (implicit). A uniform mesh density is requested. The boundary is formed by two algebraic expressions. Number of boundary segments = 1 Number of fixed points = 2 Number of holes = 0 Test number: 8 Title: "#8: Pie slice with notch and hole." P08: Strang and Persson example #8 Pie slice with notch and hole. The pie rim is a portion of a circle C1 with CENTER1 = .000000 .000000 and radius R1 = 1.00000 The interior hole is a circle C2 with CENTER2 = .600000 .000000 and radius R2 = .100000 A uniform mesh density is requested. Number of boundary segments = 2 Number of fixed points = 6 Number of holes = 1 Test number: 9 Title: "#9: Jeff Borggaard's Box with 2 hexagonal holes." P09: Jeff Borggaard's example A square with 2 hexagonal holes. The square has "center" at .500000 .500000 and "radius" R1 = .500000 Hexagon 1 has "center" at .250000 .750000 and "radius" R2 = .100000 Hexagon 2 has "center" at .600000 .400000 and "radius" R3 = .100000 A uniform mesh density is requested. Number of boundary segments = 3 Number of fixed points = 16 Number of holes = 2 Test number: 10 Title: "#10: The unit square." P10: The unit square. Number of boundary segments = 1 Number of fixed points = 4 Number of holes = 0 Test number: 11 Title: "#11: The L-shaped region." P11: The L-shaped region. The lower left corner of the indentation is P = .500000 .500000 Number of boundary segments = 1 Number of fixed points = 6 Number of holes = 0 Test number: 12 Title: "#12: John Shadid's H-shaped region." P12: John Shadid's H-shaped region. Number of boundary segments = 1 Number of fixed points = 12 Number of holes = 0 Test number: 13 Title: "#13: The Sandia Fork region." P13: The Sandia Fork region. Number of boundary segments = 1 Number of fixed points = 8 Number of holes = 0 Test number: 14 Title: "#14: Marcus Garvie's Lake Alpha with Beta Island." P14: Marcus Garvie's Lake Alpha and Beta island. Number of boundary segments = 2 Number of fixed points = 101 Number of holes = 1 Test number: 15 Title: "#15: Sangbum Kim's forward step." P15: Sangbum Kim's forward step. Number of boundary segments = 1 Number of fixed points = 8 Number of holes = 0 TEST02 P00_TEST_NUM reports the number of problems. P00_SAMPLE returns sample points from the region. Test number = 1 Title: = "#1: The unit circle." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .182719 .908686 .882655 -.248530 .168080 -.307982 -0.802429E-01 -0.327634E-01 -.439650 .258865 .873988 .385525 -0.160461E-01 .336759 .457004 .377857 .872951 .246665 -.344310 -.385464 .770117 .314621 -.712087 .233654 -.755233 .539294 0.197782E-01 -.673352 .279579 -.913133 .773199 0.893149E-02 .347621 -.260987 -.516254 .703544 .522130 .352660 .428907 0.367860E-01 Test number = 2 Title: = "#2: The circle with concentric circular hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .185100 .920532 .895921 -.252265 .245863 -.450509 -.377553 -.154156 -.530260 .312216 .880668 .388471 -0.240566E-01 .504875 .520067 .429999 .887855 .250877 -.413022 -.462389 .797059 .325628 -.755239 .247813 -.764908 .546203 0.215989E-01 -.735337 .281730 -.920157 .813733 0.939971E-02 .451474 -.338959 -.529031 .720957 .582130 .393186 .559785 0.480110E-01 Test number = 3 Title: = "#3: The unit square with circular hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col -.563163 .912635 .659018 .123391 -.169386 -.867763 -.484844 -.780086 -.912342 .267931 -.876546 -.100922 -.197387 .509347 .594574 -.996323 .795008 -.298495 -.810910 -.972766 .718194 .681695 -.753792 -.984975 -.479394 .824967 -.772672 -.296743 .645775 -.465735 .384133 .123325 .722432 -0.924124E-01 .823954 .195834 -.622091 .522984 -.206023 -.629372 Test number = 4 Title: = "#4: The unit hexagon with hexagonal hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .659018 .123391 -.484844 -.780086 -.876546 -.100922 -.197387 .509347 .795008 -.298495 -.479394 .824967 -.772672 -.296743 .645775 -.465735 .722432 -0.924124E-01 .823954 .195834 -.622091 .522984 -.206023 -.629372 -.574140 .428942 -.764586 -.401342 .490606 -.382650 .523461 -.186061 .877499 .124177 -.964360 0.220645E-02 -.456553 .717145 -.292185 -.761384 Test number = 5 Title: = "#5: The horn." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .659018 .561695 -.876546 .449539 -.197387 .754673 .594574 0.183837E-02 .795008 .350752 .645775 .267132 .384133 .561662 .722432 .453794 -.622091 .761492 .148732 .367027 .234410 .361529 -.574140 .714471 .490606 .308675 .523461 .406970 -.916181 .368851 -.456553 .858573 .225895 .809519 .174179 .215492 .536113 .723297 -.103962 .855176 Test number = 6 Title: = "#6: The superellipse with superelliptical hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col -.563163 .912635 .659018 .123391 -.169386 -.867763 -.484844 -.780086 -.912342 .267931 -.876546 -.100922 -.197387 .509347 .795008 -.298495 .718194 .681695 -.479394 .824967 -.772672 -.296743 .645775 -.465735 .722432 -0.924124E-01 .823954 .195834 -.622091 .522984 -.206023 -.629372 -.574140 .428942 -.764586 -.401342 .650006 .649320 -.876276 .421561 Test number = 7 Title: = "#7: Bicycle seat (implicit)." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col 5.17592 -1.62983 -1.33035 -4.60329 -3.80796 -4.34026 -7.16552 -1.19621 -1.55028 -.471959 6.24398 -2.89549 5.64068 0.450847E-01 -6.06855 -2.89023 5.07190 -3.39721 3.01697 -1.63003 5.67396 -2.27724 6.47132 -1.41250 -4.88589 -.431048 -1.61810 -3.88812 1.16814 -2.79784 1.84105 -2.83083 -4.50929 -.713173 -6.00505 -3.20403 5.10513 -0.520396E-01 -6.88226 -.735317 Test number = 8 Title: = "#8: Pie slice with notch and hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .829509 0.319359E-01 .897504 -0.772563E-01 .859097 .176436 .822887 -.120541 .861216 -0.239181E-01 .911977 0.506855E-01 .825003 .168056 .745303 -0.990371E-01 .899373 .136417 .761731 -0.481561E-01 0.178200E-01 0.571071E-03 .612948 .160219 .587090 -.147272 .768056 .115587 .945017 .211743 .481067 -.107507 .902640 0.868812E-01 .740790 -0.227249E-01 .775651 -0.712986E-01 .563260 .105871 Test number = 9 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .218418 .956318 .829509 .561695 .415307 0.661187E-01 .257578 .109957 0.438290E-01 .633966 0.617272E-01 .449539 .401306 .754673 .797287 0.183837E-02 .897504 .350752 0.945448E-01 0.136169E-01 .859097 .840847 .123104 0.751236E-02 .260303 .912484 .113664 .351629 .822887 .267132 .692066 .561662 .861216 .453794 .911977 .597917 .396988 .185314 .117707 .299329 Test number = 10 Title: = "#10: The unit square." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .218418 .956318 .829509 .561695 .415307 0.661187E-01 .257578 .109957 0.438290E-01 .633966 0.617272E-01 .449539 .401306 .754673 .797287 0.183837E-02 .897504 .350752 0.945448E-01 0.136169E-01 .859097 .840847 .123104 0.751236E-02 .260303 .912484 .113664 .351629 .822887 .267132 .692066 .561662 .861216 .453794 .911977 .597917 .188955 .761492 .396988 .185314 Test number = 11 Title: = "#11: The L-shaped region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .109209 .956318 .414755 .561695 .207654 0.661187E-01 .128789 .109957 0.219145E-01 .633966 .530864 .224769 .200653 .754673 .398643 0.183837E-02 .948752 .175376 .547272 0.680845E-02 .429548 .840847 .561552 0.375618E-02 .130151 .912484 .556832 .175814 .911444 .133566 .346033 .561662 .930608 .226897 .955989 .298958 .594477 .380746 .198494 .185314 Test number = 12 Title: = "#12: John Shadid's H-shaped region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .527302 .369540 .414755 .561695 .780740 0.661187E-01 .128789 .109957 0.219145E-01 .633966 .648148 .449539 .200653 .754673 .398643 0.183837E-02 .961564 .350752 .660454 0.136169E-01 .429548 .840847 .671164 0.751236E-02 .130151 .912484 .667624 .351629 .933583 .267132 .346033 .561662 .947956 .453794 .966991 .597917 .695858 .761492 .198494 .185314 Test number = 13 Title: = "#13: The Sandia Fork region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col 12.3104 .751236 57.4366 36.7027 61.7205 36.1529 74.5303 30.8675 15.2384 11.4319 20.6653 21.2924 54.0539 83.2281 53.1111 37.9977 49.0323 85.4159 84.8580 8.31430 56.5026 39.4392 53.8036 77.5243 50.9806 30.5324 52.9502 34.4293 53.0101 40.6241 46.6644 87.8359 68.4731 27.8111 67.9375 26.1012 68.8186 33.7421 42.4428 36.6481 Test number = 14 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col 321.774 188.242 237.547 216.912 314.298 638.556 239.002 741.764 160.697 374.965 539.422 319.705 469.564 512.327 559.889 441.781 311.992 266.195 406.711 385.035 429.587 381.440 162.855 340.761 497.992 346.873 506.764 411.158 245.100 706.506 288.986 223.027 210.353 284.252 356.890 336.173 320.156 247.644 495.582 443.289 Test number = 15 Title: = "#15: Sangbum Kim's forward step." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 The sample points: Row 1 2 Col .184183 .956318 .295092 -.438305 2.15307 0.661187E-01 -5.42422 -.890043 -7.56171 -.366034 -1.38273 .449539 -3.98694 -.245327 -0.271305E-01 -.998162 6.97504 .350752 -1.05455 0.136169E-01 .590969 -.159153 -.768961 0.751236E-02 -5.39697 -0.875163E-01 -.863360 .351629 6.22887 .267132 -1.07934 -.438338 6.61216 .453794 7.11977 .597917 -.110453 .761492 -4.03012 -.814686 TEST03 P00_BOUNDARY_NEAREST returns the nearest boundary point for a set of points. Test number = 1 Title: = "#1: The unit circle." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .182719 .908686 .197134 .980377 .882655 -.248530 .962570 -.271032 .168080 -.307982 .479049 -.877788 -0.802429E-01 -0.327634E-01 -.925802 -.378008 -.439650 .258865 -.861723 .507380 .873988 .385525 .914940 .403589 -0.160461E-01 .336759 -0.475947E-01 .998867 .457004 .377857 .770686 .637215 .872951 .246665 .962320 .271918 -.344310 -.385464 -.666173 -.745798 .770117 .314621 .925727 .378193 -.712087 .233654 -.950157 .311771 -.755233 .539294 -.813814 .581125 0.197782E-01 -.673352 0.293602E-01 -.999569 .279579 -.913133 .292761 -.956186 .773199 0.893149E-02 .999933 0.115506E-01 .347621 -.260987 .799700 -.600400 -.516254 .703544 -.591603 .806229 .522130 .352660 .828685 .559715 .428907 0.367860E-01 .996342 0.854531E-01 Test number = 2 Title: = "#2: The circle with concentric circular hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .185100 .920532 .197134 .980377 .895921 -.252265 .962570 -.271032 .245863 -.450509 .191620 -.351115 -.377553 -.154156 -.370321 -.151203 -.530260 .312216 -.344689 .202952 .880668 .388471 .914940 .403589 -0.240566E-01 .504875 -0.190379E-01 .399547 .520067 .429999 .308275 .254886 .887855 .250877 .962320 .271918 -.413022 -.462389 -.266469 -.298319 .797059 .325628 .925727 .378193 -.755239 .247813 -.950157 .311771 -.764908 .546203 -.813814 .581125 0.215989E-01 -.735337 0.293602E-01 -.999569 .281730 -.920157 .292761 -.956186 .813733 0.939971E-02 .999933 0.115506E-01 .451474 -.338959 .319880 -.240160 -.529031 .720957 -.591603 .806229 .582130 .393186 .828685 .559715 .559785 0.480110E-01 .398537 0.341812E-01 Test number = 3 Title: = "#3: The unit square with circular hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) -.563163 .912635 -.563163 1.00000 .659018 .123391 .393168 0.736145E-01 -.169386 -.867763 -.169386 -1.00000 -.484844 -.780086 -.484844 -1.00000 -.912342 .267931 -1.00000 .267931 -.876546 -.100922 -1.00000 -.100922 -.197387 .509347 -.144538 .372973 .594574 -.996323 .594574 -1.00000 .795008 -.298495 1.00000 -.298495 -.810910 -.972766 -.810910 -1.00000 .718194 .681695 1.00000 .681695 -.753792 -.984975 -.753792 -1.00000 -.479394 .824967 -.479394 1.00000 -.772672 -.296743 -1.00000 -.296743 .645775 -.465735 1.00000 -.465735 .384133 .123325 .380854 .122272 .722432 -0.924124E-01 1.00000 -0.924124E-01 .823954 .195834 1.00000 .195834 -.622091 .522984 -1.00000 .522984 -.206023 -.629372 -.124441 -.380150 Test number = 4 Title: = "#4: The unit hexagon with hexagonal hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .659018 .123391 .433013 .123391 -.484844 -.780086 -.533424 -.808134 -.876546 -.100922 -.925436 -.129149 -.197387 .509347 -.143993 .416865 .795008 -.298495 .819500 -.312636 -.479394 .824967 -.512627 .844155 -.772672 -.296743 -.814675 -.320993 .645775 -.465735 .709774 -.502686 .722432 -0.924124E-01 .890592 -.189500 .823954 .195834 .871190 .223105 -.622091 .522984 -.679064 .555878 -.206023 -.629372 -0.984977E-01 -.443132 -.574140 .428942 -.707797 .506110 -.764586 -.401342 -.767361 -.402943 .490606 -.382650 .418769 -.258224 .523461 -.186061 .433013 -.186061 .877499 .124177 .915605 .146177 -.964360 0.220645E-02 -.990135 0.170874E-01 -.456553 .717145 -.553605 .773178 -.292185 -.761384 -.292185 -.866025 Test number = 5 Title: = "#5: The horn." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .659018 .561695 .761067 .648673 -.876546 .449539 -.889806 .456339 -.197387 .754673 -.140705 .965798 .594574 0.183837E-02 .594574 .000000 .795008 .350752 .914912 .403653 .645775 .267132 .568890 .247493 .384133 .561662 .412964 .582314 .722432 .453794 .846798 .531914 -.622091 .761492 -.632660 .774430 .148732 .367027 .431207 .555964 .234410 .361529 .468827 .495115 -.574140 .714471 -.626399 .779503 .490606 .308675 .544859 .327478 .523461 .406970 .515079 .403275 -.916181 .368851 -.927644 .373465 -.456553 .858573 -.469505 .882930 .225895 .809519 .211666 .791116 .174179 .215492 .174179 .000000 .536113 .723297 .595469 .803378 -.103962 .855176 -0.728742E-01 .944981 Test number = 6 Title: = "#6: The superellipse with superelliptical hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) -.563163 .912635 -.575375 .971397 .659018 .123391 .499569 .121119 -.169386 -.867763 -.170035 -.999791 -.484844 -.780086 -.513765 -.982107 -.912342 .267931 -.998676 .269630 -.876546 -.100922 -.999974 -.101050 -.197387 .509347 -.196622 .496984 .795008 -.298495 .997851 -.304245 .718194 .681695 .874616 .802549 -.479394 .824967 -.500292 .983956 -.772672 -.296743 -.997885 -.303051 .645775 -.465735 .458239 -.368338 .722432 -0.924124E-01 .499862 -0.910665E-01 .823954 .195834 .999622 .197182 -.622091 .522984 -.442106 .394807 -.206023 -.629372 -.197684 -.496917 -.574140 .428942 -.458130 .368547 -.764586 -.401342 -.992248 -.418412 .650006 .649320 .841973 .839815 -.876276 .421561 -.991258 .431013 Test number = 7 Title: = "#7: Bicycle seat (implicit)." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 5.17592 -1.62983 4.19253 -.496754 -1.33035 -4.60329 -1.33527 -4.99582 -3.80796 -4.34026 -3.91862 -4.69015 -7.16552 -1.19621 -7.30254 -1.26315 -1.55028 -.471959 -1.79433 -.221678 6.24398 -2.89549 6.29605 -2.93518 5.64068 0.450847E-01 5.24306 .506111 -6.06855 -2.89023 -6.23161 -3.01842 5.07190 -3.39721 5.44951 -3.84112 3.01697 -1.63003 3.06525 -.997087 5.67396 -2.27724 6.39528 -2.80189 6.47132 -1.41250 7.06206 -1.73160 -4.88589 -.431048 -4.58468 -.127364 -1.61810 -3.88812 -1.64384 -4.99040 1.16814 -2.79784 1.18750 -4.99739 1.84105 -2.83083 2.68787 -.898822 -4.50929 -.713173 -4.26052 -.436650 -6.00505 -3.20403 -6.04844 -3.24133 5.10513 -0.520396E-01 4.88483 .171585 -6.88226 -.735317 -7.43237 -.990232 Test number = 8 Title: = "#8: Pie slice with notch and hole." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .829509 0.319359E-01 .900000 .000000 .897504 -0.772563E-01 .937380 -0.373802E-01 .859097 .176436 .845657 .226593 .822887 -.120541 .797900 -.213797 .861216 -0.239181E-01 .900000 .000000 .911977 0.506855E-01 .931331 0.313313E-01 .825003 .168056 .811752 .217508 .745303 -0.990371E-01 .682632 -0.563208E-01 .899373 .136417 .988691 .149964 .761731 -0.481561E-01 .695842 -0.285373E-01 0.178200E-01 0.571071E-03 0.167691E-01 0.449326E-02 .612948 .160219 .611943 .163970 .587090 -.147272 .584580 -.156638 .768056 .115587 .745503 .199757 .945017 .211743 .975805 .218642 .481067 -.107507 .475718 -.127468 .902640 0.868812E-01 .944761 0.447605E-01 .740790 -0.227249E-01 .698722 -0.159348E-01 .775651 -0.712986E-01 .692658 -0.376108E-01 .563260 .105871 .567215 0.944730E-01 Test number = 9 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .218418 .956318 .218418 1.00000 .829509 .561695 1.00000 .561695 .415307 0.661187E-01 .415307 .000000 .257578 .109957 .257578 .000000 0.438290E-01 .633966 .000000 .633966 0.617272E-01 .449539 .000000 .449539 .401306 .754673 .350000 .750000 .797287 0.183837E-02 .797287 .000000 .897504 .350752 1.00000 .350752 0.945448E-01 0.136169E-01 0.945448E-01 .000000 .859097 .840847 1.00000 .840847 .123104 0.751236E-02 .123104 .000000 .260303 .912484 .260303 .836603 .113664 .351629 .000000 .351629 .822887 .267132 .673188 .353561 .692066 .561662 .650000 .486603 .861216 .453794 1.00000 .453794 .911977 .597917 1.00000 .597917 .396988 .185314 .396988 .000000 .117707 .299329 .000000 .299329 Test number = 10 Title: = "#10: The unit square." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .218418 .956318 .218418 1.00000 .829509 .561695 1.00000 .561695 .415307 0.661187E-01 .415307 .000000 .257578 .109957 .257578 .000000 0.438290E-01 .633966 .000000 .633966 0.617272E-01 .449539 .000000 .449539 .401306 .754673 .350000 .750000 .797287 0.183837E-02 .797287 .000000 .897504 .350752 1.00000 .350752 0.945448E-01 0.136169E-01 0.945448E-01 .000000 .859097 .840847 1.00000 .840847 .123104 0.751236E-02 .123104 .000000 .260303 .912484 .260303 .836603 .113664 .351629 .000000 .351629 .822887 .267132 1.00000 .267132 .692066 .561662 1.00000 .561662 .861216 .453794 1.00000 .453794 .911977 .597917 1.00000 .597917 .188955 .761492 .000000 .761492 .396988 .185314 .396988 .000000 Test number = 11 Title: = "#11: The L-shaped region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .109209 .956318 .109209 1.00000 .414755 .561695 .500000 .561695 .207654 0.661187E-01 .207654 .000000 .128789 .109957 .128789 .000000 0.219145E-01 .633966 .000000 .633966 .530864 .224769 .530864 .000000 .200653 .754673 .000000 .754673 .398643 0.183837E-02 .398643 .000000 .948752 .175376 1.00000 .175376 .547272 0.680845E-02 .547272 .000000 .429548 .840847 .500000 .840847 .561552 0.375618E-02 .561552 .000000 .130151 .912484 .130151 1.00000 .556832 .175814 .556832 .000000 .911444 .133566 1.00000 .133566 .346033 .561662 .500000 .561662 .930608 .226897 1.00000 .226897 .955989 .298958 1.00000 .298958 .594477 .380746 .594477 .500000 .198494 .185314 .198494 .000000 Test number = 12 Title: = "#12: John Shadid's H-shaped region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .527302 .369540 .527302 .375000 .414755 .561695 .500000 .561695 .780740 0.661187E-01 .780740 .000000 .128789 .109957 .128789 .000000 0.219145E-01 .633966 .000000 .633966 .648148 .449539 .625000 .449539 .200653 .754673 .000000 .754673 .398643 0.183837E-02 .398643 .000000 .961564 .350752 1.00000 .350752 .660454 0.136169E-01 .660454 .000000 .429548 .840847 .500000 .840847 .671164 0.751236E-02 .671164 .000000 .130151 .912484 .130151 1.00000 .667624 .351629 .625000 .375000 .933583 .267132 1.00000 .267132 .346033 .561662 .500000 .561662 .947956 .453794 1.00000 .453794 .966991 .597917 1.00000 .597917 .695858 .761492 .625000 .761492 .198494 .185314 .198494 .000000 Test number = 13 Title: = "#13: The Sandia Fork region." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 12.3104 .751236 10.0079 .797129 57.4366 36.7027 57.9433 39.2034 61.7205 36.1529 62.3356 38.0504 74.5303 30.8675 74.8864 31.3156 15.2384 11.4319 12.0020 12.4962 20.6653 21.2924 17.6286 23.4966 54.0539 83.2281 55.0000 83.2281 53.1111 37.9977 55.0000 39.6863 49.0323 85.4159 45.0000 85.4159 84.8580 8.31430 88.9085 9.28043 56.5026 39.4392 56.5072 39.4672 53.8036 77.5243 55.0000 77.5243 50.9806 30.5324 50.9630 29.9845 52.9502 34.4293 52.5613 29.8905 53.0101 40.6241 55.0000 40.6241 46.6644 87.8359 45.0000 87.8359 68.4731 27.8111 66.5989 24.9895 67.9375 26.1012 66.9913 24.7244 68.8186 33.7421 69.4834 34.9342 42.4428 36.6481 41.9216 39.1758 Test number = 14 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 321.774 188.242 326.370 184.246 237.547 216.912 221.120 207.903 314.298 638.556 315.682 640.567 239.002 741.764 238.696 742.068 160.697 374.965 141.043 388.608 539.422 319.705 575.416 328.265 469.564 512.327 490.386 556.714 559.889 441.781 565.637 441.320 311.992 266.195 283.304 277.170 406.711 385.035 400.994 328.424 429.587 381.440 412.591 327.252 162.855 340.761 125.723 366.538 497.992 346.873 475.441 290.688 506.764 411.158 558.007 406.866 245.100 706.506 266.257 712.110 288.986 223.027 270.685 255.554 210.353 284.252 235.140 287.952 356.890 336.173 375.550 319.158 320.156 247.644 357.246 251.294 495.582 443.289 562.901 426.550 Test number = 15 Title: = "#15: Sangbum Kim's forward step." Spatial dimension M = 2 Number of samples N = 20 Initial SEED: = 123456789 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) .184183 .956318 -8.00000 .956318 .295092 -.438305 -8.00000 -.438305 2.15307 0.661187E-01 -8.00000 0.661187E-01 -5.42422 -.890043 -8.00000 -.890043 -7.56171 -.366034 -8.00000 -.366034 -1.38273 .449539 -8.00000 .449539 -3.98694 -.245327 -8.00000 -.245327 -0.271305E-01 -.998162 -8.00000 -.998162 6.97504 .350752 -8.00000 .350752 -1.05455 0.136169E-01 -8.00000 0.136169E-01 .590969 -.159153 -8.00000 -.159153 -.768961 0.751236E-02 -8.00000 0.751236E-02 -5.39697 -0.875163E-01 -8.00000 -0.875163E-01 -.863360 .351629 -8.00000 .351629 6.22887 .267132 -8.00000 .267132 -1.07934 -.438338 -8.00000 -.438338 6.61216 .453794 -8.00000 .453794 7.11977 .597917 -8.00000 .597917 -.110453 .761492 -8.00000 .761492 -4.03012 -.814686 -8.00000 -.814686 TEST04 P00_BOUNDARY_EPS makes an EPS file containing an image of the boundary. The boundary will be drawn to an absolute fineness of H = 0.601347-153 The boundary nodes will NOT be shown. Test number = 1 Title: = "#1: The unit circle." H = .100000 Boundary image file is "p01_boundary.eps". Test number = 2 Title: = "#2: The circle with concentric circular hole." H = .100000 Boundary image file is "p02_boundary.eps". Test number = 3 Title: = "#3: The unit square with circular hole." H = .100000 Boundary image file is "p03_boundary.eps". Test number = 4 Title: = "#4: The unit hexagon with hexagonal hole." H = 0.866025E-01 Boundary image file is "p04_boundary.eps". Test number = 5 Title: = "#5: The horn." H = 0.500000E-01 Boundary image file is "p05_boundary.eps". Test number = 6 Title: = "#6: The superellipse with superelliptical hole." H = .100000 Boundary image file is "p06_boundary.eps". Test number = 7 Title: = "#7: Bicycle seat (implicit)." H = .300000 Boundary image file is "p07_boundary.eps". Test number = 8 Title: = "#8: Pie slice with notch and hole." H = 0.258819E-01 Boundary image file is "p08_boundary.eps". Test number = 9 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." H = 0.500000E-01 Boundary image file is "p09_boundary.eps". Test number = 10 Title: = "#10: The unit square." H = 0.500000E-01 Boundary image file is "p10_boundary.eps". Test number = 11 Title: = "#11: The L-shaped region." H = 0.500000E-01 Boundary image file is "p11_boundary.eps". Test number = 12 Title: = "#12: John Shadid's H-shaped region." H = 0.500000E-01 Boundary image file is "p12_boundary.eps". Test number = 13 Title: = "#13: The Sandia Fork region." H = 5.00000 Boundary image file is "p13_boundary.eps". Test number = 14 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." H = 26.7000 Boundary image file is "p14_boundary.eps". Test number = 15 Title: = "#15: Sangbum Kim's forward step." H = .100000 Boundary image file is "p15_boundary.eps". TEST05 P00_POINTS_EPS makes an EPS file containing an image of a set of points, and the boundary. Test number = 1 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#1: The unit circle." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p01_points.eps". Test number = 2 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#2: The circle with concentric circular hole." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p02_points.eps". Test number = 3 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#3: The unit square with circular hole." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p03_points.eps". Test number = 4 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#4: The unit hexagon with hexagonal hole." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p04_points.eps". Test number = 5 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#5: The horn." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p05_points.eps". Test number = 6 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.500000E-01 Title: = "#6: The superellipse with superelliptical hole." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p06_points.eps". Test number = 7 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = .392699 Title: = "#7: Bicycle seat (implicit)." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p07_points.eps". Test number = 8 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.250000E-01 Title: = "#8: Pie slice with notch and hole." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p08_points.eps". Test number = 9 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.250000E-01 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p09_points.eps". Test number = 10 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.250000E-01 Title: = "#10: The unit square." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p10_points.eps". Test number = 11 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.250000E-01 Title: = "#11: The L-shaped region." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p11_points.eps". Test number = 12 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 0.250000E-01 Title: = "#12: John Shadid's H-shaped region." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p12_points.eps". Test number = 13 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 2.50000 Title: = "#13: The Sandia Fork region." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p13_points.eps". Test number = 14 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = 16.3500 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p14_points.eps". Test number = 15 Spatial dimension = 2 Relative spacing = 0.250000E-01 Actual spacing H = .400000 Title: = "#15: Sangbum Kim's forward step." DEBUG :about to allocate DEBUG :about to sample DEBUG :about to EPS Boundary image file is "p15_points.eps". TEST06 P00_POLY_WRITE creates a POLY file out of the boundary data for each problem. The number of tests available = 15 TEST07 P00_HEX_GRID returns points inside a region that lie on a hexagonal grid. P00_HEX_GRID_COUNT counts the number of hex grid points inside a region. DTABLE_WRITE writes sets of points to a file. We use a hex grid spacing of H = 0.250000E-01 P00_HEX_GRID_COUNT reports that the number of hex grid points will be 409 A few of the points: Row 1 2 Col .487500 0.216506E-01 .475000 0.138778E-16 .487500 -0.216506E-01 .500000 0.433013E-01 .475000 0.433013E-01 The hex grid points were written to the file "p08_hex_grid_points.txt". TEST_TRIANGULATION_PRB Normal end of execution. May 17 2007 11:02:36.601 AM