(pa, pb, pc, pd, pe, permanent)
| 3936 | } |
| 3937 | |
| 3938 | REAL insphereadapt(pa, pb, pc, pd, pe, permanent) |
| 3939 | REAL *pa; |
| 3940 | REAL *pb; |
| 3941 | REAL *pc; |
| 3942 | REAL *pd; |
| 3943 | REAL *pe; |
| 3944 | REAL permanent; |
| 3945 | { |
| 3946 | INEXACT REAL aex, bex, cex, dex, aey, bey, cey, dey, aez, bez, cez, dez; |
| 3947 | REAL det, errbound; |
| 3948 | |
| 3949 | INEXACT REAL aexbey1, bexaey1, bexcey1, cexbey1; |
| 3950 | INEXACT REAL cexdey1, dexcey1, dexaey1, aexdey1; |
| 3951 | INEXACT REAL aexcey1, cexaey1, bexdey1, dexbey1; |
| 3952 | REAL aexbey0, bexaey0, bexcey0, cexbey0; |
| 3953 | REAL cexdey0, dexcey0, dexaey0, aexdey0; |
| 3954 | REAL aexcey0, cexaey0, bexdey0, dexbey0; |
| 3955 | REAL ab[4], bc[4], cd[4], da[4], ac[4], bd[4]; |
| 3956 | INEXACT REAL ab3, bc3, cd3, da3, ac3, bd3; |
| 3957 | REAL abeps, bceps, cdeps, daeps, aceps, bdeps; |
| 3958 | REAL temp8a[8], temp8b[8], temp8c[8], temp16[16], temp24[24], temp48[48]; |
| 3959 | int temp8alen, temp8blen, temp8clen, temp16len, temp24len, temp48len; |
| 3960 | REAL xdet[96], ydet[96], zdet[96], xydet[192]; |
| 3961 | int xlen, ylen, zlen, xylen; |
| 3962 | REAL adet[288], bdet[288], cdet[288], ddet[288]; |
| 3963 | int alen, blen, clen, dlen; |
| 3964 | REAL abdet[576], cddet[576]; |
| 3965 | int ablen, cdlen; |
| 3966 | REAL fin1[1152]; |
| 3967 | int finlength; |
| 3968 | |
| 3969 | REAL aextail, bextail, cextail, dextail; |
| 3970 | REAL aeytail, beytail, ceytail, deytail; |
| 3971 | REAL aeztail, beztail, ceztail, deztail; |
| 3972 | |
| 3973 | INEXACT REAL bvirt; |
| 3974 | REAL avirt, bround, around; |
| 3975 | INEXACT REAL c; |
| 3976 | INEXACT REAL abig; |
| 3977 | REAL ahi, alo, bhi, blo; |
| 3978 | REAL err1, err2, err3; |
| 3979 | INEXACT REAL _i, _j; |
| 3980 | REAL _0; |
| 3981 | |
| 3982 | aex = (REAL) (pa[0] - pe[0]); |
| 3983 | bex = (REAL) (pb[0] - pe[0]); |
| 3984 | cex = (REAL) (pc[0] - pe[0]); |
| 3985 | dex = (REAL) (pd[0] - pe[0]); |
| 3986 | aey = (REAL) (pa[1] - pe[1]); |
| 3987 | bey = (REAL) (pb[1] - pe[1]); |
| 3988 | cey = (REAL) (pc[1] - pe[1]); |
| 3989 | dey = (REAL) (pd[1] - pe[1]); |
| 3990 | aez = (REAL) (pa[2] - pe[2]); |
| 3991 | bez = (REAL) (pb[2] - pe[2]); |
| 3992 | cez = (REAL) (pc[2] - pe[2]); |
| 3993 | dez = (REAL) (pd[2] - pe[2]); |
| 3994 | |
| 3995 | Two_Product(aex, bey, aexbey1, aexbey0); |
no test coverage detected