/dports/math/freefem++/FreeFem-sources-4.6/src/femlib/ |
H A D | gmres.hpp | 94 Real normb = sqrt((r,r)); in GMRES() local 101 if ( abs(normb) < 1.e-30) in GMRES() 102 normb = 1; in GMRES() 103 if( tol<0) tol = -tol/ normb; // tol < 0 => tol abs.. in GMRES() 104 if ((resid = beta / normb) <= tol) { in GMRES() 105 tol = resid*normb; in GMRES() 134 << abs(s(i+1)) / normb << " < " << tol << endl; in GMRES() 156 << beta / normb << " < " << tol << endl; in GMRES() 157 if ((resid = beta / normb) < tol) { in GMRES() 158 tol = resid*normb; in GMRES() [all …]
|
H A D | CG.cpp | 365 R relerr=1e100 , relres=1e100,normb=0.; in fgmres() local 397 normb = mysnrm2(n,(K*) vi); in fgmres() 412 … cout << " ** fgmres: " << iter << " residus 0 " << abs(g[0]) << " Cl: " << normb << endl ; in fgmres() 413 if (normb < 1.e-20 || eps < 0) normb = 1.; in fgmres() 463 … cout << " fgmres "<< iter << " Res: = " << relres << " Rel res = " << relres/normb << endl; in fgmres() 465 if(relres/normb < abs(eps)) { in fgmres() 469 << "The relative residual is " << relres/normb << endl; in fgmres() 504 cout << " ** restart fgmres iter " <<iter << " res: " << relres/normb << endl; in fgmres() 509 << "The relative residual is " << relres/normb << endl; in fgmres()
|
/dports/math/gsl/gsl-2.7/splinalg/ |
H A D | test.c | 177 double normr, normb; in test_poisson() local 183 normb = gsl_blas_dnrm2(b); in test_poisson() 185 status = (normr <= tol*normb) != 1; in test_poisson() 187 desc, N, normr, normb); in test_poisson() 266 double normr, normb; in test_toeplitz() local 272 normb = gsl_blas_dnrm2(rhs); in test_toeplitz() 274 status = (normr <= tol*normb) != 1; in test_toeplitz() 316 double normr, normb; in test_random() local 322 normb = gsl_blas_dnrm2(b); in test_random() 324 status = (normr <= tol*normb) != 1; in test_random() [all …]
|
/dports/math/xlife++/xlifepp-sources-v2.0.1-2018-05-09/src/solvers/ |
H A D | BicgStabSolver.hpp | 160 real_t normb = b.norm2(); in algorithm() local 161 if ( normb < theBreakdownThreshold ) { normb = 1.; } in algorithm() 170 residue_ = r.norm2() / normb; in algorithm() 201 residue_= s.norm2() / normb; in algorithm() 220 residue_ = r.norm2() / normb; in algorithm() 240 real_t normb = b.norm2(); in algorithm() local 241 if ( normb < theBreakdownThreshold ) { normb = 1.; } in algorithm() 250 residue_ = r.norm2() / normb; in algorithm() 282 residue_= s.norm2() / normb; in algorithm() 302 residue_ = r.norm2() / normb; in algorithm()
|
H A D | CgsSolver.hpp | 152 real_t normb = b.norm2(); in algorithm() local 153 if ( normb< theBreakdownThreshold ) { normb = 1.; } in algorithm() 154 residue_ = r.norm2() / normb; in algorithm() 183 residue_ = r.norm2() / normb; in algorithm() 207 real_t normb = b.norm2(); in algorithm() local 208 if ( normb< theBreakdownThreshold ) { normb = 1.; } in algorithm() 209 residue_ = b.norm2() / normb; in algorithm() 240 residue_ = r.norm2() / normb; in algorithm()
|
/dports/cad/opencascade/opencascade-7.6.0/src/GeomFill/ |
H A D | GeomFill_DraftTrihedron.cxx | 85 b /= normb; in D0() 86 if (normb < 1.e-12) in D0() 125 Standard_Real normT, normb; in D1() local 132 normb = b.Magnitude(); in D1() 133 if (normb < 1.e-12) in D1() 135 b /= normb; in D1() 138 db /= normb; in D1() 193 normb = b.Magnitude(); in D2() 194 if (normb < 1.e-12) in D2() 199 b /= normb; in D2() [all …]
|
/dports/science/lammps/lammps-stable_29Sep2021/lib/atc/ |
H A D | CG.h | 40 Real normb = b.norm(); in CG() local 47 if (normb == 0.0) in CG() 48 normb = 1; in CG() 50 if ((resid = r.norm() / normb) <= tol) { in CG() 74 if ((resid = r.norm() / normb) <= tol) in CG()
|
H A D | GMRES.h | 69 Real normb = p.norm(); in GMRES() local 73 if (normb == 0.0) in GMRES() 74 normb = 1; in GMRES() 76 if ((resid = r.norm() / normb) <= tol) { in GMRES() 105 if ((resid = abs(s(i+1)) / normb) < tol) { in GMRES() 116 if ((resid = beta / normb) < tol) { in GMRES()
|
/dports/science/liggghts/LIGGGHTS-PUBLIC-3.8.0-26-g6e873439/lib/atc/ |
H A D | CG.h | 40 Real normb = b.norm(); in CG() local 47 if (normb == 0.0) in CG() 48 normb = 1; in CG() 50 if ((resid = r.norm() / normb) <= tol) { in CG() 74 if ((resid = r.norm() / normb) <= tol) in CG()
|
H A D | GMRES.h | 69 Real normb = p.norm(); in GMRES() local 73 if (normb == 0.0) in GMRES() 74 normb = 1; in GMRES() 76 if ((resid = r.norm() / normb) <= tol) { in GMRES() 105 if ((resid = abs(s(i+1)) / normb) < tol) { in GMRES() 116 if ((resid = beta / normb) < tol) { in GMRES()
|
/dports/math/fplll/fplll-5.4.1/tests/ |
H A D | test_svp.cpp | 155 FP_NR<mpfr_t> normb; in test_dual_svp() local 156 if (dual_length(normb, A, b)) in test_dual_svp() 186 normb += error; in test_dual_svp() 187 if (norm_sol > normb) in test_dual_svp() 208 FP_NR<mpfr_t> normb; in test_dsvp_reduce() local 209 if (dual_length(normb, A, b)) in test_dsvp_reduce() 245 normb += error; in test_dsvp_reduce() 246 if (norm_sol > normb) in test_dsvp_reduce()
|
H A D | test_svp_gram.cpp | 178 FP_NR<mpfr_t> normb; in test_dual_svp() local 179 if (dual_length(normb, G, b)) in test_dual_svp() 221 normb += error; in test_dual_svp() 222 if (norm_sol > normb) in test_dual_svp() 243 FP_NR<mpfr_t> normb; in test_dsvp_reduce() local 244 if (dual_length(normb, G, b)) in test_dsvp_reduce() 288 normb += error; in test_dsvp_reduce() 289 if (norm_sol > normb) in test_dsvp_reduce()
|
/dports/math/freefem++/FreeFem-sources-4.6/plugin/mpi/ |
H A D | MPICG.cpp | 176 Real normb = sqrt(ReduceSum1((r,r),commworld)); in GMRES_MPI() local 183 if ( abs(normb) < 1.e-30) in GMRES_MPI() 184 normb = 1; in GMRES_MPI() 186 if ((resid = beta / normb) <= tol) { in GMRES_MPI() 215 cout << "GMRES: " << j << " " << abs(s(i+1)) << " " << normb << " " in GMRES_MPI() 216 << abs(s(i+1)) / normb << " < " << tol << endl; in GMRES_MPI() 218 if ((resid = abs(s(i+1)) / normb) < tol) { in GMRES_MPI() 221 << abs(s(i+1)) / normb << " < " << tol << endl; in GMRES_MPI() 239 << beta / normb << " < " << tol << endl; in GMRES_MPI() 240 if ((resid = beta / normb) < tol) { in GMRES_MPI() [all …]
|
/dports/editors/texstudio/texstudio-4.1.2/completion/ |
H A D | strands.cwl | 23 normb= 46 normb= 69 normb= 94 normb=
|
/dports/science/py-scipy/scipy-1.7.1/scipy/sparse/linalg/isolve/ |
H A D | lsmr.py | 237 normb = norm(b) 240 beta = normb.copy() 413 test1 = normr / normb 419 t1 = test1 / (1 + normA * normx / normb) 420 rtol = btol + atol * normA * normx / normb
|
/dports/misc/openmvg/openMVG-2.0/src/third_party/eigen/unsupported/Eigen/src/IterativeSolvers/ |
H A D | IDRS.h | 83 const RealScalar normb = b.norm(); in idrs() local 85 if (internal::isApprox(normb, RealScalar(0))) in idrs() 109 const RealScalar tolb = tol * normb; //Relative tolerance in idrs() 126 relres = normr / normb; in idrs() 243 if (trueres && normr < normb) in idrs() 268 relres=normr/normb; in idrs()
|
/dports/math/gsl/gsl-2.7/multilarge/ |
H A D | normal.c | 36 double normb; /* || b || */ member 198 state->normb = 0.0; in normal_reset() 246 state->normb = gsl_hypot(state->normb, gsl_blas_dnrm2(b)); in normal_accumulate() 469 r2 += state->normb * state->normb; in normal_calc_norms()
|
/dports/graphics/lightzone/LightZone-4.1.8-3-g36e87773/lightcrafts/extsrc/com/lightcrafts/mediax/jai/ |
H A D | IHSColorSpace.java | 452 int normb = 8 - srcComponentSize[2] ; in fromRGBByte() local 575 int normb = 16 - srcComponentSize[2] ; in fromRGBShort() local 701 int normb = 32 - srcComponentSize[2] ; in fromRGBInt() local 1129 double normr = 1.0, normg = 1.0, normb = 1.0 ; in toRGBByte() local 1140 normb = (1l << destComponentSize[2]) - 1 ; in toRGBByte() 1262 double normr = 1.0, normg = 1.0, normb = 1.0 ; in toRGBShort() local 1269 normb = (1l << destComponentSize[2]) - 1 ; in toRGBShort() 1339 double normr = 1.0, normg = 1.0, normb = 1.0 ; in toRGBInt() local 1346 normb = (1l << destComponentSize[2]) - 1 ; in toRGBInt() 1412 double normr = 1.0, normg = 1.0, normb = 1.0 ; in toRGBFloat() local [all …]
|
/dports/science/InsightToolkit/ITK-5.0.1/Modules/ThirdParty/VNL/src/vxl/v3p/netlib/linalg/ |
H A D | lsmrBase.cxx | 52 this->normb = 0.0; in lsmrBase() 350 this->normb = beta; in Solve() 517 test1 = this->normr / this->normb; in Solve() 520 double t1 = test1/(one + this->normA*this->normx/this->normb); in Solve() 521 double rtol = this->btol + this->atol*this->normA*normx/this->normb; in Solve() 597 … << " Exit LSMR. normb = " << this->normb << " ,normx = " << this->normx << std::endl in TerminationPrintOut()
|
H A D | lsmrBase.h | 503 double normb; variable
|
/dports/misc/vxl/vxl-3.3.2/v3p/netlib/linalg/ |
H A D | lsmrBase.cxx | 52 this->normb = 0.0; in lsmrBase() 350 this->normb = beta; in Solve() 517 test1 = this->normr / this->normb; in Solve() 520 double t1 = test1/(one + this->normA*this->normx/this->normb); in Solve() 521 double rtol = this->btol + this->atol*this->normA*normx/this->normb; in Solve() 597 … << " Exit LSMR. normb = " << this->normb << " ,normx = " << this->normx << std::endl in TerminationPrintOut()
|
H A D | lsmrBase.h | 503 double normb; variable
|
/dports/math/adol-c/ADOL-C-2.7.2/ADOL-C/examples/ |
H A D | luexam.cpp | 204 double *mat, *rhs, *ans, err, normx, normb; in main() local 246 normb = norm2(rhs, n); in main() 248 cout << "Norm rhs = " << normb <<"\n"; in main() 272 normb = norm2(rhs, n); in main() 274 cout << "Norm rhs = " << normb <<"\n"; in main()
|
/dports/math/dune-uggrid/dune-uggrid-bc2d1229420367563410ce9e519f5ff82b45266f/dune/uggrid/gm/ |
H A D | evm.cc | 384 DOUBLE normb, scprd; in V3_Project() local 386 normb = b[0]*b[0]+b[1]*b[1]+b[2]*b[2]; in V3_Project() 387 if (normb==0.0) in V3_Project() 392 scprd /= normb; in V3_Project()
|
/dports/math/freefem++/FreeFem-sources-4.6/idp/ |
H A D | ffddm_functions.idp | 579 real normb = 0; 587 normb = sqrt(real(pr#prfe#scalprod(wi,wi))); 590 normb = sqrt(real(pr#prfe#scalprod(vi,vi))); 592 if (normb < 1.e-20 || eps < 0) normb = 1.; 619 if(abs(g[0])/normb < abs(eps)) { 622 … cout << "["+Stringification(pr)+"] The relative residual is " + abs(g[0])/normb << endl; 720 if(relres/normb < abs(eps)) { 723 cout << "["+Stringification(pr)+"] The relative residual is " + relres/normb << endl; 731 if(relres/normb < abs(eps) || pr#nbiter > maxit) 738 if(relres/normb >= abs(eps)) [all …]
|