转换成C代码 DO I=14 IFIIDBX1KKIEQIIDBX2KKI NN=NN+1 ENDDO IFNNEQ4 THEN T1=DBX2KK-DBX1KK T2=DBY2KK-DBY1KK T3=DBZ2KK-DBZ1KK ALENGTH=DSQRTT1T1+T2T2+T3T3

int I; int NN = 0; double T1, T2, T3, ALENGTH, FIDBT; int NN1 = 0; int I1, I2, I3, I4;

for (I = 1; I <= 4; I++) { if (IIDBX1[KK][I] == IIDBX2[KK][I]) { NN++; } } if (NN == 4) { T1 = DBX2[KK] - DBX1[KK]; T2 = DBY2[KK] - DBY1[KK]; T3 = DBZ2[KK] - DBZ1[KK]; ALENGTH = sqrt(T1T1 + T2T2 + T3*T3); if (IIDBX1[KK][3] != IIDBX1[KK][4]) { FIDBT = DBF[KK] * ALENGTH / 4.0; } else { FIDBT = DBF[KK] * ALENGTH / 3.0; } for (I = 1; I < KK; I++) { NN1 += LDB[I]; } LIDB[NN1+1] = IIDBX1[KK][1]; LIDB[NN1+2] = IIDBX1[KK][2]; LIDB[NN1+3] = IIDBX1[KK][3]; if (IIDBX1[KK][3] != IIDBX1[KK][4]) { LIDB[NN1+4] = IIDBX1[KK][4]; } FIDB[NN1+1] += FIDBT; FIDB[NN1+2] += FIDBT; FIDB[NN1+3] += FIDBT; if (IIDBX1[KK][3] != IIDBX1[KK][4]) { FIDB[NN1+4] += FIDBT; } I1 = IIDBX1[KK][1]; I2 = IIDBX1[KK][2]; I3 = IIDBX1[KK][3]; I4 = IIDBX1[KK][4]; if (IXXX[I1] == 1 || IXXX[I1] == 2) { IXXX[I1] = 5; } if (IXXX[I1] == 3) { IXXX[I1] = 4; } if (IXXX[I2] == 1 || IXXX[I2] == 2) { IXXX[I2] = 5; } if (IXXX[I2] == 3) { IXXX[I2] = 4; } if (IXXX[I3] == 1 || IXXX[I3] == 2) { IXXX[I3] = 5; } if (IXXX[I3] == 3) { IXXX[I3] = 4; } if (IXXX[I4] == 1 || IXXX[I4] == 2) { IXXX[I4] = 5; } if (IXXX[I4] == 3) { IXXX[I4] = 4; }