翻译成C语言 SUBROUTINE PREBINDERFBHZIXXXFHXNYNZNXHYHIIII & IIIIHL0L0HNENEHNKNCHE IMPLICIT REAL8A-HO-Z DIMENSION FBHZL0IXXXL0FHNEHXNL0YNL0ZNL0XHL0H &

void PREBINDER(double FBHZ, int IXXX, double FH, double XN, double YN, double ZN, double XH, double YH, int IIII, int IIIIH, int L0, int L0H, int NE, int NEH, int NK, int NCHE) { double XXH[4], YYH[4], IXH[4], XX[4], YY[4], ZZ[4], IX[4]; int I00=0, IXYD=0; double F00=0.0, XA, YA, EPS=1.0e-6, XMAX, XMIN, YMAX, YMIN, DX21, DY21, DZ21, DX31, DY31, DZ31, DX41, DY41, DZ41, XD1, XD2, XD3, AREAE1, AREAE2, AREAE, ANN, BHF3, BHF4; for(int I0=0; I0<NEH; I0++) { for(int JJ=0; JJ<NK; JJ++) { IXH[JJ] = IIIIH[I0NK+JJ]; XXH[JJ] = XH[IXH[JJ]]; YYH[JJ] = YH[IXH[JJ]]; } if(IXH[2]==IXH[3] || IXH[3]==0) { XA = (XXH[0]+XXH[1]+XXH[2])/3.0; YA = (YYH[0]+YYH[1]+YYH[2])/3.0; } else { XA = (XXH[0]+XXH[1]+XXH[2]+XXH[3])/4.0; YA = (YYH[0]+YYH[1]+YYH[2]+YYH[3])/4.0; } for(int II=0; II<NE; II++) { for(int JJ=0; JJ<NK; JJ++) { IX[JJ] = IIII[IINK+JJ]; XX[JJ] = XN[IX[JJ]]; YY[JJ] = YN[IX[JJ]]; ZZ[JJ] = ZN[IX[JJ]]; } } IXYD=0; XMAX=XX[0]; XMIN=XX[0]; YMAX=YY[0]; YMIN=YY[0]; for(int I=1; I<NK; I++) { if(XMAX<XX[I]) XMAX=XX[I]; if(YMAX<YY[I]) YMAX=YY[I]; if(XMIN>XX[I]) XMIN=XX[I]; if(YMIN>YY[I]) YMIN=YY[I]; } if(XA<XMAX+EPS && XA>XMIN-EPS) { if(YA<YMAX+EPS && YA>YMIN-EPS) { IXYD=1; } } if(IXYD==1) PREXYB(&IXYD, XX, YY, &XA, &YA, NK); if(IXYD==1 || IXYD==2) goto label_600; label_200: I00++; F00+=FH[I0]; } label_600: if(IXYD==1 || IXYD==2) { DX21=XX[1]-XX[0]; DY21=YY[1]-YY[0]; DZ21=ZZ[1]-ZZ[0]; DX31=XX[2]-XX[0]; DY31=YY[2]-YY[0]; DZ31=ZZ[2]-ZZ[0]; DX41=XX[3]-XX[0]; DY41=YY[3]-YY[0]; DZ41=ZZ[3]-ZZ[0]; XD1=DY21DZ31-DY31DZ21; XD2=DZ21DX31-DZ31DX21; XD3=DX21DY31-DX31DY21; AREAE1=0.5sqrt(XD1XD1+XD2XD2+XD3XD3+1.0e-80); XD1=DY41DZ31-DY31DZ41; XD2=DZ41DX31-DZ31DX41; XD3=DX41DY31-DX31DY41; AREAE2=0.5sqrt(XD1XD1+XD2XD2+XD3XD3+1.0e-80); AREAE=AREAE1+AREAE2; if(IX[2]==IX[3] || IX[3]==0) { ANN=3.0; BHF3=FH[I0]/ANN; for(int I=0; I<3; I++) { FBHZ[IX[I]]+=BHF3; IXXX[IX[I]]=3; } } else { ANN=4.0; BHF4=FH[I0]/ANN; for(int I=0; I<NK; I++) { FBHZ[IX[I]]+=BHF4; IXXX[IX[I]]=3; } } }