buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,1+s_offset_b1)
