[РЕШЕНО] Ошибка сегментирования при вызове malloc()

        nv = 17
        n = 6*nv - 4;
	nb = 3*nv - 2;
	nr = nb - 1;
	n_ef = 2*nv - 2;

	// Set mass
	m = new double[nb];

	memset(m, 0, nb*sizeof(double));

	for (int i = 0; i < nv; i++)
	{
		int k = 3*i;

		m[k] = mv;
		m[k+1] = mt;
		m[k+2] = mt;
	}

	// Create matrixes

	A = create_matrix(nb, nb+nr);

	L = create_matrix(nr, nb);
	L1 = create_matrix(nr, nb);

	for (int i = 0; i < nr; i++)
	{
		L[i][i] = 1;
		L[i][i+1] = -1;
	}

	K = create_matrix(nb, nr);
	K1 = create_matrix(nb, nr);

	for (int i = 0; i < nv-1; i++)
	{
		int j = 3*i;

		K1[j][j] = K[j][j] = -1;
		K1[j+1][j] = K[j+1][j] = 1;

		K1[j+1][j+1] = K[j+1][j+1] = 1;
		K1[j+2][j+1] = K[j+2][j+1] = -1;

		K1[j+2][j+2] = K[j+2][j+2] = -1;
		K1[j+3][j+2] = K[j+3][j+2] = 1;
	}
	lock = new bool[nr];

	for (int i = 0; i < nr; i++)
		lock[i] = true;

	p = new double[nb];
	q = new double[nb];

	memset(p, 0, nb*sizeof(double));
	memset(q, 0, nb*sizeof(double));

	Tef = new double[n_ef];
	ef_sign = new int[n_ef];

	n_ef = 32;
	coup = new coup_t[n_ef];
	for (int i = 0; i < n_ef; i++)
	{
		coup_t tmp = new CEFcoupling();
		coup[i] = tmp;
	}

coup_t tmp = new CEFcoupling();

Program received signal SIGSEGV, Segmentation fault.
0x00007ffff707668f in _int_malloc () from /usr/lib/libc.so.6
(gdb) bt full
#0  0x00007ffff707668f in _int_malloc () from /usr/lib/libc.so.6
No symbol table info available.
#1  0x00007ffff7078030 in malloc () from /usr/lib/libc.so.6
No symbol table info available.
#2  0x00007ffff7919e9d in operator new(unsigned long) () from /usr/lib/libstdc++.so.6
No symbol table info available.
#3  0x00007ffff7919f99 in operator new[](unsigned long) () from /usr/lib/libstdc++.so.6
No symbol table info available.
#4  0x0000000000401fe9 in ode_systems_init () at train.cpp:177
No locals.
#5  0x0000000000400b69 in main () at main.cpp:23
        solver = 0x400a50 <_start>

//-------------------------------------------------------------------
//
//-------------------------------------------------------------------
matrix_t create_matrix(int n, int m)
{
	matrix_t tmp = 0;

	tmp = new matrix_str_t[n];

	for (int i = 0; i < n; i++)
	{
		tmp[i] = new double[m];
		memset(tmp[i], 0, m*sizeof(double));
	}

	return tmp;
}

typedef double* matrix_str_t;
typedef matrix_str_t* matrix_t;

for (int i = 0; i < nr; i++)
	{
		L[i][i] = 1;
		L[i][i+1] = -1;
	}

maisvendoo
Имеется участок кода

        nv = 17
        n = 6*nv - 4;
	nb = 3*nv - 2;
	nr = nb - 1;
	n_ef = 2*nv - 2;

	// Set mass
	m = new double[nb];

	memset(m, 0, nb*sizeof(double));

	for (int i = 0; i < nv; i++)
	{
		int k = 3*i;

		m[k] = mv;
		m[k+1] = mt;
		m[k+2] = mt;
	}

        m[0] = mv;
	m[1] = mt;

	for (int i = 2; i < nb-2; i += 3)
	{
		m[i] = mt;
		m[i+1] = mv;
		m[i+2] = mt;
	}

	m[nb-2] = mt;
	m[nb-1] = mv;

cucullus
оси не подписаны!

cucullus
А что есть отрицательная тяга?