|
|
Prev | Next | element_prod |
| Matlab or Octave |
= x .* y
|
| C++ |
= element_prod(x, y)
|
= x .* y
(i,j)-th element given by
@[@
z_{i,j} = x_{i,j} * y_{i,j}
@]@
function [ok] = element_prod_ok()
ok = true;
m = 3;
n = 2;
x = rand(m, n);
y = rand(m, n);
% --------------
z = x .* y;
% --------------
[m, n] = size(z);
ok = ok & (m == 3);
ok = ok & (n == 2);
for i = [1 : m]
for j = [1 : n]
zij = x(i, j) * y(i, j);
ok = ok & abs(z(i, j) - zij) < 1e-10;
end
end
return
= element_prod(x, y)
where x, y and z
ublas matrix<double> objects with size
@(@
m \times n
@)@.
# include <mat2cpp.hpp>
bool element_prod_ok(void)
{ bool ok = true;
using namespace mat2cpp;
size_t i, j, m(3), n(2);
matrix<double> x = rand(m, n);
matrix<double> y = rand(m, n);
// -----------------------------------
matrix<double> z = element_prod(x, y);
// -----------------------------------
ok &= (z.size1() == m);
ok &= (z.size2() == n);
for(i = 0; i < m; i++)
{ for(j = 0; j < n; j++)
{ double zij = x(i,j) * y(i,j);
ok &= std::fabs(z(i, j) - zij) < 1e-10;
}
}
return ok;
}