@@ -66,11 +66,11 @@ const int rungekutta_s = 4;
6666const double rungekutta_b [4 ] = { 1.0 / 6.0 , 1.0 / 3.0 , 1.0 / 3.0 , 1.0 / 6.0 };
6767const double rungekutta_c [4 ] = { 0.0 , 0.5 , 0.5 , 1.0 };
6868
69- typedef struct RK4
69+ typedef struct RK4_DATA
7070{
7171 double * * work_states ;
7272 int work_states_ndims ;
73- }RK4 ;
73+ }RK4_DATA ;
7474
7575
7676static int euler_ex_step (DATA * data , SOLVER_INFO * solverInfo );
@@ -153,7 +153,7 @@ int initializeSolverData(DATA* data, SOLVER_INFO* solverInfo)
153153 if (solverInfo -> solverMethod == 2 )
154154 {
155155 /* Allocate RK work arrays */
156- RK4 * rungeData = (RK4 * ) malloc (sizeof (RK4 ));
156+ RK4_DATA * rungeData = (RK4_DATA * ) malloc (sizeof (RK4_DATA ));
157157 rungeData -> work_states_ndims = rungekutta_s ;
158158 rungeData -> work_states = (double * * ) malloc ((rungeData -> work_states_ndims + 1 ) * sizeof (double * ));
159159 for (i = 0 ; i < rungeData -> work_states_ndims + 1 ; i ++ )
@@ -246,10 +246,10 @@ int freeSolverData(DATA* data, SOLVER_INFO* solverInfo)
246246 if (solverInfo -> solverMethod == 2 )
247247 {
248248 /* free RK work arrays */
249- for (i = 0 ; i < ((RK4 * )(solverInfo -> solverData ))-> work_states_ndims + 1 ; i ++ )
250- free (((RK4 * )(solverInfo -> solverData ))-> work_states [i ]);
251- free (((RK4 * )(solverInfo -> solverData ))-> work_states );
252- free ((RK4 * )solverInfo -> solverData );
249+ for (i = 0 ; i < ((RK4_DATA * )(solverInfo -> solverData ))-> work_states_ndims + 1 ; i ++ )
250+ free (((RK4_DATA * )(solverInfo -> solverData ))-> work_states [i ]);
251+ free (((RK4_DATA * )(solverInfo -> solverData ))-> work_states );
252+ free ((RK4_DATA * )solverInfo -> solverData );
253253 }
254254 else if (solverInfo -> solverMethod == 3 )
255255 {
@@ -620,7 +620,7 @@ static int euler_ex_step(DATA* data, SOLVER_INFO* solverInfo)
620620/*************************************** RK4 ***********************************/
621621static int rungekutta_step (DATA * data , SOLVER_INFO * solverInfo )
622622{
623- double * * k = ((RK4 * )(solverInfo -> solverData ))-> work_states ;
623+ double * * k = ((RK4_DATA * )(solverInfo -> solverData ))-> work_states ;
624624 double sum ;
625625 int i ,j ;
626626 SIMULATION_DATA * sData = (SIMULATION_DATA * )data -> localData [0 ];
0 commit comments