(pa, pb, pc, pd, pe)
| 3345 | } |
| 3346 | |
| 3347 | REAL insphereexact(pa, pb, pc, pd, pe) |
| 3348 | REAL *pa; |
| 3349 | REAL *pb; |
| 3350 | REAL *pc; |
| 3351 | REAL *pd; |
| 3352 | REAL *pe; |
| 3353 | { |
| 3354 | INEXACT REAL axby1, bxcy1, cxdy1, dxey1, exay1; |
| 3355 | INEXACT REAL bxay1, cxby1, dxcy1, exdy1, axey1; |
| 3356 | INEXACT REAL axcy1, bxdy1, cxey1, dxay1, exby1; |
| 3357 | INEXACT REAL cxay1, dxby1, excy1, axdy1, bxey1; |
| 3358 | REAL axby0, bxcy0, cxdy0, dxey0, exay0; |
| 3359 | REAL bxay0, cxby0, dxcy0, exdy0, axey0; |
| 3360 | REAL axcy0, bxdy0, cxey0, dxay0, exby0; |
| 3361 | REAL cxay0, dxby0, excy0, axdy0, bxey0; |
| 3362 | REAL ab[4], bc[4], cd[4], de[4], ea[4]; |
| 3363 | REAL ac[4], bd[4], ce[4], da[4], eb[4]; |
| 3364 | REAL temp8a[8], temp8b[8], temp16[16]; |
| 3365 | int temp8alen, temp8blen, temp16len; |
| 3366 | REAL abc[24], bcd[24], cde[24], dea[24], eab[24]; |
| 3367 | REAL abd[24], bce[24], cda[24], deb[24], eac[24]; |
| 3368 | int abclen, bcdlen, cdelen, dealen, eablen; |
| 3369 | int abdlen, bcelen, cdalen, deblen, eaclen; |
| 3370 | REAL temp48a[48], temp48b[48]; |
| 3371 | int temp48alen, temp48blen; |
| 3372 | REAL abcd[96], bcde[96], cdea[96], deab[96], eabc[96]; |
| 3373 | int abcdlen, bcdelen, cdealen, deablen, eabclen; |
| 3374 | REAL temp192[192]; |
| 3375 | REAL det384x[384], det384y[384], det384z[384]; |
| 3376 | int xlen, ylen, zlen; |
| 3377 | REAL detxy[768]; |
| 3378 | int xylen; |
| 3379 | REAL adet[1152], bdet[1152], cdet[1152], ddet[1152], edet[1152]; |
| 3380 | int alen, blen, clen, dlen, elen; |
| 3381 | REAL abdet[2304], cddet[2304], cdedet[3456]; |
| 3382 | int ablen, cdlen; |
| 3383 | REAL deter[5760]; |
| 3384 | int deterlen; |
| 3385 | int i; |
| 3386 | |
| 3387 | INEXACT REAL bvirt; |
| 3388 | REAL avirt, bround, around; |
| 3389 | INEXACT REAL c; |
| 3390 | INEXACT REAL abig; |
| 3391 | REAL ahi, alo, bhi, blo; |
| 3392 | REAL err1, err2, err3; |
| 3393 | INEXACT REAL _i, _j; |
| 3394 | REAL _0; |
| 3395 | |
| 3396 | Two_Product(pa[0], pb[1], axby1, axby0); |
| 3397 | Two_Product(pb[0], pa[1], bxay1, bxay0); |
| 3398 | Two_Two_Diff(axby1, axby0, bxay1, bxay0, ab[3], ab[2], ab[1], ab[0]); |
| 3399 | |
| 3400 | Two_Product(pb[0], pc[1], bxcy1, bxcy0); |
| 3401 | Two_Product(pc[0], pb[1], cxby1, cxby0); |
| 3402 | Two_Two_Diff(bxcy1, bxcy0, cxby1, cxby0, bc[3], bc[2], bc[1], bc[0]); |
| 3403 | |
| 3404 | Two_Product(pc[0], pd[1], cxdy1, cxdy0); |
no test coverage detected