 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2( 1 +(i-1)*kmax) * sphi_c( 1 , 7 +s_offset_c1)
 buffer1(i+imax*( 9 -1)) = buffer1(i+imax*( 9 -1)) + buffer2( 1 +(i-1)*kmax) * sphi_c( 1 , 9 +s_offset_c1)
 buffer1(i+imax*(11 -1)) = buffer1(i+imax*(11 -1)) + buffer2( 1 +(i-1)*kmax) * sphi_c( 1 ,11 +s_offset_c1)
 buffer1(i+imax*( 1 -1)) = buffer1(i+imax*( 1 -1)) + buffer2( 2 +(i-1)*kmax) * sphi_c( 2 , 1 +s_offset_c1)
 buffer1(i+imax*( 3 -1)) = buffer1(i+imax*( 3 -1)) + buffer2( 2 +(i-1)*kmax) * sphi_c( 2 , 3 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2( 2 +(i-1)*kmax) * sphi_c( 2 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2( 3 +(i-1)*kmax) * sphi_c( 3 , 6 +s_offset_c1)
 buffer1(i+imax*( 8 -1)) = buffer1(i+imax*( 8 -1)) + buffer2( 3 +(i-1)*kmax) * sphi_c( 3 , 8 +s_offset_c1)
 buffer1(i+imax*(10 -1)) = buffer1(i+imax*(10 -1)) + buffer2( 3 +(i-1)*kmax) * sphi_c( 3 ,10 +s_offset_c1)
 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2( 4 +(i-1)*kmax) * sphi_c( 4 , 7 +s_offset_c1)
 buffer1(i+imax*( 9 -1)) = buffer1(i+imax*( 9 -1)) + buffer2( 4 +(i-1)*kmax) * sphi_c( 4 , 9 +s_offset_c1)
 buffer1(i+imax*(11 -1)) = buffer1(i+imax*(11 -1)) + buffer2( 4 +(i-1)*kmax) * sphi_c( 4 ,11 +s_offset_c1)
 buffer1(i+imax*( 2 -1)) = buffer1(i+imax*( 2 -1)) + buffer2( 5 +(i-1)*kmax) * sphi_c( 5 , 2 +s_offset_c1)
 buffer1(i+imax*( 4 -1)) = buffer1(i+imax*( 4 -1)) + buffer2( 5 +(i-1)*kmax) * sphi_c( 5 , 4 +s_offset_c1)
 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2( 6 +(i-1)*kmax) * sphi_c( 6 , 7 +s_offset_c1)
 buffer1(i+imax*( 9 -1)) = buffer1(i+imax*( 9 -1)) + buffer2( 6 +(i-1)*kmax) * sphi_c( 6 , 9 +s_offset_c1)
 buffer1(i+imax*( 1 -1)) = buffer1(i+imax*( 1 -1)) + buffer2( 7 +(i-1)*kmax) * sphi_c( 7 , 1 +s_offset_c1)
 buffer1(i+imax*( 3 -1)) = buffer1(i+imax*( 3 -1)) + buffer2( 7 +(i-1)*kmax) * sphi_c( 7 , 3 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2( 7 +(i-1)*kmax) * sphi_c( 7 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2( 8 +(i-1)*kmax) * sphi_c( 8 , 6 +s_offset_c1)
 buffer1(i+imax*(10 -1)) = buffer1(i+imax*(10 -1)) + buffer2( 8 +(i-1)*kmax) * sphi_c( 8 ,10 +s_offset_c1)
 buffer1(i+imax*( 3 -1)) = buffer1(i+imax*( 3 -1)) + buffer2( 9 +(i-1)*kmax) * sphi_c( 9 , 3 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2( 9 +(i-1)*kmax) * sphi_c( 9 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2(10 +(i-1)*kmax) * sphi_c(10 , 6 +s_offset_c1)
 buffer1(i+imax*( 8 -1)) = buffer1(i+imax*( 8 -1)) + buffer2(10 +(i-1)*kmax) * sphi_c(10 , 8 +s_offset_c1)
 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2(11 +(i-1)*kmax) * sphi_c(11 , 7 +s_offset_c1)
 buffer1(i+imax*( 9 -1)) = buffer1(i+imax*( 9 -1)) + buffer2(11 +(i-1)*kmax) * sphi_c(11 , 9 +s_offset_c1)
 buffer1(i+imax*(11 -1)) = buffer1(i+imax*(11 -1)) + buffer2(11 +(i-1)*kmax) * sphi_c(11 ,11 +s_offset_c1)
 buffer1(i+imax*( 2 -1)) = buffer1(i+imax*( 2 -1)) + buffer2(12 +(i-1)*kmax) * sphi_c(12 , 2 +s_offset_c1)
 buffer1(i+imax*( 4 -1)) = buffer1(i+imax*( 4 -1)) + buffer2(12 +(i-1)*kmax) * sphi_c(12 , 4 +s_offset_c1)
 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2(13 +(i-1)*kmax) * sphi_c(13 , 7 +s_offset_c1)
 buffer1(i+imax*( 9 -1)) = buffer1(i+imax*( 9 -1)) + buffer2(13 +(i-1)*kmax) * sphi_c(13 , 9 +s_offset_c1)
 buffer1(i+imax*( 4 -1)) = buffer1(i+imax*( 4 -1)) + buffer2(14 +(i-1)*kmax) * sphi_c(14 , 4 +s_offset_c1)
 buffer1(i+imax*( 7 -1)) = buffer1(i+imax*( 7 -1)) + buffer2(15 +(i-1)*kmax) * sphi_c(15 , 7 +s_offset_c1)
 buffer1(i+imax*( 1 -1)) = buffer1(i+imax*( 1 -1)) + buffer2(16 +(i-1)*kmax) * sphi_c(16 , 1 +s_offset_c1)
 buffer1(i+imax*( 3 -1)) = buffer1(i+imax*( 3 -1)) + buffer2(16 +(i-1)*kmax) * sphi_c(16 , 3 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2(16 +(i-1)*kmax) * sphi_c(16 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2(17 +(i-1)*kmax) * sphi_c(17 , 6 +s_offset_c1)
 buffer1(i+imax*( 8 -1)) = buffer1(i+imax*( 8 -1)) + buffer2(17 +(i-1)*kmax) * sphi_c(17 , 8 +s_offset_c1)
 buffer1(i+imax*(10 -1)) = buffer1(i+imax*(10 -1)) + buffer2(17 +(i-1)*kmax) * sphi_c(17 ,10 +s_offset_c1)
 buffer1(i+imax*( 3 -1)) = buffer1(i+imax*( 3 -1)) + buffer2(18 +(i-1)*kmax) * sphi_c(18 , 3 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2(18 +(i-1)*kmax) * sphi_c(18 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2(19 +(i-1)*kmax) * sphi_c(19 , 6 +s_offset_c1)
 buffer1(i+imax*( 8 -1)) = buffer1(i+imax*( 8 -1)) + buffer2(19 +(i-1)*kmax) * sphi_c(19 , 8 +s_offset_c1)
 buffer1(i+imax*( 5 -1)) = buffer1(i+imax*( 5 -1)) + buffer2(20 +(i-1)*kmax) * sphi_c(20 , 5 +s_offset_c1)
 buffer1(i+imax*( 6 -1)) = buffer1(i+imax*( 6 -1)) + buffer2(21 +(i-1)*kmax) * sphi_c(21 , 6 +s_offset_c1)
