 buffer1(i+imax*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(1+(i-1)*kmax) * sphi_c(1,5+s_offset_c1)
 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*(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*(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*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(4+(i-1)*kmax) * sphi_c(4,5+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*(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*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(6+(i-1)*kmax) * sphi_c(6,5+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*(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*(6-1)) = buffer1(i+imax*(6-1)) + buffer2(8+(i-1)*kmax) * sphi_c(8,6+s_offset_c1)
 buffer1(i+imax*(8-1)) = buffer1(i+imax*(8-1)) + buffer2(8+(i-1)*kmax) * sphi_c(8,8+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*(6-1)) = buffer1(i+imax*(6-1)) + buffer2(10+(i-1)*kmax) * sphi_c(10,6+s_offset_c1)
 buffer1(i+imax*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(11+(i-1)*kmax) * sphi_c(11,5+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*(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*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(13+(i-1)*kmax) * sphi_c(13,5+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*(4-1)) = buffer1(i+imax*(4-1)) + buffer2(14+(i-1)*kmax) * sphi_c(14,4+s_offset_c1)
 buffer1(i+imax*(5-1)) = buffer1(i+imax*(5-1)) + buffer2(15+(i-1)*kmax) * sphi_c(15,5+s_offset_c1)
