6 October 2025 6:51:35.295 PM test_triangulation_test(): Fortran90 version Test test_triangulation(). 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 = 0.00000 0.00000 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 = 0.400000 Center = 0.00000 0.00000 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 = 0.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 = 0.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 = 0.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 = 0.00000 0.00000 and radius R1 = 1.00000 The interior hole is a circle C2 with CENTER2 = 0.600000 0.00000 and radius R2 = 0.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 0.500000 0.500000 and "radius" R1 = 0.500000 Hexagon 1 has "center" at 0.250000 0.750000 and "radius" R2 = 0.100000 Hexagon 2 has "center" at 0.600000 0.400000 and "radius" R3 = 0.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 = 0.500000 0.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_SAMPLE returns sample points from the region. Test number = 1 Title: = "#1: The unit circle." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.404550 -0.912746 -0.371602 0.738094 0.848338 0.467545 0.101632 -0.655444 0.924916 -0.234173 -0.439765 0.439653 -0.151612 -0.731669 -0.278048 0.634168 -0.324782 0.714324 0.525710 -0.439059 0.829354 0.184113 -0.448114 0.791406 -0.350767E-01 -0.640491 -0.307864 -0.897472 -0.736533 0.323853 -0.383603 -0.138981E-01 -0.856869E-01 0.125841 -0.180410 -0.150193 -0.378428 -0.840095 0.470236 0.213622 Test number = 2 Title: = "#2: The circle with concentric circular hole." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.633995 0.462722 -0.642880 0.348630 -0.232771 -0.949864 0.512716 -0.379063 0.637489 0.403928E-01 0.616746 0.694173 0.159769 0.677948 0.683661 -0.413384 -0.782627 0.440477 -0.245847 0.631595 0.717017 -0.651859 -0.971975E-01 0.894590 0.708706 -0.445269 0.700878 -0.395161 -0.843274 0.156925 -0.619114 -0.180304 -0.402763 0.135910 -0.339237 0.332691 0.387266 0.533366 -0.210656 -0.498940 Test number = 3 Title: = "#3: The unit square with circular hole." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col -0.872823 -0.574644 0.348541 -0.437857 0.535556 -0.206121 -0.645943 -0.580975E-01 -0.267318 -0.798499 -0.838875 -0.307535 0.506350 0.409232 0.152691 0.392735 -0.684995 -0.508800 -0.534642 -0.647185E-01 -0.211543E-02 0.474219 0.656188 -0.487502 -0.399286 -0.939454 -0.852583 0.127023 0.121196E-01 -0.821481 0.197914 0.822183 -0.488424 0.837592 0.872106 0.783982 0.798975 -0.873040 -0.370397 0.964862 Test number = 4 Title: = "#4: The unit hexagon with hexagonal hole." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.183639 0.836256 0.458592 0.348049 0.408687 -0.669227 0.512128 -0.447955 0.581459 0.413709 0.455230 0.314333 -0.587725 -0.121786 -0.318587 0.683110 0.588845 0.323945 0.539676 -0.557906 -0.454813 0.711226 0.760691 -0.382586E-01 -0.621804 -0.254844 -0.594312 0.169079 0.552545 0.331636 -0.814697 -0.314735 0.787786 0.298982 0.463363 -0.107048 0.316624 -0.607498 0.670026 0.490466 Test number = 5 Title: = "#5: The horn." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col -0.544032E-01 0.540720 -0.173702 0.835357 0.286966 0.512650 0.966138 0.373025E-01 -0.199231 0.900230 0.902853 0.340847 0.170143 0.819449 -0.512355E-01 0.877552 -0.358137 0.928499 -0.326224 0.781349 0.779901 0.539214 0.616009 0.596225 0.576400 0.336137 -0.397009 0.857271 -0.177381 0.715587 0.166354 0.675957 0.489769 0.774226 0.961583 0.631716E-01 -0.415186 0.686218 0.322711 0.280165 Test number = 6 Title: = "#6: The superellipse with superelliptical hole." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col -0.411693E-01 0.828653 0.224037 0.862461 0.764802 0.270910 -0.620103 -0.953262E-01 -0.644450 0.887728 -0.112031 -0.657742 -0.278254 -0.682929 0.559878 0.185195 -0.738834 0.214323 -0.600309 -0.742112 -0.322173 0.593681 0.875716 0.239936 -0.784075 0.643096 0.269915 -0.959482 -0.376026 0.943585 0.545338 -0.788990E-01 -0.755420 -0.672893 0.837063 0.769893 -0.498352 0.951474 0.801747 0.481958 Test number = 7 Title: = "#7: Bicycle seat (implicit)." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 4.68816 -2.44598 -6.37055 -1.31058 -3.53515 -1.45471 -2.34758 -1.19027 7.67207 0.280265E-01 -0.942268 -1.49734 5.23171 -0.848230 -5.64497 -1.86857 -1.36319 -3.59178 1.97898 -0.930330 -5.44218 -2.28843 4.45600 -0.855701 -5.38927 -2.72665 0.836178 -2.76738 5.85852 0.379480 -4.41022 -3.13271 5.84023 -2.80764 -4.05293 -1.50080 -6.44167 -2.29791 5.51917 -3.77164 Test number = 8 Title: = "#8: Pie slice with notch and hole." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.793609 0.331610E-01 0.499671 0.730366E-01 0.728435 0.386706E-01 0.780302 0.619959E-01 0.350289 -0.360825E-02 0.397508 0.100895 0.784300 0.177750 0.947747 -0.869511E-01 0.355696 -0.729940E-01 0.867298 -0.894365E-01 0.911805 -0.209872 0.726858 -0.140068 0.771558 -0.107978 0.869079 0.186811E-01 0.701368 -0.127068E-01 0.751802 0.156492 0.821985 0.250535E-01 0.709188 0.802867E-01 0.860640 -0.349909E-01 0.430354 -0.108950 Test number = 9 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.817486 0.476651 0.587940 0.769496 0.354987 0.108045 0.557924 0.635801 0.765150 0.429813 0.455874 0.200018 0.824302 0.122160 0.519575 0.338589 0.980132 0.595134 0.185542 0.616700 0.178278 0.139624E-01 0.877423 0.899003 0.568370E-01 0.793721 0.299160 0.328984 0.346206 0.142847 0.865167E-01 0.285984 0.400343 0.366817 0.683294 0.745694 0.204736E-03 0.399991 0.287592 0.158518 Test number = 10 Title: = "#10: The unit square." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.156113 0.247219 0.330088 0.358807 0.641408 0.949077E-01 0.174669 0.750680 0.525055E-01 0.651072 0.425654 0.674065E-01 0.988572 0.411953 0.875135 0.523584 0.493558 0.599020 0.586207 0.104815 0.129612 0.803837 0.496856 0.406968 0.112179 0.737859 0.352779 0.586799 0.565188 0.979559 0.356339 0.914803 0.693984 0.271635E-01 0.924399 0.901933 0.948306 0.433501 0.197616 0.344517 Test number = 11 Title: = "#11: The L-shaped region." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.478807 0.644395 0.460062 0.818912 0.865314 0.320853 0.716123 0.369317 0.428619 0.492048 0.168919E-01 0.135158 0.460779 0.852125 0.609317 0.437357 0.841848 0.249553 0.854585 0.212760 0.169016 0.855058 0.702833E-02 0.802718 0.231588 0.251786 0.181322 0.892930 0.333831 0.160494 0.348010 0.613929 0.824331E-01 0.181997 0.175785E-01 0.778847 0.821911E-01 0.481437 0.573528 0.448304E-01 Test number = 12 Title: = "#12: John Shadid's H-shaped region." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 0.284747 0.371938 0.312971E-01 0.865618 0.925279 0.507000 0.766266 0.510488 0.458726 0.536088 0.251869 0.314904 0.408881 0.597305 0.965471 0.725812 0.174591E-01 0.606041 0.761951 0.533318 0.723332 0.666691 0.235922 0.735811 0.604299 0.366000 0.466501E-01 0.339077 0.665577 0.495961 0.667718 0.156613 0.719900 0.411009 0.254153 0.552237E-01 0.641282 0.913885 0.471577 0.212413 Test number = 13 Title: = "#13: The Sandia Fork region." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 18.2109 17.2353 45.3279 32.2049 31.9017 34.4334 49.3826 60.3222 40.7411 36.2163 13.6416 13.8641 30.5459 23.4254 45.8786 85.9318 50.9748 57.6865 45.7842 67.9611 65.6479 35.0577 44.4563 34.9883 52.4928 38.9864 35.4484 34.8967 79.0295 25.3129 69.6236 23.4013 53.9894 74.4174 50.7077 67.4114 15.7016 17.5372 48.0050 82.1831 Test number = 14 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 326.009 338.643 334.950 313.468 617.136 268.320 312.094 278.215 395.059 557.215 446.377 545.674 137.289 299.107 210.009 250.614 177.462 349.499 197.408 396.927 352.079 438.267 368.521 560.249 245.040 363.837 265.982 398.342 200.089 336.404 339.477 283.568 270.287 742.552 292.561 287.491 493.203 484.563 341.694 539.655 Test number = 15 Title: = "#15: Sangbum Kim's forward step." Spatial dimension M = 2 Number of samples N = 20 The sample points: Row 1 2 Col 1.99764 -0.547618 6.75398 0.989584E-01 -1.64825 0.108698 1.26089 0.401966 -5.12819 -0.836942 -1.61794 -0.913194 1.34804 0.595202 -2.68969 -0.500753 1.57115 0.902693E-01 1.62316 0.455248E-01 -2.85706 -0.372415 5.15978 0.548120 7.32671 0.289115 6.60721 0.845695 -1.72340 -0.123806 -1.70470 0.966939 6.21471 0.802017 1.34967 0.926075 2.55965 0.941238 -0.127046 -0.943132 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 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.361901E-01 -0.127292 0.273470 -0.961881 -0.559532 0.644343 -0.655667 0.755050 0.207064 0.263893 0.617305 0.786724 -0.148740 0.623010 -0.232219 0.972664 0.685673 0.602830 0.751019 0.660281 0.465555 -0.867211 0.472992 -0.881066 0.110580 0.509227 0.212206 0.977225 0.521139 0.216973 0.923183 0.384361 -0.663875 -0.252147 -0.934842 -0.355064 0.619841 -0.701064 0.662375 -0.749172 -0.660784 0.285208 -0.918128 0.396284 -0.392412E-02 0.432098E-01 -0.904435E-01 0.995902 -0.307987 0.208870 -0.827627 0.561278 0.508163 -0.167552 0.949708 -0.313138 -0.333565 0.152346 -0.909620 0.415442 0.528865 -0.192093 0.939920 -0.341396 -0.209625 0.825911 -0.246010 0.969267 0.328921 -0.406453 0.629068 -0.777350 -0.896973E-01 0.731201 -0.121759 0.992560 -0.823285 -0.230102 -0.963091 -0.269177 Test number = 2 Title: = "#2: The circle with concentric circular hole." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.288443 -0.705758 0.378323 -0.925674 0.718322 -0.103419 0.989794 -0.142504 -0.449175 -0.371039E-01 -0.398642 -0.329296E-01 -0.928656 0.195543 -0.978542 0.206047 0.458983 0.420299 0.295001 0.270138 0.375649E-01 0.659717 0.227395E-01 0.399353 -0.634303 -0.977791E-01 -0.395330 -0.609410E-01 0.177882 -0.508660 0.132042 -0.377578 -0.400186 -0.328946 -0.309007 -0.253998 -0.762830 0.272267 -0.941809 0.336147 -0.713091 -0.340994 -0.902159 -0.431404 0.856615 0.407501 0.903028 0.429581 0.412454 -0.120334 0.383991 -0.112030 -0.494681E-01 -0.460222 -0.427487E-01 -0.397709 -0.144938 0.618280 -0.912936E-01 0.389443 -0.524666 -0.700523 -0.599469 -0.800398 -0.497895 0.662221 -0.600950 0.799287 0.337954 -0.894547 0.353413 -0.935467 -0.506732E-01 -0.405016 -0.496585E-01 -0.396906 0.148111 -0.949768 0.154082 -0.988058 Test number = 3 Title: = "#3: The unit square with circular hole." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.966016 -0.853212E-01 1.00000 -0.853212E-01 -0.360710 0.595541 -0.207227 0.342136 -0.366810 0.360625 -0.285237 0.280428 0.607697 -0.849274 0.607697 -1.00000 -0.973586 -0.607547E-01 -1.00000 -0.607547E-01 0.667007 0.965900 0.667007 1.00000 0.719533 -0.737005 0.719533 -1.00000 -0.419464 -0.862811 -0.419464 -1.00000 0.964976 0.446808 1.00000 0.446808 -0.343181 -0.586777 -0.201941 -0.345282 -0.726080 -0.479499 -1.00000 -0.479499 -0.748014 -0.131216E-01 -1.00000 -0.131216E-01 -0.699996 0.375714 -1.00000 0.375714 -0.680953 0.185670 -0.385912 0.105224 0.754059 -0.530677 1.00000 -0.530677 0.861599 -0.304424 1.00000 -0.304424 0.565594 0.633714 0.565594 1.00000 0.860486 -0.376247 1.00000 -0.376247 0.926689 0.582147E-01 1.00000 0.582147E-01 0.222093 0.863548 0.222093 1.00000 Test number = 4 Title: = "#4: The unit hexagon with hexagonal hole." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) -0.449243 -0.685049 -0.565676 -0.752272 0.760247 -0.506918E-02 0.937867 -0.107618 0.862825 -0.341508E-01 0.950919 -0.850116E-01 -0.926852 -0.759241E-01 -0.948837 -0.886169E-01 0.657281 -0.146967 0.850682 -0.258627 -0.807377E-01 0.516302 -0.534943E-01 0.469115 -0.630954 0.553851 -0.667914 0.575190 -0.478840 0.531205 -0.639692 0.624073 -0.699642 -0.173072 -0.849968 -0.259863 0.261995 -0.349505 0.261663 -0.348929 0.308758 -0.612607 0.182809 -0.394455 0.734311 -0.451383 0.738123 -0.453584 0.797720 0.302061 0.818634 0.314136 -0.945570E-01 0.550736 -0.489484E-01 0.471740 0.487517 0.734563 0.553804 0.772834 -0.240150 -0.603635 -0.135237 -0.421921 0.186842E-01 0.727489 0.186842E-01 0.866025 -0.368736 0.701175 -0.368736 0.866025 0.295463E-01 0.726433 0.295463E-01 0.866025 0.248563 -0.374998 0.240549 -0.361119 Test number = 5 Title: = "#5: The horn." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.181613 0.209922 0.181613 0.00000 0.905717 0.924104E-01 0.994835 0.101503 0.191763 0.780208 0.204309 0.796750 0.586099 0.526898 0.481989 0.471270 0.187784E-01 0.982495 0.191095E-01 0.999817 0.470205 0.571911 0.435678 0.549220 0.747468 0.448468 0.857499 0.514485 0.660066 0.354187 0.548460 0.316897 -0.159412 0.945851 -0.153488 0.969140 -0.168312 0.572659 -0.249501E-01 0.927005 0.207252 0.708850E-01 0.207252 0.00000 0.765446 0.243007 0.578946 0.204119 0.962983E-01 0.282627 0.150000 0.00000 0.765374 0.470128 0.852091 0.523393 0.157142 0.194432 0.157142 0.00000 0.380679 0.278791 0.541751 0.336312 0.136566 0.241314 0.150000 0.00000 -0.750157E-01 0.537668 0.117282 0.855815 0.234823 0.761601 0.240277 0.768144 -0.462875 0.797119 -0.502162 0.864774 Test number = 6 Title: = "#6: The superellipse with superelliptical hole." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.623313 -0.257946 0.492943 -0.242437 -0.425356 0.470468 -0.400471 0.437926 -0.970762 -0.230517 -0.999289 -0.230869 -0.677159 -0.161519 -0.498810 -0.156057 -0.646846 0.193175 -0.497617 0.185451 0.520281 -0.973245 0.521443 -0.980982 0.192629 0.589235 0.187702 0.497499 -0.734081 0.544282 -0.965944 0.599794 -0.384432 0.495856 -0.366014 0.459430 0.700435 0.664462 0.878287 0.797726 -0.869926 -0.375710 -0.994587 -0.382819 0.335593 0.663774 0.294961 0.484121 -0.436820 0.871088 -0.447813 0.989791 0.436682 -0.615643 0.362098 -0.461366 0.659975 0.895077 0.677606 0.942528 -0.903034 0.804244E-01 -0.999990 0.804753E-01 -0.894891 0.116689 -0.999953 0.116856 -0.711333 0.444098E-01 -0.499992 0.442630E-01 -0.382417 0.832275 -0.392365 0.994021 -0.234449 0.515481 -0.232227 0.494079 Test number = 7 Title: = "#7: Bicycle seat (implicit)." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) -4.23034 -0.489364 -4.21758 -0.474863 6.61195 -2.34578 6.67740 -2.38760 5.86563 -1.26747 6.97517 -1.88950 1.85086 -2.00994 2.53209 -0.819933 -4.30318 -2.96461 -5.10247 -4.10930 1.87302 -0.863412 2.14336 -0.541785 0.596536 -4.27786 0.597345 -4.99983 -0.656704 -1.33318 -1.76979 -0.197686 1.85807 -2.52251 2.64350 -0.878496 -6.79590 -0.174788 -7.36493 0.469793 -1.42803 -1.77401 -2.28004 -0.651257 5.76750 -1.97199 6.59153 -2.51940 3.12636 -2.99273 3.54901 -4.79153 5.12474 -0.155172 4.84236 0.129608 -1.44276 -3.41767 -1.46902 -4.99388 4.42465 -1.16608 4.01920 -0.638991 -6.04928 -2.71711 -6.31115 -2.91530 4.53924 -0.262787 4.49410 -0.216555 -4.96653 -1.13766 -4.29599 -0.404468 -5.90165 -1.64898 -6.79893 -2.19216 Test number = 8 Title: = "#8: Pie slice with notch and hole." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.527175 -0.111651 0.519774 -0.139273 0.347141 -0.804652E-01 0.344003 -0.921753E-01 0.809382 -0.472152E-01 0.900000 0.00000 0.908051 -0.223090 0.902996 -0.241957 0.137012 0.223240E-01 0.133415 0.357485E-01 0.898716 0.148533 0.986616 0.163061 0.856021 -0.118944 0.828414 -0.221973 0.829321 0.184739 0.819952 0.219705 0.388612 0.529024E-01 0.375806 0.100697 0.574539 -0.103990 0.576218 -0.971309E-01 0.934584 0.833658E-01 0.958975 0.589749E-01 0.942246 0.460872E-01 0.944167 0.441665E-01 0.445242 -0.997881E-01 0.440363 -0.117995 0.312021 -0.808416E-01 0.311330 -0.834207E-01 0.848418 0.236321E-01 0.900000 0.00000 0.877323 -0.179702 0.863479 -0.231369 0.993460 0.952180E-01 0.994339 0.943391E-01 0.389062 -0.763250E-01 0.382081 -0.102378 0.302994 -0.556934E-01 0.296621 -0.794794E-01 0.289626 0.130519E-01 0.273488 0.732809E-01 Test number = 9 Title: = "#9: Jeff Borggaard's Box with 2 hexagonal holes." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.304961 0.613550E-01 0.304961 0.00000 0.397302 0.208479 0.550000 0.313397 0.441424 0.661852 0.334687 0.723477 0.663526 0.775547 0.663526 1.00000 0.767840 0.761902E-01 0.767840 0.00000 0.844112 0.255123 1.00000 0.255123 0.165392 0.886310 0.200000 0.836603 0.763091 0.383347 0.700000 0.400000 0.599218 0.826250 0.599218 1.00000 0.476810 0.741609 0.350000 0.750000 0.872974 0.570660 1.00000 0.570660 0.767530 0.991628 0.767530 1.00000 0.532240 0.460265 0.534156 0.459159 0.261827 0.551755 0.261827 0.663397 0.779882 0.922916 0.779882 1.00000 0.857106 0.499581 1.00000 0.499581 0.302525 0.393346 0.500000 0.400000 0.135669 0.689705E-01 0.135669 0.00000 0.608228 0.648234 0.608228 0.486603 0.618240 0.124643 0.618240 0.00000 Test number = 10 Title: = "#10: The unit square." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.680358 0.564282 1.00000 0.564282 0.972613 0.562040E-01 1.00000 0.562040E-01 0.336418E-01 0.953827 0.00000 0.953827 0.800166E-01 0.700779E-01 0.800166E-01 0.00000 0.404890 0.590684 0.00000 0.590684 0.792746 0.517838 1.00000 0.517838 0.201081 0.235408 0.00000 0.235408 0.289309 0.687959 0.00000 0.687959 0.476282E-01 0.835042 0.00000 0.835042 0.720694 0.715387 1.00000 0.715387 0.489079E-01 0.914634 0.00000 0.914634 0.255240 0.377538 0.00000 0.377538 0.529803 0.935095 0.532240 0.460265 0.472113 0.278894 0.472113 0.00000 0.312889E-01 0.867579 0.00000 0.867579 0.794389E-01 0.457025 0.00000 0.457025 0.539844 0.524549 1.00000 0.524549 0.954127 0.259510 1.00000 0.259510 0.839589 0.214520 1.00000 0.214520 0.673801 0.520771E-02 0.673801 0.00000 Test number = 11 Title: = "#11: The L-shaped region." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.555880 0.385027 0.555880 0.500000 0.118533 0.156326 0.00000 0.156326 0.689691 0.253056E-01 0.689691 0.00000 0.113696 0.804324 0.00000 0.804324 0.661450 0.146903 0.661450 0.00000 0.752837 0.643549E-01 0.752837 0.00000 0.348282 0.813840 0.500000 0.813840 0.853947 0.155182 1.00000 0.155182 0.142070 0.312832 0.00000 0.312832 0.328999E-01 0.665503 0.00000 0.665503 0.167324 0.249160 0.00000 0.249160 0.584128 0.314736 0.584128 0.500000 0.377413 0.190108 0.377413 0.00000 0.179684 0.352066 0.00000 0.352066 0.146259 0.227437 0.00000 0.227437 0.906315 0.486059E-01 0.906315 0.00000 0.857451E-03 0.563925 0.00000 0.563925 0.889912 0.515964E-01 0.889912 0.00000 0.178386 0.282079 0.00000 0.282079 0.383013 0.309008 0.500000 0.500000 Test number = 12 Title: = "#12: John Shadid's H-shaped region." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 0.946237 0.272947 1.00000 0.272947 0.958082 0.957112 1.00000 0.957112 0.634506 0.929391 0.625000 0.929391 0.260098 0.438891E-02 0.260098 0.00000 0.985605 0.577451 1.00000 0.577451 0.704107 0.555961 0.625000 0.555961 0.873159 0.658663 1.00000 0.658663 0.111670 0.737150 0.00000 0.737150 0.316206 0.405290 0.500000 0.405290 0.763999 0.916562 0.763999 1.00000 0.733765 0.857437 0.625000 0.857437 0.540722E-01 0.720692 0.00000 0.720692 0.145935 0.863822 0.145935 1.00000 0.381468 0.202692 0.500000 0.202692 0.173372 0.141703 0.173372 0.00000 0.813941 0.385353 1.00000 0.385353 0.695819 0.406699 0.625000 0.406699 0.416929 0.476582E-02 0.416929 0.00000 0.311693 0.424609 0.500000 0.424609 0.656080 0.521457 0.625000 0.521457 Test number = 13 Title: = "#13: The Sandia Fork region." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 50.1589 38.6880 55.0000 39.6863 59.7022 37.0621 60.1299 38.6960 38.0912 35.5266 37.2869 37.9259 54.4595 88.5955 55.0000 88.5955 29.9119 33.4859 29.4228 34.3013 24.8207 16.4339 24.8774 16.3969 88.0367 5.94378 89.5204 6.17562 49.7354 36.6456 45.0000 39.6863 21.8478 22.2877 18.6384 24.8284 59.6167 34.8907 60.6286 38.5621 47.1357 50.2113 45.0000 50.2113 50.1067 33.6843 50.0951 29.9998 24.7916 21.8059 27.3109 19.6266 43.3896 34.6305 42.5000 39.2906 48.5075 75.1374 45.0000 75.1374 48.3390 64.4129 45.0000 64.4129 19.0345 13.5043 22.5012 11.9924 42.3059 29.0837 42.3274 29.0023 67.1080 32.0737 68.8253 35.2932 26.0199 27.1789 23.5360 29.9942 Test number = 14 Title: = "#14: Marcus Garvie's Lake Alpha with Beta Island." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 215.283 238.903 205.807 231.110 528.046 430.443 561.783 422.055 480.205 352.067 458.340 297.551 135.327 282.804 133.363 280.914 299.081 203.870 324.341 181.912 323.665 208.640 337.241 197.644 252.507 634.233 245.218 626.617 145.602 386.765 141.662 389.500 493.813 314.373 479.664 288.992 571.098 265.101 567.314 257.780 234.734 196.966 232.719 193.945 295.531 618.932 298.533 646.862 306.317 747.251 307.870 735.488 275.686 256.592 272.940 258.615 321.522 219.233 341.573 202.993 330.248 298.427 360.056 293.979 324.593 328.867 361.508 303.713 367.301 505.825 344.479 502.205 354.106 575.175 357.368 584.810 282.413 599.318 271.153 592.873 Test number = 15 Title: = "#15: Sangbum Kim's forward step." Spatial dimension M = 2 Number of samples N = 20 Pairs of Point and Nearest Boundary Point: (This data can be passed to the PLOT_POINTS routine and plotted with the DASH option.) 7.50993 0.660174 -8.00000 0.660174 0.862794 0.154004 -8.00000 0.154004 -1.22822 -0.204267E-01 -8.00000 -0.204267E-01 -2.49678 -0.822629 -8.00000 -0.822629 5.29146 0.796086 -8.00000 0.796086 -2.28003 -0.330947 -8.00000 -0.330947 5.35134 0.575449 -8.00000 0.575449 1.45098 0.909401 -8.00000 0.909401 1.14932 -0.554575 -8.00000 -0.554575 2.09638 0.404330 -8.00000 0.404330 1.62788 0.402684 -8.00000 0.402684 -4.39709 -0.246511 -8.00000 -0.246511 0.498575 0.788815 -8.00000 0.788815 -7.84407 -0.162390 -8.00000 -0.162390 0.124209E-01 -0.783864 -8.00000 -0.783864 -3.89526 -0.998846 -8.00000 -0.998846 -2.87376 -0.888921 -8.00000 -0.888921 6.64043 0.972176 -8.00000 0.972176 7.86539 0.189070 -8.00000 0.189070 6.06633 0.619690 -8.00000 0.619690 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 = NaN The boundary nodes will NOT be shown. Test number = 1 Title: = "#1: The unit circle." H = 0.100000 Boundary image file is "p01_boundary.eps". Test number = 2 Title: = "#2: The circle with concentric circular hole." H = 0.100000 Boundary image file is "p02_boundary.eps". Test number = 3 Title: = "#3: The unit square with circular hole." H = 0.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 = 0.100000 Boundary image file is "p06_boundary.eps". Test number = 7 Title: = "#7: Bicycle seat (implicit)." H = 0.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 = 0.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 = 0.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 = 0.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 0.487500 0.216506E-01 0.475000 0.138778E-16 0.487500 -0.216506E-01 0.500000 0.433013E-01 0.475000 0.433013E-01 The hex grid points were written to the file "p08_hex_grid_points.txt". test_triangulation_test(): Normal end of execution. 6 October 2025 6:51:35.311 PM