转换成C代码 IFKKL3NEKKL4 THEN DO 108 IK2=14 KKL1=KKLIK2 NN5=0 DO 107 IKK=14 IKK1=NN4+IKK IFLIDBIKK1EQKKL1 THEN FIDBIKK1=FIDBIKK1

if (KKL[2] != KKL[3]) { for (int IK2 = 0; IK2 < 4; IK2++) { int KKL1 = KKL[IK2]; int NN5 = 0; for (int IKK = 0; IKK < 4; IKK++) { int IKK1 = NN4 + IKK; if (LIDB[IKK1] == KKL1) { FIDB[IKK1] += FORCE2; } if (LIDB[IKK1] != KKL1) { NN5 += 1; } } if (NN5 == 4) { NN6 += 1; FIDB[NN4 + NN6] += FORCE2; LIDB[NN4 + NN6] = KKL1; if (IXXX[KKL1] == 1 || IXXX[KKL1] == 2) { IXXX[KKL1] = 5; } if (IXXX[KKL1] == 3) { IXXX[KKL1] = 4; } } } } for (int II = 0; II < NE; II++) { int N1 = 0; int I1 = IIII[II][0]; int I2 = IIII[II][1]; int I3 = IIII[II][2]; int I4 = IIII[II][3]; double XX[4] = {XN[I1], XN[I2], XN[I3], XN[I4]}; double YY[4] = {YN[I1], YN[I2], YN[I3], YN[I4]}; if (fabs(ALX2 - ALX1) >= EPSILON || fabs(XX[1] - XX[0]) >= EPSILON) { if (fabs(ALX2 - ALX1) < EPSILON && fabs(XX[1] - XX[0]) >= EPSILON) { double AK12 = (YY[1] - YY[0]) / (XX[1] - XX[0]); double B12 = YY[0] - AK12 * XX[0]; double XX11 = ALX1; double YY11 = AK12 * XX11 + B12; double VALUE1 = (YY11 - YY[1]) * (YY11 - YY[0]); double VALUE2 = (YY11 - ALY1) * (YY11 - ALY2); double VALUE3 = (XX11 - XX[1]) * (XX11 - XX[0]); if (VALUE1 <= 0 && VALUE2 <= 0 && VALUE3 <= 0) { N1 += 1; B[N1][0] = XX11; B[N1][1] = YY11; } } if (fabs(XX[1] - XX[0]) < EPSILON && fabs(ALX2 - ALX1) >= EPSILON) { double AK = (ALY2 - ALY1) / (ALX2 - ALX1); double BB = ALY1 - AK * ALX1; double XX12 = XX[0]; double YY12 = AK * XX12 + BB; double VALUE1 = (YY12 - YY[1]) * (YY12 - YY[0]); double VALUE2 = (YY12 - ALY1) * (YY12 - ALY2); double VALUE3 = (XX12 - ALX1) * (XX12 - ALX2); if (VALUE1 <= 0 && VALUE2 <= 0 && VALUE3 <= 0) { N1 += 1; B[N1][0] = XX12; B[N1][1] = YY12; } } if (fabs(ALX2 - ALX1) >= EPSILON && fabs(XX[1] - XX[0]) >= EPSILON) { double AK12 = (YY[1] - YY[0]) / (XX[1] - XX[0]); double AK = (ALY2 - ALY1) / (ALX2 - ALX1); if (fabs(AK12 - AK) >= EPSILON) { double B12 = YY[0] - AK12 * XX[0]; double BB = ALY1 - AK * ALX1; double XX13 = (BB - B12) / (AK12 - AK); double T = XX13; double YY13 = AK * T + BB; double VALUE1 = (XX13 - XX[1]) * (XX13 - XX[0]); double VALUE2 = (YY13 - ALY1) * (YY13 - ALY2); double VALUE3 = (XX13 - ALX1) * (XX13 - ALX2); if (VALUE1 <= 0 && VALUE2 <= 0 && VALUE3 <= 0) { N1 += 1; B[N1][0] = XX13; B[N1][1] = YY13; } } } }