-
Notifications
You must be signed in to change notification settings - Fork 2
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Implement system solver class #86
Conversation
a16d79e
to
a6ab5e0
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
will compile and run next.
When I run the |
b2f2432
to
b22cfcd
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
So both syssolver test with cuda and klu_glu test gives this (for the first matrix) which is obviously wrong
||b-Ax||_2 : 2.88790648555419e-12 (residual norm)
||b-Ax||_2 (CPU) : 2.483452223531182e-12 (residual norm)
||b-Ax||_2/||b||_2 : 1.465835083417962e-18 (scaled residual norm)
||x-x_true||_2 : 5.382658656691419e-12 (solution error)
||x-x_true||_2/||x_true||_2 : 7.898611154542518e-14 (scaled solution error)
||b-Ax_exact||_2 : 1970144.198500361 (control; residual norm with exact solution)
(I fixed this, see the newest commit)
All other tested exampleswork.
The system solver is supposed to be a high-level API for ReSolve. This class should also manage options settings.
This PR delivers main functionality for the
SystemSolver
class. It includes:SystemSolver
class.