@@ -66,23 +66,22 @@ class MeshNumerics {
6666
6767 int i=0 ;
6868#ifdef HAVE_OPENMP
69-
7069#pragma omp parallel for default(shared) private(i)
7170#endif
7271 for (i = 0 ; i < input.z_num ; ++i) {
73- for (size_t j = 0 ; j < input.x_num ; ++j) {
74- for (size_t k = 0 ; k < input.y_num ; ++k) {
72+ for (int j = 0 ; j < input.x_num ; ++j) {
73+ for (int k = 0 ; k < input.y_num ; ++k) {
7574 double neigh_sum = 0 ;
7675 double counter = 0 ;
7776
78- int stencil_z_b = std::max ((int )( -stencil_z_half+i),( int ) 0 );
79- int stencil_z_e = std::min ((int )( stencil_z_half+i),( int ) (input.z_num -1 ));
77+ int stencil_z_b = std::max ((-stencil_z_half+i),0 );
78+ int stencil_z_e = std::min ((stencil_z_half+i),(input.z_num -1 ));
8079
81- int stencil_x_b = std::max ((int )( -stencil_x_half+j),0 );
82- int stencil_x_e = std::min ((int )( stencil_x_half+j),( int ) (input.x_num -1 ));
80+ int stencil_x_b = std::max ((-stencil_x_half+j),0 );
81+ int stencil_x_e = std::min ((stencil_x_half+j),(input.x_num -1 ));
8382
84- int stencil_y_b = std::max ((int )( -stencil_y_half+k),( int ) 0 );
85- int stencil_y_e = std::min ((int )( stencil_y_half+k),( int ) (input.y_num -1 ));
83+ int stencil_y_b = std::max ((-stencil_y_half+k),0 );
84+ int stencil_y_e = std::min ((stencil_y_half+k),(input.y_num -1 ));
8685
8786
8887 for (int l = stencil_z_b; l <= stencil_z_e; ++l) {
0 commit comments