来自重点研发计划课题《面向海洋环境预测领域的高效求解技术与应用计算软件》对WRF软件的中间的分析过程。编号:2025YFB3003605。
C/C++语言由于其格式的自由,一行代码实际上并不真的是一行代码
1)类似这种
#define MAX(a,b) ( ((a) >= (b)) ?(a):(b) ) #define MIN(b,c) ( ((b) < (c)) ?(b):(c) )这个define算是一行,因为后面行跟它本身是独立的
2)
int Nfun,Njac,Nstp,Nacc,Nrej,Ndec,Nsol,Nsng;这个是真正的独立代码行
3)
#define DeltaMin (KPP_REAL)1.0e-6 /*~~~> Collect statistics: global variables */ int Nfun,Njac,Nstp,Nacc,Nrej,Ndec,Nsol,Nsng;上面这个其实真正有效的代码行只有两行
空白行和注释行在程序执行层面都没有意义。
4)
/*~~~> Function headers */ void FunTemplate(KPP_REAL, KPP_REAL [], KPP_REAL []); void JacTemplate(KPP_REAL, KPP_REAL [], KPP_REAL []) ; int Rosenbrock(KPP_REAL Y[], KPP_REAL Tstart, KPP_REAL Tend, KPP_REAL AbsTol[], KPP_REAL RelTol[], void (*ode_Fun)(KPP_REAL, KPP_REAL [], KPP_REAL []), void (*ode_Jac)(KPP_REAL, KPP_REAL [], KPP_REAL []), KPP_REAL RPAR[], int IPAR[]); int RosenbrockIntegrator( KPP_REAL Y[], KPP_REAL Tstart, KPP_REAL Tend , KPP_REAL AbsTol[], KPP_REAL RelTol[], void (*ode_Fun)(KPP_REAL, KPP_REAL [], KPP_REAL []), void (*ode_Jac)(KPP_REAL, KPP_REAL [], KPP_REAL []), int ros_S, KPP_REAL ros_M[], KPP_REAL ros_E[], KPP_REAL ros_A[], KPP_REAL ros_C[], KPP_REAL ros_Alpha[],KPP_REAL ros_Gamma[], KPP_REAL ros_ELO, char ros_NewF[], char Autonomous, char VectorTol, int Max_no_steps, KPP_REAL Roundoff, KPP_REAL Hmin, KPP_REAL Hmax, KPP_REAL Hstart, KPP_REAL FacMin, KPP_REAL FacMax, KPP_REAL FacRej, KPP_REAL FacSafe, KPP_REAL *Texit, KPP_REAL *Hexit );上面实际上只有4行代码,只是有些参数有点多。
5)
int a; double b; float c;这个其实是三行代码
6)
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/ void INTEGRATE( KPP_REAL TIN, KPP_REAL TOUT ) /*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/ { static KPP_REAL RPAR[20];这个其实是二行代码
IPAR[3] = 5; /* choice of the method */ IERR = Rosenbrock(VAR, TIN, TOUT, ATOL, RTOL, &FunTemplate, &JacTemplate, RPAR, IPAR); Ns=Ns+IPAR[12];这个其实是三行代码。
printf("\n Step=%d Acc=%d Rej=%d Singular=%d\n", Ns,Na,Nr,Ng); if (IERR < 0) printf("\n Rosenbrock: Unsucessful step at T=%g: IERR=%d\n", TIN,IERR);这个是二行代码。
7)
if (RPAR[1] == ZERO) Hmax = ABS(Tend-Tstart); else Hmax = MIN(ABS(RPAR[1]),ABS(Tend-Tstart));这个又复杂了一点。
8)
if (RPAR[2] == ZERO) Hstart = MAX(Hmin,DeltaMin); else Hstart = MIN(ABS(RPAR[2]),ABS(Tend-Tstart)); if (RPAR[2] == ZERO) Hstart = MAX(Hmin,DeltaMin); else Hstart = MIN(ABS(RPAR[2]),ABS(Tend-Tstart));这两个是正确的格式吗?可能跟编译器相关。
if (RPAR[4] == ZERO) FacMax = (KPP_REAL)6.0; //这个应该是可以的9)
switch (Method) { case 1: Ros2(&ros_S, ros_A, ros_C, ros_M, ros_E, ros_Alpha, ros_Gamma, ros_NewF, &ros_ELO, ros_Name); break;case 1这个应该算独立一行
10)总结 下,
C代码的分行,实际 上以
;
:
{}
等为分界。