|
|
Prev | Next |
| Matlab or Octave |
= x \ y
|
| C++ |
= matrix_div(x, y, rank)
|
= x \ y
function [ok] = matrix_div_ok()
ok = true;
m = 3;
n = 2;
x = rand(m, m);
y = rand(m, n);
% -------------
z = x \ y;
% -------------
[m, n] = size(z);
ok = ok & (m == 3);
ok = ok & (n == 2);
ok = ok & all ( all( abs(x * z - y) < 1e-10 ) );
return
= matrix_div(x, y, rank)
where x, y and z are
ublas matrix<double> objects.
The size of x is @(@
m \times m
@)@,
the size of y is @(@
m \times n
@)@,
and the size of z is @(@
m \times n
@)@ matrix.
The return value of the size_t argument rank is the rank
of the matrix x.
# include <boost/numeric/ublas/matrix.hpp>
# include <mat2cpp.hpp>
bool matrix_div_ok(void)
{ bool ok = true;
using namespace mat2cpp;
size_t i, j, m(3), n(2);
matrix<double> x = rand(m, m);
matrix<double> y = rand(m, n);
// ---------------------------------------
size_t rank;
matrix<double> z = matrix_div(x, y, rank);
// ---------------------------------------
ok &= (rank == m);
ok &= (z.size1() == m);
ok &= (z.size2() == n);
matrix<double> x_z = prod(x, z);
for(i = 0; i < m; i++)
{ for(j = 0; j < n; j++)
ok &= std::fabs( x_z(i,j) - y(i,j) ) < 1e-10;
}
return ok;
}