将Fortran代码转化为C代码 SUBROUTINE AnglecosVectorCosangleVector1Vector2 IMPLICIT NONE REAL8 CosangleVector13Vector23abdjaabb abdj=Vector11Vector21+Vector12Vector22+ + Vector13Vector23 aa=DSQRTVecto

void AnglecosVector(double *Cosangle, double Vector1[3], double Vector2[3]) { double abdj, aa, bb; abdj = Vector1[0]*Vector2[0] + Vector1[1]*Vector2[1] + Vector1[2]*Vector2[2]; aa = sqrt(Vector1[0]*Vector1[0] + Vector1[1]*Vector1[1] + Vector1[2]*Vector1[2]); bb = sqrt(Vector2[0]*Vector2[0] + Vector2[1]*Vector2[1] + Vector2[2]*Vector2[2]); Cosangle = abdj / (aabb + 1.0e-50); }