/dports/misc/visp/visp-3.4.0/modules/visual_features/src/feature-builder/ |
H A D | vpFeatureBuilderVanishingPoint.cpp | 127 double rho_1 = L1.getRho(); in create() local 134 double denom = sqrt( rho_1*rho_1 + rho_2*rho_2 - 2*rho_1*rho_2 * cos(theta_diff) ); in create() 136 …double alpha = atan2( rho_1 * cos(theta_2) - rho_2 * cos(theta_1), rho_2*sin(theta_1) - rho_1 * si… in create() 202 double rho_1 = ( (x1+x2)*cos(theta_1) + (y1+y2)*sin(theta_1) ) / 2.; in create() local 219 double x = (sin(theta_1) * rho_2 - sin(theta_2)*rho_1 ) / sin(theta_diff); in create() 220 double y = (cos(theta_2) * rho_1 - cos(theta_1)*rho_2 ) / sin(theta_diff); in create() 227 double denom = sqrt( rho_1*rho_1 + rho_2*rho_2 - 2*rho_1*rho_2 * cos(theta_diff) ); in create() 232 …double alpha = atan2( rho_1 * cos(theta_2) - rho_2 * cos(theta_1), rho_2*sin(theta_1) - rho_1 * si… in create() 239 double denom = sqrt( rho_1*rho_1 + rho_2*rho_2 - 2*rho_1*rho_2 * cos(theta_diff) ); in create() 240 …double alpha = atan2( rho_1 * cos(theta_2) - rho_2 * cos(theta_1), rho_2*sin(theta_1) - rho_1 * si… in create()
|
/dports/math/xlife++/xlifepp-sources-v2.0.1-2018-05-09/src/solvers/ |
H A D | CgsSolver.hpp | 149 K rho_1(1), rho_2(1), alpha(1), beta(1); in algorithm() local 165 assign(rho_1, dotRC(rtilde, r)); in algorithm() 166 if (std::abs(rho_1) < theBreakdownThreshold) breakdown(std::abs(rho_1),"rho_1"); in algorithm() 171 beta = rho_1 / rho_2; in algorithm() 176 assign(alpha, rho_1 / dotRC(rtilde, v)); in algorithm() 182 rho_2 = rho_1; in algorithm() 204 K rho_1(1), rho_2(1), alpha(1), beta(1); in algorithm() local 220 assign(rho_1, dotRC(rtilde, r)); in algorithm() 221 if (std::abs(rho_1) < theBreakdownThreshold) breakdown(std::abs(rho_1),"rho_1"); in algorithm() 226 beta = rho_1 / rho_2; in algorithm() [all …]
|
H A D | BicgStabSolver.hpp | 164 K rho_1, rho_2(0), alpha(0), beta, omega(0), sigma; in algorithm() local 175 assign(rho_1, hermitianProduct(r,rtilde)); in algorithm() 176 if (std::abs(rho_1) < theBreakdownThreshold) { breakdown(std::abs(rho_1), "rho_1"); } in algorithm() 183 beta = (rho_1 / rho_2) * (alpha / omega); in algorithm() 195 alpha = rho_1 / sigma; in algorithm() 221 rho_2 = rho_1; in algorithm() 255 assign(rho_1, hermitianProduct(r,rtilde)); in algorithm() 256 if (std::abs(rho_1) < theBreakdownThreshold) { breakdown(std::abs(rho_1), "rho_1"); } in algorithm() 263 beta = (rho_1 / rho_2) * (alpha / omega); in algorithm() 276 alpha = rho_1 / sigma; in algorithm() [all …]
|
/dports/science/quantum-espresso/q-e-qe-6.7.0/TDDFPT/src/ |
H A D | lr_calc_dens.f90 | 95 rho_1(:,:) = 0.0d0 108 IF ( doublegrid ) CALL fft_interpolate(dffts, rho_1(:,1), dfftp, rho_1(:,1)) 111 CALL mp_sum(rho_1, inter_bgrp_comm) 152 rho_1(:,is) = rho_1(:,is) + DBLE(psic(:)) 169 CALL mp_sum(rho_1, inter_pool_comm) 182 rho_sum = SUM(rho_1(:,is)) 232 rho_sum_resp_x(i)=rho_sum_resp_x(i)+rho_1(ir,1) 233 rho_sum_resp_y(j)=rho_sum_resp_y(j)+rho_1(ir,1) 234 rho_sum_resp_z(k)=rho_sum_resp_z(k)+rho_1(ir,1) 450 rho_1(ir,1) = rho_1(ir,1) & [all …]
|
/dports/math/gmm++/gmm-5.4/include/gmm/ |
H A D | gmm_solver_bicgstab.h | 94 T rho_1, rho_2(0), alpha(0), beta, omega(0); in bicgstab() local 108 rho_1 = gmm::vect_sp(rtilde, r); in bicgstab() 109 if (rho_1 == T(0)) { in bicgstab() 124 beta = (rho_1 / rho_2) * (alpha / omega); in bicgstab() 131 alpha = rho_1 / gmm::vect_sp(v, rtilde); in bicgstab() 145 rho_2 = rho_1; in bicgstab()
|
H A D | gmm_solver_cg.h | 92 T rho, rho_1(0), a; in cg() 110 add(z, scaled(p, rho / rho_1), p); in cg() 117 rho_1 = rho; in cg() 132 T rho, rho_1(0), a; in cg() 147 add(r, scaled(p, rho / rho_1), p); in cg() 153 rho_1 = rho; in cg()
|
H A D | gmm_solver_constrained_cg.h | 62 value_type rho, rho_1, alpha; in pseudo_inverse() local 79 rho_1 = rho; in pseudo_inverse() 81 add(r, scaled(p, rho / rho_1), p); in pseudo_inverse() 105 value_type rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 136 rho_1 = rho; rho = vect_sp(PS, r, z); // ... in constrained_cg() 142 else gamma = std::max(0.0, (rho - vect_sp(PS, old_z, z) ) / rho_1); in constrained_cg()
|
/dports/cad/gmsh/gmsh-4.9.2-source/contrib/gmm/ |
H A D | gmm_solver_bicgstab.h | 94 T rho_1, rho_2(0), alpha(0), beta, omega(0); in bicgstab() local 108 rho_1 = gmm::vect_sp(rtilde, r); in bicgstab() 109 if (rho_1 == T(0)) { in bicgstab() 124 beta = (rho_1 / rho_2) * (alpha / omega); in bicgstab() 131 alpha = rho_1 / gmm::vect_sp(v, rtilde); in bicgstab() 145 rho_2 = rho_1; in bicgstab()
|
H A D | gmm_solver_cg.h | 92 T rho, rho_1(0), a; in cg() 114 add(z, scaled(p, rho / rho_1), p); in cg() 128 rho_1 = rho; in cg() 147 T rho, rho_1(0), a; in cg() 166 add(r, scaled(p, rho / rho_1), p); in cg() 180 rho_1 = rho; in cg()
|
H A D | gmm_solver_constrained_cg.h | 62 value_type rho, rho_1, alpha; in pseudo_inverse() local 79 rho_1 = rho; in pseudo_inverse() 81 add(r, scaled(p, rho / rho_1), p); in pseudo_inverse() 106 value_type rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 137 rho_1 = rho; rho = vect_sp(PS, r, z); // ... in constrained_cg() 143 else gamma = std::max(0.0, (rho - vect_sp(PS, old_z, z) ) / rho_1); in constrained_cg()
|
/dports/science/mbdyn/mbdyn-1.7.3/mbdyn/base/ |
H A D | bicg.cc | 111 doublereal rho_1; in Solve() local 226 rho_1 = dErr*dErr; /*rhat.InnerProd(r); */ in Solve() 264 std::cerr << "rho_1 " << rho_1 << std::endl; in Solve() 272 rho_1 = rHat.InnerProd(*pr); in Solve() 275 std::cerr << "rho_1 " << rho_1 << std::endl; in Solve() 278 if (fabs(rho_1) < std::numeric_limits<doublereal>::epsilon()) { in Solve() 285 beta = (rho_1/rho_2) * (alpha/omega); in Solve() 308 alpha = rho_1 / alpha; in Solve() 340 rho_2 = rho_1; in Solve()
|
/dports/science/dynare/dynare-4.6.4/matlab/ |
H A D | bicgstab_.m | 30 rho_1 = rh_0'*r; variable 35 beta = (rho_1/rho_0)*(alpha/w); 38 alpha = rho_1/(rh_0'*v); 42 rho_0 = rho_1; 43 rho_1 = -w*(rh_0'*t); variable
|
/dports/science/dynare/dynare-4.6.4/tests/deterministic_simulations/multiple_lead_lags/ |
H A D | AR2.mod | 3 parameters rho_1 rho_2; 5 rho_1=0.2; 10 y_backward=rho_1*y_backward(-1)+rho_2*y_backward(-2); 44 …junk(ii)=M_.params(strmatch('rho_1',M_.param_names,'exact'))*junk(ii-1)+M_.params(strmatch('rho_2'…
|
H A D | AR2_forward.mod | 3 parameters rho_1 rho_2; 5 rho_1=0.2; 10 y_backward=rho_1*y_backward(-1)+rho_2*y_backward(-2); 46 …junk(ii)=M_.params(strmatch('rho_1',M_.param_names,'exact'))*junk(ii-1)+M_.params(strmatch('rho_2'…
|
/dports/math/eigen3/eigen-3.3.9/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/math/stanmath/math-4.2.0/lib/eigen_3.3.9/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/math/R-cran-RcppEigen/RcppEigen/inst/include/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/devel/upp/upp/uppsrc/plugin/Eigen/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/math/libmesh/libmesh-1.6.2/contrib/eigen/3.2.9/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/science/pcmsolver/pcmsolver-1.3.0/external/eigen3/include/eigen3/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/math/libsemigroups/libsemigroups-1.3.7/extern/eigen-3.3.7/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/graphics/wdune/wdune-1.926/vcglib/eigenlib/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/misc/opennn/opennn-5.0.5/eigen/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/misc/openmvg/openMVG-2.0/src/third_party/eigen/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 162 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|
/dports/devel/taskflow/taskflow-3.2.0/3rd-party/eigen-3.3.7/unsupported/Eigen/src/IterativeSolvers/ |
H A D | ConstrainedConjGrad.h | 58 Scalar rho, rho_1, alpha; in pseudo_inverse() local 80 rho_1 = rho; in pseudo_inverse() 82 p = (rho/rho_1) * p + r; in pseudo_inverse() 114 Scalar rho = 1.0, rho_1, lambda, gamma; in constrained_cg() local 157 rho_1 = rho; in constrained_cg() 164 else gamma = (std::max)(0.0, (rho - old_z.dot(z)) / rho_1); in constrained_cg()
|