void DifferenceScheme::solveOMP()
{
    n = round((xMax - xMin) / xStep);
    m = round((tMax - tMin) / tStep);
    u = dMatrix(n + 1, dVector(m + 1, initValue));
    TridiagonalMatrix tm(n + 1);
    dVector F(n + 1);
    Progonka progonka(tm, F);
    double error;

    workTime = clock();
#pragma omp parallel for
    for (int i = 0; i <= n; i++)
    {
        u[i][0] = w(xMin + i * xStep, tMin);
    }
#pragma omp parallel for
    for (int i = 0; i <= m; i++)
    {
        u[0][i] = w(xMin, tMin + i * tStep);
        u[n][i] = w(xMax, tMin + i * tStep);
    }

    for (int j = 1; j <= m; j++)
    {
        do
        {
#pragma omp parallel for
            for (int i = 1; i < n; i++)
            {
                tm.lower[i - 1] = dLeft(i, j);
                tm.main[i] = dCenter(i, j);
                tm.upper[i] = dRight(i, j);
                F[i] = -1 * f(i, j);
            }
            tm.main[0] = tm.main[n] = 1;
            F[0] = -1 * (u[0][j] - w(xMin, tMin + j * tStep));
            F[n] = -1 * (u[n][j] - w(xMax, tMin + j * tStep));
            progonka.solveOMP();
#pragma omp parallel for
            for (int i = 1; i < n; i++)
            {
                u[i][j] += progonka.x[i];
            }
            error = 0;
            for (int i = 1; i < n; i++)
                if (fabs(progonka.x[i]) > error)
                    error = fabs(progonka.x[i]);
        }
        while (error > accuracy);
    }
    workTime = clock() - workTime;
}