|
Ublas : |
Subject: Re: [ublas] todo list
From: Jesse Perla (jesseperla_at_[hidden])
Date: 2010-07-24 17:14:32
Matwey V. Kornilov <matwey.kornilov <at> gmail.com> writes:
> This is for memory economy. When we want to solve linear set Ax=y by means
> of SVD decomposition we write
Let me throw out a crazy suggestion. One cool thing about matlab is
overloading of the '\' operator. You can then solve: A x = B
by writing:
x = B\A;
Can we do something like this overloading the '/' operator?
matrix<double> A;
vector<double> B;
auto x = B / A;
And to solve in place:
B /= A;
To specify different algorithms, you could add in tags around the rhs, which
could then be detected in the '/' overload and the appropriate algorithm used.
auto x = B / ublas::svd(A);