buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(1+(i-1)*kmax) * sphi_b(1,21+s_offset_b1)
buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(2+(i-1)*kmax) * sphi_b(2,1+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(2+(i-1)*kmax) * sphi_b(2,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(2+(i-1)*kmax) * sphi_b(2,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(2+(i-1)*kmax) * sphi_b(2,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(2+(i-1)*kmax) * sphi_b(2,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(3+(i-1)*kmax) * sphi_b(3,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(3+(i-1)*kmax) * sphi_b(3,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(3+(i-1)*kmax) * sphi_b(3,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(3+(i-1)*kmax) * sphi_b(3,18+s_offset_b1)
buffer2(i+imax*(20-1)) = buffer2(i+imax*(20-1)) + buffer1(3+(i-1)*kmax) * sphi_b(3,20+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(4+(i-1)*kmax) * sphi_b(4,21+s_offset_b1)
buffer2(i+imax*(2-1)) = buffer2(i+imax*(2-1)) + buffer1(5+(i-1)*kmax) * sphi_b(5,2+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(5+(i-1)*kmax) * sphi_b(5,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(5+(i-1)*kmax) * sphi_b(5,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(5+(i-1)*kmax) * sphi_b(5,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(5+(i-1)*kmax) * sphi_b(5,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(6+(i-1)*kmax) * sphi_b(6,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(6+(i-1)*kmax) * sphi_b(6,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(6+(i-1)*kmax) * sphi_b(6,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(6+(i-1)*kmax) * sphi_b(6,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(6+(i-1)*kmax) * sphi_b(6,19+s_offset_b1)
buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(7+(i-1)*kmax) * sphi_b(7,1+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(7+(i-1)*kmax) * sphi_b(7,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(7+(i-1)*kmax) * sphi_b(7,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(7+(i-1)*kmax) * sphi_b(7,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(7+(i-1)*kmax) * sphi_b(7,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(8+(i-1)*kmax) * sphi_b(8,12+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(8+(i-1)*kmax) * sphi_b(8,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(8+(i-1)*kmax) * sphi_b(8,18+s_offset_b1)
buffer2(i+imax*(20-1)) = buffer2(i+imax*(20-1)) + buffer1(8+(i-1)*kmax) * sphi_b(8,20+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(9+(i-1)*kmax) * sphi_b(9,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(9+(i-1)*kmax) * sphi_b(9,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(9+(i-1)*kmax) * sphi_b(9,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(9+(i-1)*kmax) * sphi_b(9,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(10+(i-1)*kmax) * sphi_b(10,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(10+(i-1)*kmax) * sphi_b(10,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(10+(i-1)*kmax) * sphi_b(10,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(10+(i-1)*kmax) * sphi_b(10,18+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(11+(i-1)*kmax) * sphi_b(11,21+s_offset_b1)
buffer2(i+imax*(2-1)) = buffer2(i+imax*(2-1)) + buffer1(12+(i-1)*kmax) * sphi_b(12,2+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(12+(i-1)*kmax) * sphi_b(12,4+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(12+(i-1)*kmax) * sphi_b(12,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(12+(i-1)*kmax) * sphi_b(12,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(13+(i-1)*kmax) * sphi_b(13,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(13+(i-1)*kmax) * sphi_b(13,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(13+(i-1)*kmax) * sphi_b(13,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(13+(i-1)*kmax) * sphi_b(13,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(13+(i-1)*kmax) * sphi_b(13,19+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(14+(i-1)*kmax) * sphi_b(14,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(14+(i-1)*kmax) * sphi_b(14,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(14+(i-1)*kmax) * sphi_b(14,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(14+(i-1)*kmax) * sphi_b(14,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(15+(i-1)*kmax) * sphi_b(15,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(15+(i-1)*kmax) * sphi_b(15,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(15+(i-1)*kmax) * sphi_b(15,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(15+(i-1)*kmax) * sphi_b(15,17+s_offset_b1)
buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(16+(i-1)*kmax) * sphi_b(16,1+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(16+(i-1)*kmax) * sphi_b(16,5+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(16+(i-1)*kmax) * sphi_b(16,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(17+(i-1)*kmax) * sphi_b(17,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(17+(i-1)*kmax) * sphi_b(17,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(17+(i-1)*kmax) * sphi_b(17,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(17+(i-1)*kmax) * sphi_b(17,18+s_offset_b1)
buffer2(i+imax*(20-1)) = buffer2(i+imax*(20-1)) + buffer1(17+(i-1)*kmax) * sphi_b(17,20+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(18+(i-1)*kmax) * sphi_b(18,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(18+(i-1)*kmax) * sphi_b(18,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(18+(i-1)*kmax) * sphi_b(18,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(18+(i-1)*kmax) * sphi_b(18,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(19+(i-1)*kmax) * sphi_b(19,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(19+(i-1)*kmax) * sphi_b(19,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(19+(i-1)*kmax) * sphi_b(19,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(19+(i-1)*kmax) * sphi_b(19,18+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(20+(i-1)*kmax) * sphi_b(20,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(20+(i-1)*kmax) * sphi_b(20,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(20+(i-1)*kmax) * sphi_b(20,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(21+(i-1)*kmax) * sphi_b(21,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(21+(i-1)*kmax) * sphi_b(21,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(21+(i-1)*kmax) * sphi_b(21,16+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(22+(i-1)*kmax) * sphi_b(22,21+s_offset_b1)
buffer2(i+imax*(2-1)) = buffer2(i+imax*(2-1)) + buffer1(23+(i-1)*kmax) * sphi_b(23,2+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(23+(i-1)*kmax) * sphi_b(23,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(23+(i-1)*kmax) * sphi_b(23,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(23+(i-1)*kmax) * sphi_b(23,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(23+(i-1)*kmax) * sphi_b(23,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(24+(i-1)*kmax) * sphi_b(24,11+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(24+(i-1)*kmax) * sphi_b(24,15+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(24+(i-1)*kmax) * sphi_b(24,19+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(25+(i-1)*kmax) * sphi_b(25,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(25+(i-1)*kmax) * sphi_b(25,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(25+(i-1)*kmax) * sphi_b(25,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(25+(i-1)*kmax) * sphi_b(25,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(26+(i-1)*kmax) * sphi_b(26,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(26+(i-1)*kmax) * sphi_b(26,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(26+(i-1)*kmax) * sphi_b(26,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(26+(i-1)*kmax) * sphi_b(26,17+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(27+(i-1)*kmax) * sphi_b(27,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(27+(i-1)*kmax) * sphi_b(27,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(27+(i-1)*kmax) * sphi_b(27,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(28+(i-1)*kmax) * sphi_b(28,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(28+(i-1)*kmax) * sphi_b(28,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(28+(i-1)*kmax) * sphi_b(28,15+s_offset_b1)
buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(29+(i-1)*kmax) * sphi_b(29,1+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(29+(i-1)*kmax) * sphi_b(29,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(29+(i-1)*kmax) * sphi_b(29,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(29+(i-1)*kmax) * sphi_b(29,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(29+(i-1)*kmax) * sphi_b(29,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(30+(i-1)*kmax) * sphi_b(30,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(30+(i-1)*kmax) * sphi_b(30,14+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(30+(i-1)*kmax) * sphi_b(30,18+s_offset_b1)
buffer2(i+imax*(20-1)) = buffer2(i+imax*(20-1)) + buffer1(30+(i-1)*kmax) * sphi_b(30,20+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(31+(i-1)*kmax) * sphi_b(31,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(31+(i-1)*kmax) * sphi_b(31,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(31+(i-1)*kmax) * sphi_b(31,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(31+(i-1)*kmax) * sphi_b(31,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(32+(i-1)*kmax) * sphi_b(32,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(32+(i-1)*kmax) * sphi_b(32,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(32+(i-1)*kmax) * sphi_b(32,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(32+(i-1)*kmax) * sphi_b(32,18+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(33+(i-1)*kmax) * sphi_b(33,5+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(33+(i-1)*kmax) * sphi_b(33,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(34+(i-1)*kmax) * sphi_b(34,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(34+(i-1)*kmax) * sphi_b(34,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(34+(i-1)*kmax) * sphi_b(34,16+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(35+(i-1)*kmax) * sphi_b(35,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(35+(i-1)*kmax) * sphi_b(35,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(36+(i-1)*kmax) * sphi_b(36,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(36+(i-1)*kmax) * sphi_b(36,14+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(37+(i-1)*kmax) * sphi_b(37,21+s_offset_b1)
buffer2(i+imax*(2-1)) = buffer2(i+imax*(2-1)) + buffer1(38+(i-1)*kmax) * sphi_b(38,2+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(38+(i-1)*kmax) * sphi_b(38,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(38+(i-1)*kmax) * sphi_b(38,6+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(38+(i-1)*kmax) * sphi_b(38,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(39+(i-1)*kmax) * sphi_b(39,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(39+(i-1)*kmax) * sphi_b(39,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(39+(i-1)*kmax) * sphi_b(39,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(39+(i-1)*kmax) * sphi_b(39,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(39+(i-1)*kmax) * sphi_b(39,19+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(40+(i-1)*kmax) * sphi_b(40,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(40+(i-1)*kmax) * sphi_b(40,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(40+(i-1)*kmax) * sphi_b(40,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(40+(i-1)*kmax) * sphi_b(40,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(41+(i-1)*kmax) * sphi_b(41,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(41+(i-1)*kmax) * sphi_b(41,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(41+(i-1)*kmax) * sphi_b(41,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(41+(i-1)*kmax) * sphi_b(41,17+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(42+(i-1)*kmax) * sphi_b(42,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(42+(i-1)*kmax) * sphi_b(42,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(42+(i-1)*kmax) * sphi_b(42,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(43+(i-1)*kmax) * sphi_b(43,11+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(43+(i-1)*kmax) * sphi_b(43,15+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(44+(i-1)*kmax) * sphi_b(44,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(44+(i-1)*kmax) * sphi_b(44,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(45+(i-1)*kmax) * sphi_b(45,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(45+(i-1)*kmax) * sphi_b(45,13+s_offset_b1)
buffer2(i+imax*(1-1)) = buffer2(i+imax*(1-1)) + buffer1(46+(i-1)*kmax) * sphi_b(46,1+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(46+(i-1)*kmax) * sphi_b(46,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(46+(i-1)*kmax) * sphi_b(46,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(46+(i-1)*kmax) * sphi_b(46,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(46+(i-1)*kmax) * sphi_b(46,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(47+(i-1)*kmax) * sphi_b(47,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(47+(i-1)*kmax) * sphi_b(47,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(47+(i-1)*kmax) * sphi_b(47,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(47+(i-1)*kmax) * sphi_b(47,18+s_offset_b1)
buffer2(i+imax*(20-1)) = buffer2(i+imax*(20-1)) + buffer1(47+(i-1)*kmax) * sphi_b(47,20+s_offset_b1)
buffer2(i+imax*(3-1)) = buffer2(i+imax*(3-1)) + buffer1(48+(i-1)*kmax) * sphi_b(48,3+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(48+(i-1)*kmax) * sphi_b(48,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(48+(i-1)*kmax) * sphi_b(48,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(48+(i-1)*kmax) * sphi_b(48,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(49+(i-1)*kmax) * sphi_b(49,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(49+(i-1)*kmax) * sphi_b(49,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(49+(i-1)*kmax) * sphi_b(49,16+s_offset_b1)
buffer2(i+imax*(18-1)) = buffer2(i+imax*(18-1)) + buffer1(49+(i-1)*kmax) * sphi_b(49,18+s_offset_b1)
buffer2(i+imax*(5-1)) = buffer2(i+imax*(5-1)) + buffer1(50+(i-1)*kmax) * sphi_b(50,5+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(50+(i-1)*kmax) * sphi_b(50,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(50+(i-1)*kmax) * sphi_b(50,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(51+(i-1)*kmax) * sphi_b(51,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(51+(i-1)*kmax) * sphi_b(51,14+s_offset_b1)
buffer2(i+imax*(16-1)) = buffer2(i+imax*(16-1)) + buffer1(51+(i-1)*kmax) * sphi_b(51,16+s_offset_b1)
buffer2(i+imax*(7-1)) = buffer2(i+imax*(7-1)) + buffer1(52+(i-1)*kmax) * sphi_b(52,7+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(52+(i-1)*kmax) * sphi_b(52,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(53+(i-1)*kmax) * sphi_b(53,12+s_offset_b1)
buffer2(i+imax*(14-1)) = buffer2(i+imax*(14-1)) + buffer1(53+(i-1)*kmax) * sphi_b(53,14+s_offset_b1)
buffer2(i+imax*(9-1)) = buffer2(i+imax*(9-1)) + buffer1(54+(i-1)*kmax) * sphi_b(54,9+s_offset_b1)
buffer2(i+imax*(12-1)) = buffer2(i+imax*(12-1)) + buffer1(55+(i-1)*kmax) * sphi_b(55,12+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,19+s_offset_b1)
buffer2(i+imax*(21-1)) = buffer2(i+imax*(21-1)) + buffer1(56+(i-1)*kmax) * sphi_b(56,21+s_offset_b1)
buffer2(i+imax*(2-1)) = buffer2(i+imax*(2-1)) + buffer1(57+(i-1)*kmax) * sphi_b(57,2+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(57+(i-1)*kmax) * sphi_b(57,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(57+(i-1)*kmax) * sphi_b(57,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(57+(i-1)*kmax) * sphi_b(57,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(57+(i-1)*kmax) * sphi_b(57,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(58+(i-1)*kmax) * sphi_b(58,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(58+(i-1)*kmax) * sphi_b(58,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(58+(i-1)*kmax) * sphi_b(58,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(58+(i-1)*kmax) * sphi_b(58,17+s_offset_b1)
buffer2(i+imax*(19-1)) = buffer2(i+imax*(19-1)) + buffer1(58+(i-1)*kmax) * sphi_b(58,19+s_offset_b1)
buffer2(i+imax*(4-1)) = buffer2(i+imax*(4-1)) + buffer1(59+(i-1)*kmax) * sphi_b(59,4+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(59+(i-1)*kmax) * sphi_b(59,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(59+(i-1)*kmax) * sphi_b(59,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(59+(i-1)*kmax) * sphi_b(59,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(60+(i-1)*kmax) * sphi_b(60,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(60+(i-1)*kmax) * sphi_b(60,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(60+(i-1)*kmax) * sphi_b(60,15+s_offset_b1)
buffer2(i+imax*(17-1)) = buffer2(i+imax*(17-1)) + buffer1(60+(i-1)*kmax) * sphi_b(60,17+s_offset_b1)
buffer2(i+imax*(6-1)) = buffer2(i+imax*(6-1)) + buffer1(61+(i-1)*kmax) * sphi_b(61,6+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(61+(i-1)*kmax) * sphi_b(61,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(61+(i-1)*kmax) * sphi_b(61,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(62+(i-1)*kmax) * sphi_b(62,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(62+(i-1)*kmax) * sphi_b(62,13+s_offset_b1)
buffer2(i+imax*(15-1)) = buffer2(i+imax*(15-1)) + buffer1(62+(i-1)*kmax) * sphi_b(62,15+s_offset_b1)
buffer2(i+imax*(8-1)) = buffer2(i+imax*(8-1)) + buffer1(63+(i-1)*kmax) * sphi_b(63,8+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(63+(i-1)*kmax) * sphi_b(63,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(64+(i-1)*kmax) * sphi_b(64,11+s_offset_b1)
buffer2(i+imax*(13-1)) = buffer2(i+imax*(13-1)) + buffer1(64+(i-1)*kmax) * sphi_b(64,13+s_offset_b1)
buffer2(i+imax*(10-1)) = buffer2(i+imax*(10-1)) + buffer1(65+(i-1)*kmax) * sphi_b(65,10+s_offset_b1)
buffer2(i+imax*(11-1)) = buffer2(i+imax*(11-1)) + buffer1(66+(i-1)*kmax) * sphi_b(66,11+s_offset_b1)
