Баг в коде Reaction Diffusion Gray-Scott model

#include <gl/glut.h>
#include <iostream> 
using namespace std;
 
const int imSiz=300; // размер картинки
 
float grid[imSiz+2][imSiz+2][2]={0};  // основной массив
float next[imSiz+2][imSiz+2][2]={0}; // буферный массив
bool render_complete=0;
float
Da=1.0,
Db=0.5,
k=0.062,
f=0.085;
 
 
//==============================================
//==============================================
//рисуем пиксели с окрасом
void drawPixelArray() 
{
    glClear(GL_COLOR_BUFFER_BIT);
    glBegin(GL_POINTS);
    for(int i = 1; i < imSiz+1; i++)
    {
        for (int j = 1; j < imSiz+1; j++)
        {
            float a=grid[i][j][1];
            glColor3f(a, a, a);
            glVertex2i(i,j);
        }
    }
    glEnd();
    render_complete=1;
}
//рисуем пиксели с окрасом
//==============================================
//==============================================
 
 
 
 
//==============================================
//==============================================
//====расчет  новой концетрации 
void conectratCorrection()
{
 
//------------------------------------------
//очистка буфера next
for(int i = 1; i < imSiz+1; i++)
{
    for (int j = 1; j < imSiz+1; j++)
    {
    next[i][j][0]=0;
    next[i][j][1]=0;
    }
}
//очистка буфера next
//------------------------------------------
 
 
 
//------------------------------------------
//свертка массива greed с ядром 3 на 3, лапласиан    
for(int i = 1; i < imSiz+1; i++)
{
    for (int j = 1; j < imSiz+1; j++)
    {
    float sumA=
    grid[i][j][0]*-1.0+
    grid[i-1][j][0]*0.2+
    grid[i+1][j][0]*0.2+
    grid[i][j+1][0]*0.2+
    grid[i][j-1][0]*0.2+
    grid[i-1][j-1][0]*0.05+
    grid[i+1][j-1][0]*0.05+
    grid[i+1][j+1][0]*0.05+
    grid[i-1][j+1][0]*0.05;
 
    float sumB=
    grid[i][j][1]*-1.0+
    grid[i-1][j][1]*0.2+
    grid[i+1][j][1]*0.2+
    grid[i][j+1][1]*0.2+
    grid[i][j-1][1]*0.2+
    grid[i-1][j-1][1]*0.05+
    grid[i+1][j-1][1]*0.05+
    grid[i+1][j+1][1]*0.05+
    grid[i-1][j+1][1]*0.05; 
 
    float a=grid[i][j][0];
    float b=grid[i][j][1];
    //next[i][j][0]=a+(Da*sumA)-(a*b*b)+(f*(1.0-a));
    //next[i][j][1]=b+(Db*sumB)+(a*b*b)-((k+f)*b);
    next[i][j][0]=(Da*sumA)-(a*b*b)+(f*(1.0-a));
    next[i][j][1]=(Db*sumB)+(a*b*b)-((k+f)*b);
    }
}
//свертка массива greed с ядром 3 на 3, лапласиан
//------------------------------------------
 
 
//------------------------------------------
//сумма старых значений  grid и новых next.
for(int i = 1; i < imSiz+1; i++)
{
    for (int j = 1; j < imSiz+1; j++)
    {
        grid[i][j][0]+=next[i][j][0];
        grid[i][j][1]+=next[i][j][1];
    }
}
 
//сумма старых значений  grid и новых next.
//------------------------------------------
drawPixelArray(); 
}
//====расчет  новой концетрации
//==============================================
//==============================================
 
 
//==============================================
//==============================================
//нарисовать квадрат 10 пикселей веществом B
void genArr()
{
for(int i = 0; i < imSiz+2; i++)
{
    for (int j = 0; j < imSiz+2; j++)
    {
        if(i>120 && i<180 && j>120 && j<180)
        {
            grid[i][j][0]=1;
            grid[i][j][1]=1;
        }
        else
        {
            grid[i][j][0]=1;
            grid[i][j][1]=0;
        }
    }
}
//drawPixelArray();
conectratCorrection();
}
//нарисовать квадрат 10 пикселей веществом B
//==============================================
//==============================================
 
 
 
//==============================================
//==============================================
//проверяет готов ли рендер, если да то обновляет экран и запуск таймера
void TimerFunction(int value)
{
    if(render_complete) 
    {
        glutSwapBuffers();
        render_complete=0;
        glutTimerFunc(33, TimerFunction, 1);
        cout << "New frame"<<endl;
        conectratCorrection();
    }
    else
    {
        cout << "Frame lost!----------"<<endl;
    }
}
//проверяет готов ли рендер, если да то обновляет экран и запуск таймера
//==============================================
//==============================================
 
 
void setupGL() 
{
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    glOrtho(0, imSiz , imSiz, 0, -1, 1);
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
    glClearColor(0, 0, 0, 0);
    glPointSize(1);
    genArr();
}
 
int main(int argc, char** argv) {
    glutInit(&argc, argv);
    glutInitDisplayMode(GLUT_DOUBLE | GLUT_RGB);
    glutInitWindowSize(imSiz, imSiz); 
    glutInitWindowPosition(500, 200);
    glutCreateWindow("Reaction Diffusion Gray-Scott model");
    setupGL();
    glutTimerFunc(33, TimerFunction, 1);
    glutMainLoop();
    return 0;
}

float smallConcentration=0.7f;  
float setConcentrationUp=1.0f; // концентрации растворов
float setConcentrationDn=0.7f; //