49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
AA (row,col) = 0.0;
AA (col,row) = 0.0;
}
}
}
#ifdef OCTAVE_DEV
x = static_cast<ColumnVector> (xleftdiv (AA, static_cast<Matrix>(AAx + dx * t)));
#else // !OCTAVE_DEV
Matrix tmp = xleftdiv (AA, static_cast<Matrix>(static_cast<ColumnVector>(AAx + dx * t)));
for (row = 0; row < Nx; row++)
{
x (row) = tmp (row,0);
}
#endif // OCTAVE_DEV
|
|
|
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
AA (row,col) = 0.0;
AA (col,row) = 0.0;
}
}
}
#ifdef OCTAVE_DEV
x = xleftdiv (AA, Matrix(AAx + dx * t)).column (0);
#else // !OCTAVE_DEV
Matrix tmp = xleftdiv (AA, static_cast<Matrix>(static_cast<ColumnVector>(AAx + dx * t)));
for (row = 0; row < Nx; row++)
{
x (row) = tmp (row,0);
}
#endif // OCTAVE_DEV
|