Home
last modified time | relevance | path

Searched refs:normb (Results 1 – 25 of 69) sorted by relevance

123

/dports/math/freefem++/FreeFem-sources-4.6/src/femlib/
H A Dgmres.hpp94 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 DCG.cpp365 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 Dtest.c177 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 DBicgStabSolver.hpp160 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 DCgsSolver.hpp152 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 DGeomFill_DraftTrihedron.cxx85 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 DCG.h40 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 DGMRES.h69 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 DCG.h40 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 DGMRES.h69 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 Dtest_svp.cpp155 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 Dtest_svp_gram.cpp178 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 DMPICG.cpp176 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 Dstrands.cwl23 normb=
46 normb=
69 normb=
94 normb=
/dports/science/py-scipy/scipy-1.7.1/scipy/sparse/linalg/isolve/
H A Dlsmr.py237 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 DIDRS.h83 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 Dnormal.c36 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 DIHSColorSpace.java452 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 DlsmrBase.cxx52 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 DlsmrBase.h503 double normb; variable
/dports/misc/vxl/vxl-3.3.2/v3p/netlib/linalg/
H A DlsmrBase.cxx52 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 DlsmrBase.h503 double normb; variable
/dports/math/adol-c/ADOL-C-2.7.2/ADOL-C/examples/
H A Dluexam.cpp204 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 Devm.cc384 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 Dffddm_functions.idp579 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 …]

123