September 8 2005 8:26:28.041 AM GEOMPACK3_PRB: Test the routines in GEOMPACK3. TEST_ANGLE ANGLE computes the angles of a polygon. I X Y 1 1.00000 0.00000 2 0.500000 0.866025 3 -0.500000 0.866025 4 -1.00000 0.00000 5 -0.500000 -0.866025 6 0.500000 -0.866025 I Angle 1 2.09440 2 2.09440 3 2.09440 4 2.09440 5 2.09440 6 2.09440 TEST_AREA AREAPG computes the area of a polygon. AREATR computes the area of a triangle. I, (X,Y) 1 1.00000 0.00000 2 0.500000 0.866025 3 -0.500000 0.866025 4 -1.00000 0.00000 5 -0.500000 -0.866025 6 0.500000 -0.866025 I, Triangle 1 0.00000 2 0.866025 3 1.73205 4 1.73205 5 0.866025 6 0.00000 Area computation by: AREAPG AREATR 5.19615 5.19615 TEST_DHPSRT DHPSRT sorts multidimensional double precision vectors. Unsorted array produced by RANDPT: 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 11 0.2717472 0.6328177 12 0.5525799 0.3104528 13 0.0130475 0.0333018 14 0.8271604 0.9233513 15 0.6098558 0.9847704 16 0.0812091 0.6316019 17 0.4682860 0.4678636 18 0.4022023 0.0043875 19 0.8238210 0.7603937 20 0.5267480 0.1333443 After sorting by DHPSRT: 13 0.0130475 0.0333018 5 0.0330602 0.8666114 16 0.0812091 0.6316019 6 0.1252102 0.0976546 9 0.2092338 0.2415745 4 0.2141440 0.3698621 11 0.2717472 0.6328177 18 0.4022023 0.0043875 17 0.4682860 0.4678636 2 0.4706911 0.7599232 10 0.5248115 0.2677833 20 0.5267480 0.1333443 12 0.5525799 0.3104528 3 0.5826738 0.8983039 15 0.6098558 0.9847704 7 0.6590629 0.7871201 1 0.7599364 0.5054622 8 0.8119425 0.8964329 19 0.8238210 0.7603937 14 0.8271604 0.9233513 TEST_DTRIS3 : DTRIS3 computes the Delaunay "triangulation" of a set of points in 3D. Node coordinates: Row 1 2 3 Col 1 0.00000 0.00000 -1.00000 2 0.00000 -1.00000 0.00000 3 1.00000 0.00000 0.00000 4 0.00000 1.00000 0.00000 5 -1.00000 0.00000 0.00000 6 0.00000 0.00000 1.00000 Tetrahedron indices: Row 1 2 3 4 Col 1 1 2 5 6 2 1 4 5 6 3 1 2 3 6 4 1 3 4 6 BF_NUM = 9 FC_NUM = 12 FACE_NUM = 12 TETRA_NUM = 4 TETRA_BAD_NUM = 0 HT_NUM = 17 min ( ETA ) = 0.857143 min ( RHO ) = 0.803848 min ( SIGMA ) = 0.621320 vm 6 5 2 1 6 4 3 ht 17 9 5 8 0 0 0 0 3 12 4 11 10 0 0 6 0 0 FC 12 1 1 2 3 4 -1 0 5 2 1 2 4 3 -2 0 0 3 1 3 4 2 5 0 -1 4 2 3 4 1 6 0 -1 5 3 4 5 1 6 2 -1 6 1 4 5 3 -6 0 -1 7 1 3 5 4 -7 0 -1 8 3 4 6 2 5 0 -1 9 2 4 6 3 -3 1 -1 10 2 3 6 4 -8 0 -1 11 4 5 6 3 -4 0 -1 12 3 5 6 4 -9 7 -1 BF 9 1 10 7 2 2 9 6 1 3 11 10 2 4 12 9 6 5 0 7 4 6 11 7 2 7 12 6 1 8 12 9 1 9 11 10 7 TEST_IHPSRT IHPSRT sorts multidimensional integer vectors. Unsorted array produced by RANDPT: 21 0 0 After sorting by IHPSRT: 13 0 0 5 0 17 16 1 12 6 2 1 9 4 4 4 4 7 11 5 12 18 8 0 17 9 9 2 9 15 20 10 2 10 10 5 12 11 6 3 11 17 15 12 19 7 13 15 1 15 10 19 16 15 8 16 17 14 16 18 TEST_LU LUFAC computes the LU factorization of a matrix. LUSOL solves a linear system. Matrix A and right hand side 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 tf, ts = 0.000 0.000 TEST_MEAS Compute the three tetrahedron quality measures. EMNRTH: eigenvalue measurement; RADRTH: inradius / circumradius; SANGMN: solid angle measurement A = 0.0000 0.0000 0.0000 B = 1.0000 0.0000 0.0000 C = 0.0000 0.0255 0.3525 D = 0.6669 0.9631 0.8383 Sigma = 0.3965910E-01 Rho = 0.3636671E+00 Eta = 0.4619737E+00 Q: 0.3688522E+01 0.1263040E+00 0.2998713E+00 U: 0.6489947E+00 0.7077502E-01 0.6576440E-01 A = 0.0000 0.0000 0.0000 B = 1.0000 0.0000 0.0000 C = 0.3354 0.9153 0.7959 D = 0.8327 0.3450 0.8712 Sigma = 0.1197557E+00 Rho = 0.5724350E+00 Eta = 0.6495986E+00 Q: 0.2088291E+01 0.2287329E+00 0.3654634E+00 U: 0.7911200E+00 0.1655055E+00 0.1582826E+00 A = 0.0000 0.0000 0.0000 B = 1.0000 0.0000 0.0000 C = 0.0899 0.8883 0.7010 D = 0.7346 0.3002 0.0497 Sigma = 0.2674363E-01 Rho = 0.1341967E+00 Eta = 0.3583679E+00 Q: 0.2915778E+01 0.1246599E+00 0.1485035E+01 U: 0.2897310E+00 0.5773958E-01 0.7300448E-01 A = 0.0000 0.0000 0.0000 B = 1.0000 0.0000 0.0000 C = 0.9082 0.0977 0.0403 D = 0.0850 0.5588 0.9265 Sigma = 0.9959059E-02 Rho = 0.6250437E-01 Eta = 0.1874096E+00 Q: 0.9495877E+01 0.1227526E+00 0.2549163E+01 U: 0.2194404E+00 0.3496427E-01 0.3983485E-01 A = 0.0000 0.0000 0.0000 B = 1.0000 0.0000 0.0000 C = 0.0756 0.9118 0.0918 D = 0.6378 0.8523 0.1211 Sigma = 0.8748222E-02 Rho = 0.1125599E-01 Eta = 0.1299517E+00 Q: 0.5129066E+01 0.1867442E+00 0.6904817E+02 U: 0.5200528E-01 0.4041880E-01 0.8245707E-01 Lower bounds on Q: 0.1000000E+01 0.1227526E+00 0.2998713E+00 Upper bounds on U: 0.7911200E+00 0.1655055E+00 0.1582826E+00 TEST_PRIME PRIME returns a prime number greater than a given value I. However, its largest stored prime is 14011. I PRIME(I) 14919 14011 11674 11813 295 307 2548 2557 14922 14011 11306 11411 3270 3299 9890 10037 8086 8219 7201 7211 TEST_ROTIAR ROTIAR shifts an integer array. Initial vector: 1 2 3 4 5 6 7 8 9 10 Vector after applying a shift of  4 5 6 7 8 9 10 1 2 3 GEOMPACK3_PRB: Normal end of execution. September 8 2005 8:26:28.049 AM