Changeset 3993889 in tol
- Timestamp:
- Feb 4, 2009, 7:55:16 PM (18 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- f4016fe
- Parents:
- 46184cc
- File:
-
- 1 edited
-
tol/stdlib/math/linalg/_solve_precond_sym.tol (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
tol/stdlib/math/linalg/_solve_precond_sym.tol
r46184cc r3993889 128 128 { 129 129 //As operator is symmetric parameter traspose is noy used in this case 130 CholeskiSolve( SPDP::_.Lc,_.S*u,"PtLLtP")130 CholeskiSolve(_.Lc,_.S*u,"PtLLtP") 131 131 }; 132 //Solves linear system M*xy 132 133 //Solves linear system M*xy 133 134 NameBlock Solve(VMatrix y, Real tolerance, Real forze.lanczos) … … 136 137 VMatrix _.y = y; 137 138 Real _.tol.solver = tolerance; 138 Real _.tol.mrs = (tolerance ^2)/_.m;139 VMatrix _.v = If(!_.do.normalize, _.y, _.Di*_.y);140 VMatrix _.u0 = CholeskiSolve( SPDP::_.Lc,_.v,"PtLLtP");139 Real _.tol.mrs = (tolerance)^2/_.m; 140 VMatrix _.v = If(!_.do.normalize,y, _.Di*y); 141 VMatrix _.u0 = CholeskiSolve(_.Lc,_.v,"PtLLtP"); 141 142 VMatrix _.u = If(And(!forze.lanczos, !_.c), _.u0, 142 143 {
Note: See TracChangeset
for help on using the changeset viewer.
