31 August 2011 9:50:04.105 AM GEOMPACK2_PRB: FORTRAN90 version Test the GEOMPACK2 library. TEST01 AREAPG computes the area of a polygon; AREATR computes the area of a triangle; ANGLE computes the polygonal angle at any vertex. Number of polygonal vertices is N = 9 I, XC(I), YC(I) 1 0.5000000D+01 0.0000000D+00 2 0.5000000D+01 0.2000000D+01 3 0.7000000D+01 0.2000000D+01 4 0.9000000D+01 0.4000000D+01 5 0.5000000D+01 0.8000000D+01 6 0.1000000D+01 0.7000000D+01 7 0.2000000D+01 0.5000000D+01 8 0.0000000D+00 0.3000000D+01 9 0.3000000D+01 0.3000000D+01 Area computed directly by AREAPG = 68.0000 I, AREATR(I), ANGLE(I) 1 0.0000000D+00 0.5880026D+00 2 -0.4000000D+01 0.4712389D+01 3 0.0000000D+00 0.2356194D+01 4 0.3200000D+02 0.1570796D+01 5 0.3200000D+02 0.2111216D+01 6 0.1000000D+01 0.1352127D+01 7 0.1600000D+02 0.4390638D+01 8 -0.9000000D+01 0.7853982D+00 9 -0.0000000D+00 0.4124386D+01 Area computed indirectly by summing AREATR = 68.0000 TEST02 CMCIRC determines if a point lies in, on or outside a circle given by 3 points. The points defining the circle: X1,Y1 = 5.00000 0.00000 X2,Y2 = 0.00000 5.00000 X3,Y3 = -5.00000 0.00000 The point to be tested: X0,Y0 = 3.00000 3.00000 Test results: The point is inside the circle. TEST03 Test CVDEC. Test FNDSEP. Reading input file: "cmos.in". CMOS REGION TOLIN = 0.0000000D+00 ANGSPC = 0.3000000D+02 ANGTOL = 0.2000000D+02 Calling DSPGDC to set data structures. MSGLVL = 0 NVC = 30 1 0.0000000 -0.1100000 2 1.8000000 -0.1200000 3 2.4000000 -0.3100000 4 2.5500000 0.2500000 5 3.1500000 0.3700000 6 4.0500000 0.4100000 7 7.3500000 0.4000000 8 7.8000000 0.2100000 9 8.2500000 0.1850000 10 11.8500000 0.1400000 11 12.6000000 0.0100000 12 12.7500000 -0.5500000 13 13.2000000 -0.3900000 14 15.0000000 -0.3500000 15 15.0000000 -0.6300000 16 12.6000000 -0.6300000 17 11.8500000 -0.8100000 18 8.2500000 -0.8300000 19 7.8000000 -0.8000000 20 7.3500000 -0.5800000 21 3.1500000 -0.5600000 22 2.5500000 -0.3900000 23 7.3500000 -1.5000000 24 10.0000000 -1.5000000 25 0.0000000 -0.3900000 26 0.0000000 -7.5000000 27 10.0000000 -7.5000000 28 15.0000000 -7.5000000 29 0.0000000 -10.0000000 30 15.0000000 -10.0000000 6 0 1 6 11 28 36 46 1 2 3 4 5 6 50 1 1 1 2 0 1.5652408 2 25 1 3 35 1.5707963 3 22 1 4 11 0.4899573 4 3 1 5 0 2.9583116 5 2 1 1 0 2.8404718 6 12 2 7 18 2.9932508 7 16 2 8 41 0.4899573 8 15 2 9 0 1.5707963 9 14 2 10 0 1.5485778 10 13 2 6 0 2.8221957 11 3 3 12 3 1.7990402 12 22 3 13 34 2.9277323 13 21 3 14 33 2.8702575 14 20 3 15 45 3.5915501 15 19 3 16 44 2.7534415 16 18 3 17 43 3.0694690 17 17 3 18 42 2.9116032 18 16 3 19 6 2.8871803 19 12 3 20 0 1.7990402 20 11 3 21 0 2.0041379 21 10 3 22 0 2.9824639 22 9 3 23 0 3.1845918 23 8 3 24 0 3.4856096 24 7 3 25 0 2.7451075 25 6 3 26 0 3.0941471 26 5 3 27 0 2.9886123 27 4 3 11 0 2.0299054 28 25 4 29 0 1.5707963 29 26 4 30 50 1.5707963 30 27 4 31 38 1.5707963 31 24 4 32 37 1.5707963 32 23 4 33 36 4.7123890 33 20 4 34 13 1.5755582 34 21 4 35 12 3.4129278 35 22 4 28 2 2.8654956 36 20 5 37 32 1.1160771 37 23 5 38 31 1.5707963 38 24 5 39 30 4.7123890 39 27 5 40 49 1.5707963 40 28 5 41 0 1.5707963 41 15 5 42 7 1.5707963 42 16 5 43 17 2.9060477 43 17 5 44 16 3.3715821 44 18 5 45 15 3.2137163 45 19 5 36 14 3.5297438 46 26 6 47 0 1.5707963 47 29 6 48 0 1.5707963 48 30 6 49 0 1.5707963 49 28 6 50 39 1.5707963 50 27 6 46 29 3.1415927 0 0 initds: nvc= 30 npolg= 6 nvert= 50 nhole= 0 nhola= 0 nrfv= 9 angmin= 28.072 Calling SPDEC2 Calling CVDEC2 NVC = 36 31 5.9125797 -7.5000000 32 0.0000000 -2.4713916 33 10.0029395 -0.8202614 34 15.0000000 -2.6846492 35 8.3424692 -1.5000000 36 7.8308007 -1.5000000 14 1 6 14 32 44 46 22 53 34 60 43 67 45 75 1 2 3 4 5 6 3 3 4 4 5 5 5 5 76 1 1 1 2 0 1.5652408 2 25 1 3 35 1.5707963 3 22 1 4 11 0.4899573 4 3 1 5 0 2.9583116 5 2 1 1 0 2.8404718 6 12 2 7 18 2.9932508 7 16 2 8 41 0.4899573 8 15 2 9 0 1.5707963 9 14 2 10 0 1.5485778 10 13 2 6 0 2.8221957 11 3 3 12 3 1.7990402 12 22 3 13 60 2.9277323 13 21 3 14 33 2.8702575 14 20 3 52 23 2.0838163 15 19 8 16 71 2.7534415 16 18 8 53 22 1.5042282 17 17 7 18 42 2.9116032 18 16 7 19 6 2.8871803 19 12 7 20 0 1.7990402 20 11 7 21 0 2.0041379 21 10 7 22 0 2.9824639 22 9 7 54 16 1.5582970 23 8 8 51 14 2.0330797 24 7 3 25 0 2.7451075 25 6 3 26 0 3.0941471 26 5 3 27 0 2.9886123 27 4 3 11 0 2.0299054 28 25 10 59 0 1.5707963 29 26 9 55 56 1.5707963 30 27 4 31 64 1.5707963 31 24 4 70 72 1.5707963 32 23 4 58 55 1.8059347 33 20 9 34 13 1.5755582 34 21 9 61 59 2.5914329 35 22 10 28 2 2.8654956 36 20 14 37 57 1.1160771 37 23 14 73 74 1.5707963 38 24 5 65 62 1.5751208 39 27 12 40 49 1.5707963 40 28 12 66 0 1.5707963 41 15 11 42 7 1.5707963 42 16 11 43 17 2.9060477 43 17 11 68 66 0.7723735 44 18 5 72 69 1.4392044 45 19 13 76 73 1.4602556 46 26 6 47 0 1.5707963 47 29 6 48 0 1.5707963 48 30 6 49 0 1.5707963 49 28 6 50 39 1.5707963 50 27 6 56 58 3.1415927 51 20 8 15 75 1.5077337 52 8 3 24 0 1.4525300 53 9 8 23 0 1.6262948 54 18 7 63 65 1.5652408 55 31 9 57 32 1.8059347 56 31 6 46 29 3.1415927 57 23 9 33 36 2.9064542 58 31 4 30 50 1.3356579 59 32 10 60 34 1.0253984 60 21 10 35 12 0.8214949 61 32 9 29 0 2.1161942 62 33 12 64 38 1.5806763 63 33 7 17 67 3.1415927 64 24 12 39 30 3.1372682 65 33 5 44 54 1.5609164 66 34 12 67 43 2.1076248 67 17 12 62 63 2.5992086 68 34 11 41 0 1.0339678 69 35 13 71 44 1.4336489 70 35 4 74 76 3.1415927 71 18 13 45 15 1.7745119 72 35 5 38 31 1.7079437 73 36 14 75 45 1.5268238 74 36 4 32 37 3.1415927 75 19 14 36 51 2.0694882 76 36 13 69 70 1.6147689 DECOMP: NVC = 36 NPOLG = 14 NVERT = 76 ANGMIN = 28.0725 TEST04 DIAEDG determines which diagonal of a quadrilateral is to be preferred, based on the circumcircle criterion. The points defining the quadrilateral: P0: X0,Y0 = 0.00000 0.00000 P1: X1,Y1 = 5.00000 0.00000 P2: X2,Y2 = 6.00000 1.00000 P3: X3,Y3 = 1.00000 1.00000 DIAEDG results: Use diagonal P1--P3. TEST05 Reading input file: annulus.in Calling DSPGDC: NVC = 24 1 1.0000000 0.0000000 2 0.9238795 0.3826834 3 0.7071068 0.7071068 4 0.3826834 0.9238795 5 0.0000000 1.0000000 6 -0.3826834 0.9238795 7 -0.7071068 0.7071068 8 -0.9238795 0.3826834 9 -1.0000000 0.0000000 10 -0.9238795 -0.3826834 11 -0.7071068 -0.7071068 12 -0.3826834 -0.9238795 13 0.0000000 -1.0000000 14 0.3826834 -0.9238795 15 0.7071068 -0.7071068 16 0.9238795 -0.3826834 17 0.7500000 0.0000000 18 0.6767767 0.1767767 19 0.5000000 0.2500000 20 0.3232233 0.1767767 21 0.2500000 0.0000000 22 0.3232233 -0.1767767 23 0.5000000 -0.2500000 24 0.6767767 -0.1767767 1 0 1 1 26 1 1 1 2 0 1.3744467 2 2 1 3 0 2.7488939 3 3 1 4 0 2.7488932 4 4 1 5 0 2.7488939 5 5 1 6 0 2.7488934 6 6 1 7 0 2.7488939 7 7 1 8 0 2.7488932 8 8 1 9 0 2.7488939 9 9 1 10 0 2.7488934 10 10 1 11 0 2.7488939 11 11 1 12 0 2.7488932 12 12 1 13 0 2.7488939 13 13 1 14 0 2.7488934 14 14 1 15 0 2.7488939 15 15 1 16 0 2.7488932 16 16 1 17 0 2.7488939 17 1 1 18 26 1.3744467 18 17 1 19 0 1.9634954 19 24 1 20 0 3.9269909 20 23 1 21 0 3.9269908 21 22 1 22 0 3.9269909 22 21 1 23 0 3.9269908 23 20 1 24 0 3.9269909 24 19 1 25 0 3.9269908 25 18 1 26 0 3.9269909 26 17 1 1 17 1.9634954 0 0 TEST06 DTRIS2 constructs the Delaunay triangulation of a set of points in the plane. The number of points to triangulate is 24 The coordinates of the points are: 1 1.0000000 0.0000000 2 0.9238795 0.3826834 3 0.7071068 0.7071068 4 0.3826834 0.9238795 5 0.0000000 1.0000000 6 -0.3826834 0.9238795 7 -0.7071068 0.7071068 8 -0.9238795 0.3826834 9 -1.0000000 0.0000000 10 -0.9238795 -0.3826834 11 -0.7071068 -0.7071068 12 -0.3826834 -0.9238795 13 0.0000000 -1.0000000 14 0.3826834 -0.9238795 15 0.7071068 -0.7071068 16 0.9238795 -0.3826834 17 0.7500000 0.0000000 18 0.6767767 0.1767767 19 0.5000000 0.2500000 20 0.3232233 0.1767767 21 0.2500000 0.0000000 22 0.3232233 -0.1767767 23 0.5000000 -0.2500000 24 0.6767767 -0.1767767 The number of triangles = 30 NLO = 0 I, NOD_TRI(1:3,I), TNBR(1:3,I) 1 9 10 21 -10 3 4 2 12 22 11 6 9 -19 3 21 10 11 1 -8 9 4 21 8 9 5 -3 1 5 7 8 21 -13 4 7 6 22 12 13 2 -37 12 7 20 7 21 8 5 16 8 20 6 7 11 -15 7 9 21 11 22 3 2 10 10 23 21 22 14 9 13 11 5 6 20 -25 8 15 12 22 13 14 6 -55 13 13 23 22 14 10 12 18 14 19 21 23 16 10 20 15 5 20 4 11 17 -33 16 19 20 21 17 7 14 17 4 20 19 15 16 25 18 23 14 15 13 -68 22 19 17 23 24 20 23 28 20 17 19 23 24 14 19 21 2 19 18 29 24 30 22 16 23 15 23 18 -83 23 16 24 23 27 19 22 24 18 19 17 21 20 26 25 4 19 3 17 29 -47 26 1 18 17 30 24 28 27 1 24 16 28 23 -92 28 1 17 24 26 19 27 29 3 19 2 25 21 -77 30 2 18 1 21 26 -89 TRIANGULATION_PLOT_EPS has created an EPS file containing an image of the triangulation, in test06_triangulation_plot.eps TEST07 MSGLVL = 0 NPT = 24 BINEXP = 0.500000 The number of points to triangulate is 24 The coordinates of the points are: 1 1.0000000 0.0000000 2 0.9238795 0.3826834 3 0.7071068 0.7071068 4 0.3826834 0.9238795 5 0.0000000 1.0000000 6 -0.3826834 0.9238795 7 -0.7071068 0.7071068 8 -0.9238795 0.3826834 9 -1.0000000 0.0000000 10 -0.9238795 -0.3826834 11 -0.7071068 -0.7071068 12 -0.3826834 -0.9238795 13 0.0000000 -1.0000000 14 0.3826834 -0.9238795 15 0.7071068 -0.7071068 16 0.9238795 -0.3826834 17 0.7500000 0.0000000 18 0.6767767 0.1767767 19 0.5000000 0.2500000 20 0.3232233 0.1767767 21 0.2500000 0.0000000 22 0.3232233 -0.1767767 23 0.5000000 -0.2500000 24 0.6767767 -0.1767767 ALG = 2 NLO = 0 DELAUNAY_PRINT Information defining a Delaunay triangulation. The number of points is 24 Point coordinates (transpose of internal array) 1 2 1 1.00000 0.00000 2 0.923879 0.382683 3 0.707107 0.707107 4 0.382683 0.923879 5 0.00000 1.00000 6 -0.382683 0.923879 7 -0.707107 0.707107 8 -0.923879 0.382683 9 -1.00000 0.00000 10 -0.923879 -0.382683 11 -0.707107 -0.707107 12 -0.382683 -0.923879 13 0.00000 -1.00000 14 0.382683 -0.923879 15 0.707107 -0.707107 16 0.923879 -0.382683 17 0.750000 0.00000 18 0.676777 0.176777 19 0.500000 0.250000 20 0.323223 0.176777 21 0.250000 0.00000 22 0.323223 -0.176777 23 0.500000 -0.250000 24 0.676777 -0.176777 The number of triangles is 30 Sets of three points are used as vertices of the triangles. For each triangle, the points are listed in counterclockwise order. Nodes that make up triangles (transpose of internal array) 1 2 3 1 4 19 3 2 5 6 20 3 6 7 20 4 21 8 9 5 19 4 20 6 12 22 11 7 13 14 22 8 14 15 23 9 22 23 21 10 19 21 17 11 20 21 19 12 24 1 17 13 23 17 21 14 1 2 18 15 2 3 19 16 18 19 17 17 17 1 18 18 4 5 20 19 7 8 21 20 18 2 19 21 9 10 21 22 10 11 21 23 20 7 21 24 12 13 22 25 15 16 23 26 21 11 22 27 16 1 24 28 22 14 23 29 17 23 24 30 23 16 24 On each side of a given triangle, there is either another triangle, or a piece of the convex hull. For each triangle, we list the indices of the three neighbors, or (if negative) the codes of the segments of the convex hull. Indices of neighboring triangles (transpose of internal array) 1 2 3 1 5 15 -54 2 -9 3 18 3 -57 23 2 4 19 -63 21 5 1 18 11 6 24 26 -72 7 -24 28 24 8 -75 25 28 9 28 13 26 10 11 13 16 11 23 10 5 12 27 17 29 13 29 10 9 14 -45 20 17 15 -5 1 20 16 20 10 17 17 12 14 16 18 -6 2 5 19 -13 4 23 20 14 15 16 21 -66 22 4 22 -20 26 21 23 3 19 11 24 -21 7 6 25 -81 30 8 26 22 6 9 27 -42 12 30 28 7 8 9 29 13 30 12 30 25 27 29 The number of boundary points (and segments) is 16 The segments that make up the convex hull can be determined from the negative entries of the triangle neighbor list. # Tri Side N1 N2 1 18 1 4 5 2 3 1 6 7 3 19 1 7 8 4 21 1 9 10 5 24 1 12 13 6 8 1 14 15 7 25 1 15 16 8 15 1 2 3 9 1 3 3 4 10 2 1 5 6 11 4 2 8 9 12 22 1 10 11 13 6 3 11 12 14 7 1 13 14 15 27 1 16 1 16 14 1 1 2 ALG = 3 NLO = 0 DELAUNAY_PRINT Information defining a Delaunay triangulation. The number of points is 27 Point coordinates (transpose of internal array) 1 2 1 1.00000 0.00000 2 0.923879 0.382683 3 0.707107 0.707107 4 0.382683 0.923879 5 0.00000 1.00000 6 -0.382683 0.923879 7 -0.707107 0.707107 8 -0.923879 0.382683 9 -1.00000 0.00000 10 -0.923879 -0.382683 11 -0.707107 -0.707107 12 -0.382683 -0.923879 13 0.00000 -1.00000 14 0.382683 -0.923879 15 0.707107 -0.707107 16 0.923879 -0.382683 17 0.750000 0.00000 18 0.676777 0.176777 19 0.500000 0.250000 20 0.323223 0.176777 21 0.250000 0.00000 22 0.323223 -0.176777 23 0.500000 -0.250000 24 0.676777 -0.176777 25 -1000.00 -1000.00 26 1000.00 -1000.00 27 0.00000 1000.00 The number of triangles is 49 Sets of three points are used as vertices of the triangles. For each triangle, the points are listed in counterclockwise order. Nodes that make up triangles (transpose of internal array) 1 2 3 1 23 21 22 2 26 27 2 3 27 25 8 4 3 4 19 5 1 26 2 6 4 5 20 7 2 27 3 8 18 2 19 9 3 27 4 10 9 21 8 11 4 27 5 12 10 11 21 13 5 27 6 14 25 26 13 15 6 27 7 16 12 22 11 17 7 27 8 18 13 14 22 19 8 25 9 20 14 15 23 21 9 25 10 22 26 1 16 23 10 25 11 24 19 21 17 25 11 25 12 26 20 21 19 27 12 25 13 28 24 1 17 29 13 26 14 30 23 17 21 31 14 26 15 32 1 2 18 33 15 26 16 34 2 3 19 35 18 19 17 36 5 6 20 37 17 1 18 38 6 7 20 39 7 8 21 40 19 4 20 41 9 10 21 42 20 7 21 43 12 13 22 44 15 16 23 45 21 11 22 46 16 1 24 47 22 14 23 48 17 23 24 49 23 16 24 On each side of a given triangle, there is either another triangle, or a piece of the convex hull. For each triangle, we list the indices of the three neighbors, or (if negative) the codes of the segments of the convex hull. Indices of neighboring triangles (transpose of internal array) 1 2 3 1 30 45 47 2 -9 7 5 3 -42 19 17 4 9 40 34 5 22 2 32 6 11 36 40 7 2 9 34 8 32 34 35 9 7 11 4 10 41 39 19 11 9 13 6 12 23 45 41 13 11 15 36 14 -6 29 27 15 13 17 38 16 43 45 25 17 15 3 39 18 29 47 43 19 3 21 10 20 31 44 47 21 19 23 41 22 5 46 33 23 21 25 12 24 26 30 35 25 23 27 16 26 42 24 40 27 25 14 43 28 46 37 48 29 14 31 18 30 48 24 1 31 29 33 20 32 5 8 37 33 31 22 44 34 7 4 8 35 8 24 37 36 13 38 6 37 28 32 35 38 15 42 36 39 17 10 42 40 4 6 26 41 21 12 10 42 38 39 26 43 27 18 16 44 33 49 20 45 12 16 1 46 22 28 49 47 18 20 1 48 30 49 28 49 44 46 48 The number of boundary points (and segments) is 3 The segments that make up the convex hull can be determined from the negative entries of the triangle neighbor list. # Tri Side N1 N2 1 3 1 27 25 2 14 1 25 26 3 2 1 26 27 ALG = 4 NLO = 0 DELAUNAY_PRINT Information defining a Delaunay triangulation. The number of points is 24 Point coordinates (transpose of internal array) 1 2 1 1.00000 0.00000 2 0.923879 0.382683 3 0.707107 0.707107 4 0.382683 0.923879 5 0.00000 1.00000 6 -0.382683 0.923879 7 -0.707107 0.707107 8 -0.923879 0.382683 9 -1.00000 0.00000 10 -0.923879 -0.382683 11 -0.707107 -0.707107 12 -0.382683 -0.923879 13 0.00000 -1.00000 14 0.382683 -0.923879 15 0.707107 -0.707107 16 0.923879 -0.382683 17 0.750000 0.00000 18 0.676777 0.176777 19 0.500000 0.250000 20 0.323223 0.176777 21 0.250000 0.00000 22 0.323223 -0.176777 23 0.500000 -0.250000 24 0.676777 -0.176777 The number of triangles is 30 Sets of three points are used as vertices of the triangles. For each triangle, the points are listed in counterclockwise order. Nodes that make up triangles (transpose of internal array) 1 2 3 1 13 14 22 2 14 23 22 3 22 11 12 4 8 9 21 5 8 21 7 6 6 7 20 7 21 10 11 8 5 6 20 9 19 20 21 10 4 19 3 11 21 11 22 12 19 4 20 13 2 19 18 14 17 19 23 15 3 19 2 16 5 20 4 17 1 18 17 18 15 16 23 19 20 7 21 20 23 16 24 21 14 15 23 22 24 17 23 23 12 13 22 24 23 19 21 25 9 10 21 26 22 23 21 27 2 18 1 28 18 19 17 29 24 1 17 30 1 24 16 On each side of a given triangle, there is either another triangle, or a piece of the convex hull. For each triangle, we list the indices of the three neighbors, or (if negative) the codes of the segments of the convex hull. Indices of neighboring triangles (transpose of internal array) 1 2 3 1 -63 2 23 2 21 26 1 3 11 -69 23 4 -75 25 5 5 4 19 -12 6 -17 19 8 7 25 -10 11 8 -18 6 16 9 12 19 24 10 12 15 -50 11 7 3 26 12 10 16 9 13 15 28 27 14 28 24 22 15 10 13 -32 16 8 12 -24 17 27 28 29 18 -92 20 21 19 6 5 9 20 18 30 22 21 -54 18 2 22 29 14 20 23 -3 1 3 24 14 9 26 25 -22 7 4 26 2 24 11 27 13 17 -47 28 13 14 17 29 30 17 22 30 29 20 -83 The number of boundary points (and segments) is 16 The segments that make up the convex hull can be determined from the negative entries of the triangle neighbor list. # Tri Side N1 N2 1 21 1 14 15 2 23 1 12 13 3 25 1 9 10 4 4 1 8 9 5 5 3 7 8 6 3 2 11 12 7 6 1 6 7 8 16 3 4 5 9 10 3 3 4 10 8 1 5 6 11 30 3 16 1 12 18 1 15 16 13 1 1 13 14 14 7 2 10 11 15 15 3 2 3 16 27 3 1 2 TEST08 Reading input file: annulus.in Region: ANNULUS REGION input : tol= 0.2220446D-13 angspc= 30.000 angtol= 20.000 kappa= 0.250 dmin= 0.500 nmin= 10 ntrid= 500 2 24 1 1.0000000 0.0000000 2 0.9238795 0.3826834 3 0.7071068 0.7071068 4 0.3826834 0.9238795 5 0.0000000 1.0000000 6 -0.3826834 0.9238795 7 -0.7071068 0.7071068 8 -0.9238795 0.3826834 9 -1.0000000 0.0000000 10 -0.9238795 -0.3826834 11 -0.7071068 -0.7071068 12 -0.3826834 -0.9238795 13 0.0000000 -1.0000000 14 0.3826834 -0.9238795 15 0.7071068 -0.7071068 16 0.9238795 -0.3826834 17 0.7500000 0.0000000 18 0.6767767 0.1767767 19 0.5000000 0.2500000 20 0.3232233 0.1767767 21 0.2500000 0.0000000 22 0.3232233 -0.1767767 23 0.5000000 -0.2500000 24 0.6767767 -0.1767767 1 0 1 1 26 1 1 1 2 0 1.3744467 2 2 1 3 0 2.7488939 3 3 1 4 0 2.7488932 4 4 1 5 0 2.7488939 5 5 1 6 0 2.7488934 6 6 1 7 0 2.7488939 7 7 1 8 0 2.7488932 8 8 1 9 0 2.7488939 9 9 1 10 0 2.7488934 10 10 1 11 0 2.7488939 11 11 1 12 0 2.7488932 12 12 1 13 0 2.7488939 13 13 1 14 0 2.7488934 14 14 1 15 0 2.7488939 15 15 1 16 0 2.7488932 16 16 1 17 0 2.7488939 17 1 1 18 26 1.3744467 18 17 1 19 0 1.9634954 19 24 1 20 0 3.9269909 20 23 1 21 0 3.9269908 21 22 1 22 0 3.9269909 22 21 1 23 0 3.9269908 23 20 1 24 0 3.9269909 24 19 1 25 0 3.9269908 25 18 1 26 0 3.9269909 26 17 1 1 17 1.9634954 0 0 initds: nvc= 24 npolg= 1 nvert= 26 nhole= 0 nhola= 0 nrfv= 7 angmin= 78.750 TEST08: Before call to EQDIS2: NVC = 26 NPOLG = 8 NVERT = 42 ANGMIN = 1.01446 DEBUG: Call EQDIS2 0 0 0.0000000 0.0000000 0.0000000 0.0000000 26 25 0.9593608 -0.2043071 26 0.9593608 0.2043071 8 19 20 21 22 23 24 25 41 1 1 1 1 1 1 1 1 42 1 1 8 40 0 1.3744467 2 2 7 3 0 2.7488939 3 3 7 38 24 1.0144623 4 4 6 36 23 1.0611681 5 5 5 6 0 2.7488934 6 6 5 7 0 2.7488939 7 7 5 34 22 1.2253404 8 8 4 9 0 2.7488939 9 9 4 10 0 2.7488934 10 10 4 32 21 1.6035142 11 11 3 12 0 2.7488932 12 12 3 13 0 2.7488939 13 13 3 14 0 2.7488934 14 14 3 15 0 2.7488939 15 15 3 30 20 1.7344309 16 16 2 27 0 2.7488939 17 1 1 18 26 1.3744467 18 17 1 19 0 1.9634954 19 24 1 29 27 1.2752147 20 23 2 31 15 1.5380815 21 22 3 33 10 2.5852618 22 21 4 35 7 2.5997874 23 20 5 37 4 2.8283138 24 19 6 39 3 2.3889092 25 18 7 42 40 2.6517762 26 17 8 1 17 1.9634954 27 25 2 28 19 1.4715641 28 24 2 20 0 2.6517762 29 25 1 17 0 1.6700285 30 23 3 21 0 2.3889092 31 15 2 16 0 1.0144623 32 22 4 22 0 1.3417291 33 10 3 11 0 1.1453797 34 21 5 23 0 1.3272034 35 7 4 8 0 1.5235528 36 20 6 24 0 1.0986771 37 4 5 5 0 1.6877258 38 19 7 25 0 1.5380815 39 3 6 4 0 1.7344309 40 26 8 41 25 1.6700285 41 18 8 26 0 1.2752147 42 26 7 2 0 1.4715641 1 0.5152347E-01 0.5643481E+00 0.8418921E-01 15 2 0.1360774E+00 0.5635911E+00 0.8424574E-01 38 3 0.8977994E+00 0.3185055E+00 0.1120654E+00 143 4 0.8280168E+00 0.2794903E+00 0.1196318E+00 116 5 0.6232188E+00 0.3176521E+00 0.1122158E+00 99 6 0.1604539E+00 0.4556956E+00 0.9368986E-01 37 7 0.1360774E+00 0.5635911E+00 0.8424574E-01 38 8 0.5152347E-01 0.5643481E+00 0.8418921E-01 15 TEST08: After call to EQDIS2: NVC = 26 NPOLG = 8 NVERT = 42 NTRIE = 501 ANGMIN = 58.1244 TEST09 LRLINE determines if a point is to the right, left, or on a directed line that is a directed distance away from a directed line from one point to another. The directed base line goes from 0.00000 -1.00000 to 1.00000 0.00000 The directed line distance is 0.00000 The point to be located is 0.00000 0.00000 The point is to the right of the line. TEST10 LUFAC factors a linear system; LUSOL solves a factored linear system. a, b 0.5198728 0.1653477 -0.9338796 0.3181258 0.0694666 0.0109244 0.7966077 0.7332229 0.5742402 2.1149952 -0.0586178 -0.5717121 -0.7495796 0.6238850 -0.7560245 0.5198465 -0.2602758 -0.8046908 0.7928658 0.2477457 ipvt, lu 1 0.5198728 0.1653477 -0.9338796 0.3181258 2 0.0210136 0.7931332 0.7528470 0.5675552 4 -0.1127541 -0.6973210 0.5331380 0.7793204 4 0.9999493 -0.5366249 -0.6187935 1.5377616 x = 1.0000000 1.0000000 1.0000000 1.0000000 emax,esum = 0.8881784E-15 0.1776357E-14 TEST11 Reading input file: annulus.in 8 1 1 0.0515235 0.0417414 0.0330770 3 2 5 0.1360774 0.1816106 0.0330770 4 3 10 0.8977994 0.5769510 0.0804384 8 4 18 0.8280168 1.0106952 0.1414331 6 5 24 0.6232188 0.6755580 0.0793917 5 6 30 0.1604539 0.1160075 0.0366117 1 7 34 0.1360774 0.1816106 0.0330770 5 8 39 0.0515235 0.0417414 0.0330770 2 42 26 1 24 0.0417414 0.0417414 2 25 0.0417414 0.0330770 3 1 0.0417414 0.1160075 4 17 0.0366117 0.1160075 5 23 0.1816106 0.1522409 6 15 0.1522409 0.1522409 7 16 0.0330770 0.1522409 8 25 0.0417414 0.1522409 9 24 0.0366117 0.1522409 10 22 0.5769510 0.1522409 11 10 0.1522409 0.1522409 12 11 0.1522409 0.1522409 13 12 0.1522409 0.1522409 14 13 0.1522409 0.1522409 15 14 0.1522409 0.1522409 16 15 0.1816106 0.0330770 17 23 0.0366117 0.0366117 18 21 0.6755580 0.0366117 19 7 0.1522409 0.0366117 20 8 0.1522409 0.0366117 21 9 0.1522409 0.0366117 22 10 0.5769510 0.0366117 23 22 0.0366117 0.0366117 24 20 0.1160075 0.0366117 25 4 0.1522409 0.0330770 26 5 0.1522409 0.0330770 27 6 0.1522409 28 7 0.6755580 29 21 0.0366117 30 19 0.1160075 31 3 0.1160075 32 4 0.1160075 33 20 0.0366117 34 18 0.0417414 35 26 0.0330770 36 2 0.1522409 37 3 0.1160075 38 19 0.0366117 39 18 0.0366117 40 17 0.0417414 41 1 0.0417414 42 26 0.0417414 TEST12 PRIME finds the smallest prime bigger than a given value. I, PRIME(I) 100 113 200 211 300 307 400 401 500 503 TEST13 Reading input file: ptpg.in DIM = 3 A = 1.00000 B = -2.00000 1 8.0000000 4.0000000 0.0000000 2 10.0000000 4.0000000 2.0000000 3 12.0000000 4.0000000 4.0000000 4 14.0000000 2.0000000 10.0000000 5 15.0000000 4.0000000 7.0000000 6 17.0000000 4.0000000 9.0000000 7 19.0000000 8.0000000 3.0000000 8 18.0000000 10.0000000 -2.0000000 9 16.0000000 8.0000000 0.0000000 10 13.0000000 8.0000000 -3.0000000 11 15.0000000 10.0000000 -5.0000000 12 17.0000000 12.0000000 -7.0000000 13 12.0000000 12.0000000 -12.0000000 14 4.0000000 12.0000000 -20.0000000 15 4.0000000 9.0000000 -14.0000000 16 1.0000000 12.0000000 -23.0000000 17 2.0000000 2.0000000 -2.0000000 18 9.0000000 2.0000000 5.0000000 19 9.0000000 1.0000000 7.0000000 20 7.0000000 1.0000000 5.0000000 21 7.0000000 0.0000000 7.0000000 22 9.0000000 0.0000000 9.0000000 23 11.0000000 0.0000000 11.0000000 24 11.0000000 -2.0000000 15.0000000 25 9.0000000 -2.0000000 13.0000000 26 9.0000000 -3.0000000 15.0000000 27 9.0000000 -4.0000000 17.0000000 28 16.0000000 -2.0000000 20.0000000 inout= 1 pt= 6.0000000 4.0000000 -2.0000000 inout= 1 pt= 7.0000000 6.0000000 -5.0000000 inout= 1 pt= 9.0000000 6.0000000 -3.0000000 inout= 1 pt= 11.0000000 6.0000000 -1.0000000 inout= 1 pt= 9.0000000 3.0000000 3.0000000 inout= 1 pt= 12.0000000 -1.0000000 14.0000000 inout= -1 pt= 0.0000000 4.0000000 -8.0000000 inout= -1 pt= 3.0000000 11.0000000 -19.0000000 inout= -1 pt= 4.0000000 -2.0000000 8.0000000 inout= -1 pt= 8.0000000 14.0000000 -20.0000000 inout= -1 pt= 16.0000000 0.0000000 16.0000000 inout= 0 pt= 12.0000000 4.0000000 4.0000000 inout= 0 pt= 18.0000000 6.0000000 6.0000000 inout= 0 pt= 9.0000000 1.5000000 6.0000000 TEST14 ROTIAR "rotates" an array. Using N = 10 SHIFT = 3 1 2 3 4 5 6 7 8 9 10 Shifted array: 4 5 6 7 8 9 10 1 2 3 TEST15 Reading data file: "shr1.in". N = 4 0.0000000 0.0000000 5.0000000 0.0000000 5.0000000 5.0000000 0.0000000 5.0000000 0.0000000 0.0000000 Calling SHRNK2 NSHR = 4 1.0000000 0.5000000 3.6000000 0.5000000 3.6000000 2.6000000 1.0000000 2.6000000 1.0000000 0.5000000 1 3 50.0000000 2 4 25.0000000 1 3 11.1700000 2 4 4.4100000 TEST16 1 0.7599364 0.5054622 2 0.4706911 0.7599232 3 0.5826738 0.8983039 4 0.2141440 0.3698621 5 0.0330602 0.8666114 6 0.1252102 0.0976546 7 0.6590629 0.7871201 8 0.8119425 0.8964329 9 0.2092338 0.2415745 10 0.5248115 0.2677833 5 0.0330602 0.8666114 6 0.1252102 0.0976546 9 0.2092338 0.2415745 4 0.2141440 0.3698621 2 0.4706911 0.7599232 10 0.5248115 0.2677833 3 0.5826738 0.8983039 7 0.6590629 0.7871201 1 0.7599364 0.5054622 8 0.8119425 0.8964329 TEST17 Reading input file: "annulus.in". ANNULUS REGION input : tol= 0.2220446D-13 angspc= 30.000 angtol= 20.000 kappa= 0.250 dmin= 0.500 nmin= 10 ntrid= 500 MSGLVL = 2 24 1 1.0000000 0.0000000 2 0.9238795 0.3826834 3 0.7071068 0.7071068 4 0.3826834 0.9238795 5 0.0000000 1.0000000 6 -0.3826834 0.9238795 7 -0.7071068 0.7071068 8 -0.9238795 0.3826834 9 -1.0000000 0.0000000 10 -0.9238795 -0.3826834 11 -0.7071068 -0.7071068 12 -0.3826834 -0.9238795 13 0.0000000 -1.0000000 14 0.3826834 -0.9238795 15 0.7071068 -0.7071068 16 0.9238795 -0.3826834 17 0.7500000 0.0000000 18 0.6767767 0.1767767 19 0.5000000 0.2500000 20 0.3232233 0.1767767 21 0.2500000 0.0000000 22 0.3232233 -0.1767767 23 0.5000000 -0.2500000 24 0.6767767 -0.1767767 1 0 1 1 26 1 1 1 2 0 1.3744467 2 2 1 3 0 2.7488939 3 3 1 4 0 2.7488932 4 4 1 5 0 2.7488939 5 5 1 6 0 2.7488934 6 6 1 7 0 2.7488939 7 7 1 8 0 2.7488932 8 8 1 9 0 2.7488939 9 9 1 10 0 2.7488934 10 10 1 11 0 2.7488939 11 11 1 12 0 2.7488932 12 12 1 13 0 2.7488939 13 13 1 14 0 2.7488934 14 14 1 15 0 2.7488939 15 15 1 16 0 2.7488932 16 16 1 17 0 2.7488939 17 1 1 18 26 1.3744467 18 17 1 19 0 1.9634954 19 24 1 20 0 3.9269909 20 23 1 21 0 3.9269908 21 22 1 22 0 3.9269909 22 21 1 23 0 3.9269908 23 20 1 24 0 3.9269909 24 19 1 25 0 3.9269908 25 18 1 26 0 3.9269909 26 17 1 1 17 1.9634954 0 0 TEST17: INITDS: NVC = 24 NPOLG = 1 NVERT = 26 NHOLE = 0 NHOLA = 0 NRFV = 7 ANGMIN = 78.7500 decomp: nvc= 26 npolg= 8 nvert= 42 angmin= 58.124 0 0 0.0000000 0.0000000 0.0000000 0.0000000 26 25 0.9593608 -0.2043071 26 0.9593608 0.2043071 8 19 20 21 22 23 24 25 41 1 1 1 1 1 1 1 1 eqdist: nvc= 26 npolg= 8 nvert= 42 angmin= 58.124 Calling TRIPR2 Calling TRINBR 42 1 1 8 40 0 1.3744467 120 2 2 2 7 3 0 2.7488939 114 4 3 3 7 38 24 1.0144623 102 -5 4 4 6 36 23 1.0611681 86 -6 5 5 5 6 0 2.7488934 95 3 6 6 5 7 0 2.7488939 98 3 7 7 5 34 22 1.2253404 70 -9 8 8 4 9 0 2.7488939 81 2 9 9 4 10 0 2.7488934 83 2 10 10 4 32 21 1.6035142 44 -10 11 11 3 12 0 2.7488932 57 3 12 12 3 13 0 2.7488939 60 3 13 13 3 14 0 2.7488934 63 3 14 14 3 15 0 2.7488939 66 3 15 15 3 30 20 1.7344309 34 -4 16 16 2 27 0 2.7488939 42 1 17 1 1 18 26 1.3744467 31 2 18 17 1 19 0 1.9634954 33 1 19 24 1 29 27 1.2752147 27 2 20 23 2 31 15 1.5380815 34 4 21 22 3 33 10 2.5852618 44 10 22 21 4 35 7 2.5997874 70 9 23 20 5 37 4 2.8283138 86 6 24 19 6 39 3 2.3889092 102 5 25 18 7 42 40 2.6517762 111 2 26 17 8 1 17 1.9634954 31 -2 27 25 2 28 19 1.4715641 27 -2 28 24 2 20 0 2.6517762 43 1 29 25 1 17 0 1.6700285 29 2 30 23 3 21 0 2.3889092 69 1 31 15 2 16 0 1.0144623 38 4 32 22 4 22 0 1.3417291 85 1 33 10 3 11 0 1.1453797 54 3 34 21 5 23 0 1.3272034 101 1 35 7 4 8 0 1.5235528 79 2 36 20 6 24 0 1.0986771 110 1 37 4 5 5 0 1.6877258 92 3 38 19 7 25 0 1.5380815 118 1 39 3 6 4 0 1.7344309 107 3 40 26 8 41 25 1.6700285 111 -2 41 18 8 26 0 1.2752147 119 1 42 26 7 2 0 1.4715641 113 1 289 27 0.7709714 -0.1859535 28 0.8651661 -0.1951303 29 0.9729072 -0.1362048 30 0.9864536 -0.0681024 31 0.9166667 0.0000000 32 0.8333333 0.0000000 33 0.7133884 -0.0883884 34 0.5414214 -0.3414214 35 0.5828427 -0.4328427 36 0.6242641 -0.5242641 37 0.6656854 -0.6156854 38 0.7504613 -0.6422221 39 0.7938159 -0.5773374 40 0.8371704 -0.5124528 41 0.8805250 -0.4475681 42 0.9416201 -0.2934953 43 0.5883884 -0.2133884 44 0.2098503 -0.1954955 45 0.0964773 -0.2142143 46 -0.0168956 -0.2329331 47 -0.1302686 -0.2516519 48 -0.2436416 -0.2703707 49 -0.3570146 -0.2890894 50 -0.4703876 -0.3078082 51 -0.5837606 -0.3265270 52 -0.6971335 -0.3452458 53 -0.8105065 -0.3639646 54 -0.8696863 -0.4637893 55 -0.8154931 -0.5448951 56 -0.7613000 -0.6260010 57 -0.6260010 -0.7613000 58 -0.5448951 -0.8154931 59 -0.4637892 -0.8696863 60 -0.2870126 -0.9429096 61 -0.1913417 -0.9619397 62 -0.0956709 -0.9809699 63 0.0956709 -0.9809699 64 0.1913417 -0.9619397 65 0.2870126 -0.9429096 66 0.4637893 -0.8696863 67 0.5448951 -0.8154931 68 0.6260010 -0.7613000 69 0.4116116 -0.2133884 70 0.1542893 0.0707107 71 0.0585786 0.1414214 72 -0.0371320 0.2121320 73 -0.1328427 0.2828427 74 -0.2285534 0.3535534 75 -0.3242641 0.4242641 76 -0.4199748 0.4949748 77 -0.5156854 0.5656854 78 -0.6113961 0.6363961 79 -0.7793644 0.5989657 80 -0.8516219 0.4908245 81 -0.9492530 0.2551223 82 -0.9746265 0.1275611 83 -0.9746265 -0.1275611 84 -0.9492530 -0.2551223 85 0.2866116 -0.0883884 86 0.3317176 0.2835057 87 0.3402119 0.3902346 88 0.3487062 0.4969636 89 0.3572005 0.6036926 90 0.3656948 0.7104216 91 0.3741891 0.8171505 92 0.2870126 0.9429096 93 0.1913417 0.9619397 94 0.0956709 0.9809699 95 -0.0956709 0.9809699 96 -0.1913417 0.9619397 97 -0.2870126 0.9429096 98 -0.4637893 0.8696863 99 -0.5448951 0.8154931 100 -0.6260010 0.7613000 101 0.2866116 0.0883884 102 0.5345178 0.3261845 103 0.5690356 0.4023689 104 0.6035534 0.4785534 105 0.6380712 0.5547379 106 0.6725890 0.6309223 107 0.6260010 0.7613000 108 0.5448951 0.8154931 109 0.4637892 0.8696863 110 0.4116116 0.2133884 111 0.7709714 0.1859535 112 0.8651661 0.1951303 113 0.9416201 0.2934953 114 0.8805250 0.4475681 115 0.8371704 0.5124528 116 0.7938159 0.5773374 117 0.7504613 0.6422221 118 0.5883884 0.2133884 119 0.7133884 0.0883884 120 0.9864536 0.0681024 121 0.9729072 0.1362048 122 0.9241641 -0.0636696 123 0.8469267 -0.0972690 124 0.7663617 -0.1225261 125 0.8647757 -0.2634361 126 0.6872287 -0.2658199 127 0.7619184 -0.3047915 128 0.8366081 -0.3437632 129 0.5948904 -0.3126640 130 0.6695802 -0.3516356 131 0.7442699 -0.3906073 132 0.8189596 -0.4295789 133 0.6742322 -0.4490873 134 0.7489219 -0.4880589 135 0.7162290 -0.5660248 136 -0.7562841 -0.4595738 137 -0.6619026 -0.5729536 138 -0.6431015 -0.4624766 139 -0.5630494 -0.6600576 140 -0.5442484 -0.5495805 141 -0.5254473 -0.4391035 142 -0.4629334 -0.7397408 143 -0.4441323 -0.6292638 144 -0.4253312 -0.5187868 145 -0.4065302 -0.4083097 146 -0.3628174 -0.8194241 147 -0.3440163 -0.7089470 148 -0.3252152 -0.5984700 149 -0.3064141 -0.4879930 150 -0.2876130 -0.3775160 151 -0.2499178 -0.8239901 152 -0.2311167 -0.7135131 153 -0.2123156 -0.6030360 154 -0.1935146 -0.4925590 155 -0.1747135 -0.3820820 156 -0.1458595 -0.8805085 157 -0.1270585 -0.7700315 158 -0.1082574 -0.6595544 159 -0.0894563 -0.5490774 160 -0.0706552 -0.4386004 161 -0.0518541 -0.3281234 162 -0.0324007 -0.8817884 163 -0.0135997 -0.7713114 164 0.0052014 -0.6608343 165 0.0240025 -0.5503573 166 0.0428036 -0.4398803 167 0.0616047 -0.3294033 168 0.0835967 -0.8681510 169 0.1023978 -0.7576740 170 0.1211988 -0.6471970 171 0.1399999 -0.5367199 172 0.1588010 -0.4262429 173 0.1776021 -0.3157659 174 0.2009080 -0.8467933 175 0.2197091 -0.7363162 176 0.2385101 -0.6258392 177 0.2573112 -0.5153622 178 0.2761123 -0.4048852 179 0.2949134 -0.2944081 180 0.3134400 -0.8535189 181 0.3322411 -0.7430419 182 0.3510422 -0.6325648 183 0.3698433 -0.5220878 184 0.3886443 -0.4116108 185 0.4074454 -0.3011338 186 0.4395650 -0.7803708 187 0.4583661 -0.6698938 188 0.4771672 -0.5594168 189 0.4959683 -0.4489398 190 0.5732160 -0.6629993 191 -0.8532903 0.2438889 192 -0.8051954 0.3534272 193 -0.7571004 0.4629656 194 -0.8640779 -0.0782536 195 -0.8159830 0.0312848 196 -0.7678880 0.1408231 197 -0.7197931 0.2503614 198 -0.6716982 0.3598997 199 -0.6236032 0.4694380 200 -0.8060544 -0.2436757 201 -0.7579595 -0.1341374 202 -0.7098645 -0.0245991 203 -0.6617696 0.0849393 204 -0.6136746 0.1944776 205 -0.5655797 0.3040159 206 -0.5174848 0.4135542 207 -0.6622464 -0.2137199 208 -0.6141514 -0.1041816 209 -0.5660565 0.0053568 210 -0.5179616 0.1148951 211 -0.4698666 0.2244334 212 -0.4217717 0.3339717 213 -0.5184383 -0.1837641 214 -0.4703434 -0.0742257 215 -0.4222485 0.0353126 216 -0.3741535 0.1448509 217 -0.3260586 0.2543892 218 -0.3986778 -0.2085774 219 -0.3505828 -0.0990391 220 -0.3024879 0.0104993 221 -0.2543930 0.1200376 222 -0.2062980 0.2295759 223 -0.2548697 -0.1786216 224 -0.2067748 -0.0690832 225 -0.1586799 0.0404551 226 -0.1105849 0.1499934 227 -0.1110617 -0.1486657 228 -0.0629668 -0.0391274 229 -0.0148718 0.0704109 230 0.0327463 -0.1187099 231 0.0808412 -0.0091716 232 0.1750791 -0.0921139 233 -0.5056907 0.7242229 234 -0.4207117 0.6504001 235 -0.3540311 0.7406558 236 -0.2873505 0.8309115 237 -0.3271470 0.5881985 238 -0.2604664 0.6784542 239 -0.1937858 0.7687099 240 -0.1271052 0.8589656 241 -0.2421194 0.5144414 242 -0.1754388 0.6046971 243 -0.1087582 0.6949528 244 -0.0420776 0.7852085 245 0.0246029 0.8754642 246 -0.1695103 0.4238751 247 -0.1028298 0.5141308 248 -0.0361492 0.6043865 249 0.0305314 0.6946422 250 0.0972120 0.7848979 251 0.1638925 0.8751536 252 -0.0635610 0.3784367 253 0.0031196 0.4686924 254 0.0698001 0.5589481 255 0.1364807 0.6492038 256 0.2031613 0.7394595 257 0.2698418 0.8297152 258 0.0094131 0.2883645 259 0.0760937 0.3786202 260 0.1427742 0.4688759 261 0.2094548 0.5591316 262 0.2761354 0.6493873 263 0.1127993 0.2394568 264 0.1794799 0.3297125 265 0.2461604 0.4199682 266 0.2161855 0.1905490 267 0.4640927 0.7740895 268 0.4766008 0.6791079 269 0.5699953 0.6716749 270 0.4523336 0.5870532 271 0.5457282 0.5796201 272 0.4280664 0.4949984 273 0.5214610 0.4875654 274 0.4504965 0.3992272 275 0.4262294 0.3071724 276 0.7162290 0.5660248 277 0.6742322 0.4490873 278 0.7489219 0.4880589 279 0.5948904 0.3126640 280 0.6695802 0.3516356 281 0.7442699 0.3906073 282 0.8189596 0.4295789 283 0.6872287 0.2658199 284 0.7619184 0.3047915 285 0.8366081 0.3437632 286 0.8647757 0.2634361 287 0.7663617 0.1225261 288 0.8469267 0.0972690 289 0.9241641 0.0636696 1 16 53 199 313 410 445 483 497 1 1 31 30 488 2 0 2 31 122 30 3 15 1 3 32 122 31 4 2 487 4 32 123 122 5 14 3 5 17 123 32 6 4 485 6 17 124 123 7 10 5 7 33 124 17 8 6 0 8 24 124 33 9 7 0 9 27 124 24 10 8 28 10 27 123 124 11 6 9 11 28 123 27 12 10 25 12 28 29 123 13 14 11 13 25 29 28 0 12 52 14 29 122 123 15 4 12 15 30 122 29 2 14 0 16 125 127 128 26 20 50 17 126 129 130 30 34 18 18 127 126 130 27 17 19 19 127 130 131 18 22 20 20 128 127 131 16 19 21 21 128 131 132 20 24 47 22 131 130 133 19 35 23 23 131 133 134 22 38 24 24 132 131 134 21 23 45 25 27 125 28 26 52 11 26 27 127 125 27 16 25 27 27 126 127 28 18 26 28 24 126 27 29 27 9 29 43 126 24 30 28 0 30 43 129 126 31 17 29 31 23 129 43 32 30 0 32 34 129 23 33 31 173 33 35 129 34 34 32 170 34 35 130 129 35 17 33 35 35 133 130 36 22 34 36 36 133 35 37 35 169 37 36 135 133 39 38 36 38 133 135 134 37 43 23 39 37 135 36 41 37 166 40 15 38 37 0 41 165 41 38 135 37 42 39 40 42 39 135 38 43 41 0 43 39 134 135 44 38 42 44 40 134 39 45 43 0 45 40 132 134 46 24 44 46 41 132 40 47 45 0 47 41 128 132 48 21 46 48 16 128 41 49 47 0 49 42 128 16 50 48 0 50 42 125 128 51 16 49 51 25 125 42 52 50 0 52 28 125 25 25 51 13 53 136 137 138 134 55 194 54 137 139 140 137 58 55 55 138 137 140 53 54 56 56 138 140 141 55 60 192 57 139 142 143 140 63 58 58 140 139 143 54 57 59 59 140 143 144 58 65 60 60 141 140 144 56 59 61 61 141 144 145 60 67 190 62 142 146 147 142 69 63 63 143 142 147 57 62 64 64 143 147 148 63 71 65 65 144 143 148 59 64 66 66 144 148 149 65 73 67 67 145 144 149 61 66 68 68 145 149 150 67 75 188 69 147 146 151 62 145 70 70 147 151 152 69 78 71 71 148 147 152 64 70 72 72 148 152 153 71 80 73 73 149 148 153 66 72 74 74 149 153 154 73 82 75 75 150 149 154 68 74 76 76 150 154 155 75 84 186 77 151 156 157 147 86 78 78 152 151 157 70 77 79 79 152 157 158 78 88 80 80 153 152 158 72 79 81 81 153 158 159 80 90 82 82 154 153 159 74 81 83 83 154 159 160 82 92 84 84 155 154 160 76 83 85 85 155 160 161 84 94 184 86 157 156 162 77 149 87 87 157 162 163 86 96 88 88 158 157 163 79 87 89 89 158 163 164 88 98 90 90 159 158 164 81 89 91 91 159 164 165 90 100 92 92 160 159 165 83 91 93 93 160 165 166 92 102 94 94 161 160 166 85 93 95 95 161 166 167 94 104 182 96 163 162 168 87 151 97 97 163 168 169 96 106 98 98 164 163 169 89 97 99 99 164 169 170 98 108 100 100 165 164 170 91 99 101 101 165 170 171 100 110 102 102 166 165 171 93 101 103 103 166 171 172 102 112 104 104 167 166 172 95 103 105 105 167 172 173 104 114 180 106 169 168 174 97 154 107 107 169 174 175 106 116 108 108 170 169 175 99 107 109 109 170 175 176 108 118 110 110 171 170 176 101 109 111 111 171 176 177 110 120 112 112 172 171 177 103 111 113 113 172 177 178 112 122 114 114 173 172 178 105 113 115 115 173 178 179 114 124 178 116 175 174 180 107 156 117 117 175 180 181 116 126 118 118 176 175 181 109 117 119 119 176 181 182 118 128 120 120 177 176 182 111 119 121 121 177 182 183 120 130 122 122 178 177 183 113 121 123 123 178 183 184 122 132 124 124 179 178 184 115 123 125 125 179 184 185 124 172 175 126 181 180 186 117 159 127 127 181 186 187 126 161 128 128 182 181 187 119 127 129 129 182 187 188 128 133 130 130 183 182 188 121 129 131 131 183 188 189 130 168 132 132 184 183 189 123 131 171 133 188 187 190 129 162 167 134 55 137 136 135 53 198 135 56 137 55 136 134 0 136 11 137 56 137 135 0 137 11 139 137 138 54 136 138 57 139 11 139 137 0 139 58 139 57 140 138 0 140 58 142 139 141 57 139 141 59 142 58 142 140 0 142 59 146 142 143 62 141 143 12 146 59 144 142 0 144 60 146 12 145 143 0 145 60 151 146 146 69 144 146 61 151 60 147 145 0 147 61 156 151 148 77 146 148 62 156 61 149 147 0 149 62 162 156 150 86 148 150 13 162 62 151 149 0 151 13 168 162 152 96 150 152 63 168 13 153 151 0 153 64 168 63 154 152 0 154 64 174 168 155 106 153 155 65 174 64 156 154 0 156 65 180 174 157 116 155 157 14 180 65 158 156 0 158 66 180 14 159 157 0 159 66 186 180 160 126 158 160 67 186 66 161 159 0 161 67 187 186 162 127 160 162 67 190 187 163 133 161 163 68 190 67 164 162 0 164 15 190 68 165 163 0 165 37 190 15 166 164 40 166 36 190 37 167 165 39 167 36 188 190 168 133 166 168 36 189 188 169 131 167 169 35 189 36 170 168 36 170 34 189 35 171 169 33 171 34 184 189 172 132 170 172 34 185 184 173 125 171 173 23 185 34 174 172 32 174 69 185 23 175 173 0 175 69 179 185 176 125 174 176 22 179 69 177 175 0 177 44 179 22 178 176 285 178 44 173 179 179 115 177 179 45 173 44 180 178 284 180 45 167 173 181 105 179 181 46 167 45 182 180 282 182 46 161 167 183 95 181 183 47 161 46 184 182 280 184 47 155 161 185 85 183 185 48 155 47 186 184 278 186 48 150 155 187 76 185 187 49 150 48 188 186 277 188 49 145 150 189 68 187 189 50 145 49 190 188 275 190 50 141 145 191 61 189 191 51 141 50 192 190 273 192 51 138 141 193 56 191 193 52 138 51 194 192 271 194 52 136 138 195 53 193 195 53 136 52 197 194 269 196 10 54 53 0 197 267 197 54 136 53 198 195 196 198 55 136 54 134 197 0 199 191 196 197 259 209 200 200 192 191 197 311 199 201 201 192 197 198 200 211 202 202 193 192 198 309 201 203 203 193 198 199 202 213 307 204 194 201 202 266 216 205 205 195 194 202 262 204 206 206 195 202 203 205 218 207 207 196 195 203 260 206 208 208 196 203 204 207 220 209 209 197 196 204 199 208 210 210 197 204 205 209 222 211 211 198 197 205 201 210 212 212 198 205 206 211 224 213 213 199 198 206 203 212 303 214 201 200 207 266 270 215 215 201 207 208 214 225 216 216 202 201 208 204 215 217 217 202 208 209 216 227 218 218 203 202 209 206 217 219 219 203 209 210 218 229 220 220 204 203 210 208 219 221 221 204 210 211 220 231 222 222 205 204 211 210 221 223 223 205 211 212 222 233 224 224 206 205 212 212 223 301 225 208 207 213 215 272 226 226 208 213 214 225 234 227 227 209 208 214 217 226 228 228 209 214 215 227 236 229 229 210 209 215 219 228 230 230 210 215 216 229 238 231 231 211 210 216 221 230 232 232 211 216 217 231 240 233 233 212 211 217 223 232 299 234 214 213 218 226 274 235 235 214 218 219 234 242 236 236 215 214 219 228 235 237 237 215 219 220 236 244 238 238 216 215 220 230 237 239 239 216 220 221 238 246 240 240 217 216 221 232 239 241 241 217 221 222 240 248 297 242 219 218 223 235 276 243 243 219 223 224 242 249 244 244 220 219 224 237 243 245 245 220 224 225 244 251 246 246 221 220 225 239 245 247 247 221 225 226 246 253 248 248 222 221 226 241 247 295 249 224 223 227 243 279 250 250 224 227 228 249 254 251 251 225 224 228 245 250 252 252 225 228 229 251 256 253 253 226 225 229 247 252 293 254 228 227 230 250 281 255 255 228 230 231 254 257 256 256 229 228 231 252 255 291 257 231 230 232 255 283 289 258 82 191 81 259 312 0 259 82 196 191 260 199 258 260 82 195 196 261 207 259 261 9 195 82 262 260 0 262 9 194 195 263 205 261 263 83 194 9 264 262 0 264 84 194 83 265 263 0 265 84 200 194 268 266 264 266 194 200 201 265 214 204 267 10 53 84 196 268 0 268 53 200 84 269 265 267 269 52 200 53 270 268 195 270 52 207 200 271 214 269 271 51 207 52 272 270 193 272 51 213 207 273 225 271 273 50 213 51 274 272 191 274 50 218 213 275 234 273 275 49 218 50 276 274 189 276 49 223 218 277 242 275 277 48 223 49 278 276 187 278 47 223 48 279 277 185 279 47 227 223 280 249 278 280 46 227 47 281 279 183 281 46 230 227 282 254 280 282 45 230 46 283 281 181 283 45 232 230 284 257 282 284 44 232 45 286 283 179 285 22 85 44 0 286 177 286 85 232 44 287 284 285 287 21 232 85 288 286 0 288 70 232 21 289 287 377 289 70 231 232 290 257 288 290 71 231 70 291 289 376 291 71 229 231 292 256 290 292 72 229 71 293 291 374 293 72 226 229 294 253 292 294 73 226 72 295 293 372 295 73 222 226 296 248 294 296 74 222 73 297 295 369 297 74 217 222 298 241 296 298 75 217 74 299 297 368 299 75 212 217 300 233 298 300 76 212 75 301 299 365 301 76 206 212 302 224 300 302 77 206 76 303 301 363 303 77 199 206 304 213 302 304 78 199 77 305 303 361 305 78 79 199 306 307 304 306 7 79 78 0 305 359 307 79 193 199 308 203 305 308 80 193 79 309 307 0 309 80 192 193 310 202 308 310 8 192 80 311 309 0 311 8 191 192 312 200 310 312 81 191 8 258 311 0 313 233 234 235 362 314 408 314 235 234 237 313 364 315 315 235 237 238 314 319 316 316 236 235 238 407 315 317 317 236 238 239 316 321 404 318 237 241 242 366 326 319 319 238 237 242 315 318 320 320 238 242 243 319 328 321 321 239 238 243 317 320 322 322 239 243 244 321 330 323 323 240 239 244 404 322 324 324 240 244 245 323 332 400 325 241 246 247 367 334 326 326 242 241 247 318 325 327 327 242 247 248 326 336 328 328 243 242 248 320 327 329 329 243 248 249 328 338 330 330 244 243 249 322 329 331 331 244 249 250 330 340 332 332 245 244 250 324 331 333 333 245 250 251 332 342 397 334 247 246 252 325 370 335 335 247 252 253 334 345 336 336 248 247 253 327 335 337 337 248 253 254 336 347 338 338 249 248 254 329 337 339 339 249 254 255 338 349 340 340 250 249 255 331 339 341 341 250 255 256 340 351 342 342 251 250 256 333 341 343 343 251 256 257 342 390 394 344 252 258 259 371 352 345 345 253 252 259 335 344 346 346 253 259 260 345 354 347 347 254 253 260 337 346 348 348 254 260 261 347 356 349 349 255 254 261 339 348 350 350 255 261 262 349 387 351 351 256 255 262 341 350 389 352 259 258 263 344 373 353 353 259 263 264 352 357 354 354 260 259 264 346 353 355 355 260 264 265 354 382 356 356 261 260 265 348 355 385 357 264 263 266 353 375 381 358 100 233 99 360 409 0 359 7 78 100 306 360 0 360 78 233 100 361 358 359 361 77 233 78 362 360 304 362 77 234 233 363 313 361 363 76 234 77 364 362 302 364 76 237 234 365 314 363 365 75 237 76 366 364 300 366 75 241 237 367 318 365 367 75 246 241 368 325 366 368 74 246 75 369 367 298 369 73 246 74 370 368 296 370 73 252 246 371 334 369 371 73 258 252 372 344 370 372 72 258 73 373 371 294 373 72 263 258 374 352 372 374 71 263 72 375 373 292 375 71 266 263 376 357 374 376 70 266 71 378 375 290 377 21 101 70 0 378 288 378 101 266 70 379 376 377 379 20 266 101 380 378 0 380 86 266 20 381 379 428 381 86 264 266 382 357 380 382 86 265 264 383 355 381 383 87 265 86 384 382 427 384 88 265 87 385 383 424 385 88 261 265 386 356 384 386 89 261 88 387 385 423 387 89 262 261 388 350 386 388 90 262 89 389 387 420 389 90 256 262 390 351 388 390 90 257 256 391 343 389 391 91 257 90 392 390 418 392 4 257 91 393 391 416 393 92 257 4 394 392 0 394 92 251 257 395 343 393 395 93 251 92 396 394 0 396 94 251 93 397 395 0 397 94 245 251 398 333 396 398 5 245 94 399 397 0 399 95 245 5 400 398 0 400 95 240 245 401 324 399 401 96 240 95 403 400 0 402 97 236 96 405 403 0 403 96 236 240 402 404 401 404 240 236 239 403 317 323 405 6 236 97 406 402 0 406 98 236 6 407 405 0 407 98 235 236 408 316 406 408 98 233 235 409 313 407 409 99 233 98 358 408 0 410 268 270 271 421 413 411 411 269 268 271 444 410 438 412 270 272 273 422 414 413 413 271 270 273 410 412 436 414 273 272 274 412 425 434 415 109 267 108 417 443 0 416 4 91 109 392 417 0 417 91 267 109 418 415 416 418 90 267 91 419 417 391 419 90 268 267 420 444 418 420 89 268 90 421 419 388 421 89 270 268 422 410 420 422 89 272 270 423 412 421 423 88 272 89 424 422 386 424 87 272 88 425 423 384 425 87 274 272 426 414 424 426 87 275 274 427 432 425 427 86 275 87 429 426 383 428 20 110 86 0 429 380 429 110 275 86 430 427 428 430 19 275 110 431 429 0 431 102 275 19 432 430 462 432 102 274 275 433 426 431 433 103 274 102 434 432 461 434 103 273 274 435 414 433 435 104 273 103 436 434 458 436 104 271 273 437 413 435 437 105 271 104 438 436 457 438 105 269 271 439 411 437 439 106 269 105 441 438 454 440 3 107 106 0 441 481 441 107 269 106 442 439 440 442 108 269 107 443 441 0 443 108 267 269 415 444 442 444 269 267 268 443 419 411 445 277 280 281 459 450 446 446 278 277 281 456 445 447 447 278 281 282 446 452 477 448 280 279 283 460 464 449 449 280 283 284 448 467 450 450 281 280 284 445 449 451 451 281 284 285 450 453 452 452 282 281 285 447 451 475 453 285 284 286 451 468 472 454 105 276 106 455 482 439 455 105 277 276 457 456 454 456 276 277 278 455 446 479 457 104 277 105 458 455 437 458 103 277 104 459 457 435 459 103 280 277 460 445 458 460 103 279 280 461 448 459 461 102 279 103 462 460 433 462 19 279 102 463 461 431 463 118 279 19 464 462 0 464 118 283 279 465 448 463 465 18 283 118 466 464 0 466 111 283 18 467 465 496 467 111 284 283 468 449 466 468 111 286 284 469 453 467 469 112 286 111 470 468 494 470 26 286 112 471 469 492 471 113 286 26 472 470 0 472 113 285 286 473 453 471 473 2 285 113 474 472 0 474 114 285 2 475 473 0 475 114 282 285 476 452 474 476 115 282 114 477 475 0 477 115 278 282 478 447 476 478 116 278 115 479 477 0 479 116 276 278 480 456 478 480 117 276 116 482 479 0 481 3 106 117 440 482 0 482 106 276 117 454 480 481 483 17 287 119 484 497 0 484 17 288 287 485 495 483 485 32 288 17 486 484 5 486 32 289 288 487 491 485 487 31 289 32 489 486 3 488 1 120 31 0 489 1 489 120 289 31 490 487 488 490 121 289 120 491 489 0 491 121 288 289 493 486 490 492 26 112 121 470 493 0 493 112 288 121 494 491 492 494 111 288 112 495 493 469 495 111 287 288 496 484 494 496 18 287 111 497 495 466 497 119 287 18 483 496 0 triang: nvc= 289 ntri= 497 angmin= 30.099 angmax= 103.877 relative frequency of triangle angles 0.0 0.00000 15.0 0.00000 30.0 0.10731 45.0 0.49497 60.0 0.21395 75.0 0.14554 90.0 0.03823 105.0 0.00000 120.0 0.00000 135.0 0.00000 150.0 0.00000 165.0 0.00000 TEST18 Test being SKIPPED FOR NOW. TEST19 TEST BEING SKIPPED FOR NOW. TEST20 XEDGE determines whether two edges, or an edge and a ray, intersect. (An edge is a finite line segment.) (A ray is a semi-infinite line segment.) Edge 1 is from ( 3.0000000000000000 , 0.0000000000000000 ) to ( 3.0000000000000000 , 2.0000000000000000 ). Edge 2 is from ( 0.0000000000000000 , 0.0000000000000000 ) to ( 6.0000000000000000 , 2.0000000000000000 ). The objects intersect at 3.00000 1.00000 (Expecting the answer (3,1) ). TEST21 XLINE finds the intersection of two lines. Each line is defined as the line a given distance to the left of a line through two points. Line 1 is -6.0000000000000000 units left of the line through ( 0.0000000000000000 , 0.0000000000000000 ) and ( 0.0000000000000000 , 1.0000000000000000 ). Line 2 is 0.0000000000000000 units left of the line through ( 0.0000000000000000 , 0.0000000000000000 ) and ( 3.0000000000000000 , 1.0000000000000000 ). The lines intersect at 6.00000 2.00000 (Expecting the answer (6,2) ). GEOMPACK2_PRB: Normal end of execution. 31 August 2011 9:50:04.132 AM