FORTRAN SUBROUTINE FORMJB3U to C Code Conversion

#include <stdio.h>

void FORMJB3U(double *BBP, double *CCP, double *AJF, double *DAJ, double *XT, double *YT) { double AJ[2][2]; AJ[0][0] = -0.25 * XT[1]; AJ[0][1] = 0.0; AJ[1][0] = 0.25 * (XT[1] - 2.0 * XT[2]); AJ[1][1] = -0.5 * YT[2]; *DAJ = AJ[0][0] * AJ[1][1]; AJF[0] = AJ[1][1] / (*DAJ); AJF[1] = 0.0; AJF[2] = -AJ[1][0] / (*DAJ); AJF[3] = AJ[0][0] / (*DAJ); BBP[0] = 0.25 * AJF[0]; BBP[1] = -0.25 * AJF[0]; BBP[2] = 0.0; CCP[0] = 0.25 * (AJF[2] + AJF[3]); CCP[1] = 0.25 * (-AJF[2] + AJF[3]); CCP[2] = -0.5 * AJF[3]; }

int main() { double BBP[3], CCP[3], AJF[4], DAJ, XT[3], YT[3]; FORMJB3U(BBP, CCP, AJF, &DAJ, XT, YT); printf("BBP: %f %f %f\n", BBP[0], BBP[1], BBP[2]); printf("CCP: %f %f %f\n", CCP[0], CCP[1], CCP[2]); printf("AJF: %f %f %f %f\n", AJF[0], AJF[1], AJF[2], AJF[3]); printf("DAJ: %f\n", DAJ); return 0; }