Debug asseption failed! is_block_type_valid(header->_block_use)

#include <iostream>
#include <cstdlib>
#include "array.h"
using namespace std;
 
class SimplexMethod {
    int row, col, vedrow = 0, vedcol = 0;
    int * minarr;
    double ** arr;
    double * farr;
    double f = 0.0;
public:
    SimplexMethod();
    ~SimplexMethod();
    double vedelement();
    void iteration();
    void showmatrix();
    void reshenie();
};
 
double SimplexMethod::vedelement()
{
    int i, j, maxel = -1;
    for (i = 0; i<col-1; ++i)
    {
        minarr[i] = -1;
        for (j = 0; j < row; ++j)
            if (arr[j][i] > 0)
                if (minarr[i] == -1)
                    minarr[i] = j;
                else if (arr[j][col - 1] / arr[j][i] * farr[i]<arr[minarr[i]][col - 1] / arr[minarr[i]][i] * farr[i])
                    minarr[i] = j;
    }
    for (i = 0; i<col-1; ++i)
    {
 
        if (farr[i] > 0)
            if (maxel == -1)
                maxel = i;
            else if (arr[minarr[i]][col - 1] / arr[minarr[i]][i] * farr[i]>arr[minarr[maxel]][col - 1] / arr[minarr[maxel]][maxel] * farr[maxel])
                maxel = i;
    }
    vedrow = minarr[maxel];
    vedcol = maxel;
    cout<<"ved. element= "<< arr[minarr[maxel]][maxel]<<" row= "<< minarr[maxel]+1 <<" col= "<<maxel+1<<endl;
    return arr[minarr[maxel]][maxel];
}
 
void SimplexMethod::iteration()
{
    int i, j;
    double temp;
    temp = arr[vedrow][vedcol];
    for (i = 0; i<col; ++i)
        arr[vedrow][i] /= temp;
    f -= arr[vedrow][col - 1] * farr[vedcol];
    for (i = 0; i<row; ++i)
    {
        if (i == vedrow)
            continue;
        temp = arr[i][vedcol];
        for (j = 0; j<col; ++j)
            arr[i][j] -= temp*arr[vedrow][j];
    }
    temp = farr[vedcol];
    for (i = 0; i<col - 1; ++i)
        farr[i] -= temp*arr[vedrow][i];
}
 
void SimplexMethod::reshenie()
{
    int maxiter, j, k=1;
    cout << "input maxiteration" << endl;
    cin >> maxiter;
    while (positivv(farr, col-1) && maxiter)
    {
        if (vedelement() != 0.0)
        {
            iteration();
            cout << k++ << " iteration" << endl;
            showmatrix();
        }
        else
            break;
        --maxiter;
    }
    for (int i = 0; i<col - 1; ++i)
    {
        cout << "x[" << i + 1 << "]= ";
        if (farr[i]<0)
            cout << 0 << endl;
        else
        {
            for (j = row - 1; j>-1; --j)
                if (arr[j][i] == 1)
                {
                    cout << arr[j][col - 1] << endl;;
                    break;
                }
            if (j == -1)
                cout << 0 << endl;
        }
    }
    cout << "max F= " << f*(-1) << endl;
}
 
SimplexMethod::SimplexMethod() {
    cout << "input kol row for coefficient matrix" << endl;
    cin >> row;
    cout << "input kol x for coefficient matrix" << endl;;
    cin >> col;
    ++col;
    arr = newarr(row, col);
    farr = newv(col - 1);
    minarr = newiv(col - 1);
    showmatrix();
}
 
SimplexMethod::~SimplexMethod() {
    for (int i = 0; i<row; ++i)
        free(*(arr+i));
    free(arr);
    free(farr);
    free(minarr);
}
 
void SimplexMethod::showmatrix()
{
    specprintarr(arr, farr, f, row, col);
}
 
int main()
{
    {
        SimplexMethod m1;
        m1.reshenie();
    }
    system("pause");
    return 0;
}

#include <cmath>
#include "array.h"
#include <iostream>
using namespace std;
double** newarr(int row, int col)
{
    double **rowarr = (double**)malloc(row * sizeof *rowarr);
    double *fullarr = (double*)malloc(row*col * sizeof(double));
    for (int i = 0, j; i<row; ++i)
    {
        rowarr[i] = fullarr + i*col;
        for (j = 0; j<col; ++j)
        {
            printf("input a[%i][%i]\n", i, j);
            cin >> (*(*(rowarr + i) + j));
        }
    }
 
    return rowarr;
}
 
double* newv(int row)
{
    double* varr = (double*)malloc(row * sizeof(double));
    for (int i = 0; i<row; ++i)
    {
        printf("input a[%i]\n", i);
        cin >> (*(varr + i));
    }
    return varr;
}
 
int* newiv(int row)
{
    return (int*)malloc(row * sizeof(int));
}
int positivv(double *far, int col)
{
    for (--col; col>-1; --col)
        if (far[col]>0)
            return 1;
    return 0;
}
void specprintarr(double **ar, double *far, double f, int row, int col)
{
    int i, j;
    for (i = 0; i<row; ++i)
        for (j = 0; j<col; ++j)
            printf("%8.5f%c", ar[i][j], (j == col - 1) ? '\n' : ' ');
    for (j = 0; j<col - 1; ++j)
        printf("%8.5f ", far[j]);
    printf("%8.5f\n", f);
 
}
 
void printarr(double **arr, int row, int col)
{
    for (int i = 0, j; i<row; ++i)
        for (j = 0; j<col; ++j)
            printf("a[%i][%i]= %.5f\n", i, j, arr[i][j]);
}
 
void printv(double *arr, int row)
{
    for (int i = 0; i<row; ++i)
        printf("x[%i]= %.5f\n", i, arr[i]);
}