__host__ void solveMatrixGaussCPU(float* matrix, int width, int height)
{
    float coeff;
    for (int i = 0; i < height; i++)
            for (int k = i; k + 1 < height; k++)
            {
                coeff = matrix[i * width + i] / matrix[(k + 1) * width + i];
                for (int j = 0; j < width; j++)
                    matrix[(k + 1) * width + j] = matrix[(k + 1) * width + j] - (matrix[i * width + j] / coeff);
            }
}