Home
last modified time | relevance | path

Searched refs:eigvect (Results 1 – 16 of 16) sorted by relevance

/dports/math/R-cran-expm/expm/src/
H A Dlogm-eigen.c46 …Rcomplex *eigvect = (Rcomplex *) R_alloc(nsqr, sizeof(Rcomplex)); /* (right) eigenvectors matrix */ in logm_eigen() local
86 eigvect[i * n + j].r = right[i * n + j]; in logm_eigen()
87 eigvect[i * n + j].i = right[(i+1) * n + j]; in logm_eigen()
95 eigvect[i * n + j].r = right[(i-1) * n + j]; in logm_eigen()
96 eigvect[i * n + j].i = -right[i * n + j]; in logm_eigen()
102 eigvect[i * n + j].r = right[i * n + j]; in logm_eigen()
103 eigvect[i * n + j].i = 0.0; in logm_eigen()
112 Memcpy(ctmp, eigvect, nsqr); in logm_eigen()
117 F77_CALL(zgesv)(&n, &n, eigvect, &n, ipiv, eigvectinv, &n, &info); in logm_eigen()
171 Memcpy(eigvect, ctmp, nsqr); in logm_eigen()
[all …]
H A Dexpm-eigen.c42 …Rcomplex *eigvect = (Rcomplex *) R_alloc(nsqr, sizeof(Rcomplex)); /* (right) eigenvectors matrix */ in expm_eigen() local
82 eigvect[i * n + j].r = right[i * n + j]; in expm_eigen()
83 eigvect[i * n + j].i = right[(i+1) * n + j]; in expm_eigen()
91 eigvect[i * n + j].r = right[(i-1) * n + j]; in expm_eigen()
92 eigvect[i * n + j].i = -right[i * n + j]; in expm_eigen()
98 eigvect[i * n + j].r = right[i * n + j]; in expm_eigen()
99 eigvect[i * n + j].i = 0.0; in expm_eigen()
108 Memcpy(ctmp, eigvect, nsqr); in expm_eigen()
113 F77_CALL(zgesv)(&n, &n, eigvect, &n, ipiv, eigvectinv, &n, &info); in expm_eigen()
168 Memcpy(eigvect, ctmp, nsqr); in expm_eigen()
[all …]
/dports/math/eigen3/eigen-3.3.9/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/math/stanmath/math-4.2.0/lib/eigen_3.3.9/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/math/libmesh/libmesh-1.6.2/contrib/eigen/3.2.9/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/math/libsemigroups/libsemigroups-1.3.7/extern/eigen-3.3.7/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/misc/opennn/opennn-5.0.5/eigen/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/devel/taskflow/taskflow-3.2.0/3rd-party/eigen-3.3.7/bench/
H A DbenchEigenSolver.cpp96 gmm::dense_matrix<Scalar> eigvect(covMat.rows(),covMat.cols()); in benchEigenSolver() local
104 gmm::symmetric_qr_algorithm(gmmCovMat, eigval, eigvect); in benchEigenSolver()
105 acc += eigvect(r,c); in benchEigenSolver()
135 gsl_matrix* eigvect = gsl_matrix_alloc(covMat.rows(),covMat.cols()); in benchEigenSolver() local
150 gsl_eigen_symmv(gslCopy, eigval, eigvect, eisymm); in benchEigenSolver()
151 acc += gsl_matrix_get(eigvect,r,c); in benchEigenSolver()
173 gsl_matrix_free(eigvect); in benchEigenSolver()
/dports/science/dakota/dakota-6.13.0-release-public.src-UI/packages/surfpack/src/surfaces/nkm/
H A DNKM_SurfMat.hpp94 MtxDbl& eig_sym(MtxDbl& eigvect, MtxDbl& eigval, const MtxDbl& A, char jobz='V');
98 MtxDbl eigvect; in eig_sym() local
99 return (eig_sym(eigvect,eigval,A,'N')); in eig_sym()
H A DNKM_SurfMat.cpp765 MtxDbl& eig_sym(MtxDbl& eigvect, MtxDbl& eigval, const MtxDbl& A, char jobz) { in eig_sym() argument
767 eigvect.copy(A); in eig_sym()
768 int nrowsA=static_cast<int>(eigvect.getNRows()); in eig_sym()
771 assert((0<nrowsA)&&(nrowsA==eigvect.getNCols())); in eig_sym()
777 DSYEV_F77(&jobz, &uplo, &nrowsA, eigvect.ptr(0,0), &lda, eigval.ptr(0,0), in eig_sym()
784 DSYEV_F77(&jobz, &uplo, &nrowsA, eigvect.ptr(0,0), &lda, eigval.ptr(0,0), in eig_sym()
813 MtxDbl eigvect(nrowsA,nrowsA), eigval(nrowsA,1); in pseudo_inverse_sym() local
814 eig_sym(eigvect, eigval, A); in pseudo_inverse_sym()
833 eigvect_eigval_inv(i,j)=eigvect(i,j)*eigval_inv; in pseudo_inverse_sym()
837 eigvect_eigval_inv(i,j)=eigvect(i,j)=0.0; in pseudo_inverse_sym()
[all …]
/dports/math/R-cran-robustbase/robustbase/R/
H A DtolEllipse.R61 eigvect <- eA$vectors[, order(ev)[2]] functionVar
65 alfa <- atan(eigvect[2]/eigvect[1])
/dports/math/gmm++/gmm-5.4/include/gmm/
H A Dgmm_range_basis.h58 MAT1 &eigvect = const_cast<MAT1 &>(eigvect_); variable
61 if (compvect) gmm::copy(identity_matrix(), eigvect);
70 sub_interval SUBI(p, n-p-q), SUBJ(0, mat_ncols(eigvect)), SUBK(p, n-p-q);
75 sub_matrix(eigvect, SUBJ, SUBK), compvect);
105 dense_matrix<T> eigvect(restart, restart);
139 …sdiag.resize(restart); eigval.resize(restart); diag.resize(restart); gmm::resize(eigvect, restart,…
163 gmm::resize(eigvect, eff_restart, eff_restart); gmm::resize(lv, nc_r, eff_restart);
165 tridiag_qr_algorithm(diag, sdiag, eigval, eigvect, true);
174 gmm::mult(lv, mat_col(eigvect, num), v);
H A Dgmm_dense_qr.h393 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
401 Hessenberg_reduction(A1, eigvect, compvect);
407 if (compvect) { gmm::mult(eigvect, Q, R); gmm::copy(R, eigvect); }
706 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
710 if (compvect) gmm::copy(identity_matrix(), eigvect);
715 Householder_tridiagonalization(Tri, eigvect, compvect);
724 sub_matrix(eigvect, SUBJ, SUBK), compvect);
741 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
746 if (compvect) gmm::copy(identity_matrix(), eigvect);
752 Householder_tridiagonalization(Tri, eigvect, compvect);
[all …]
/dports/cad/gmsh/gmsh-4.9.2-source/contrib/gmm/
H A Dgmm_dense_qr.h392 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
400 Hessenberg_reduction(A1, eigvect, compvect);
406 if (compvect) { gmm::mult(eigvect, Q, R); gmm::copy(R, eigvect); }
705 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
709 if (compvect) gmm::copy(identity_matrix(), eigvect);
714 Householder_tridiagonalization(Tri, eigvect, compvect);
723 sub_matrix(eigvect, SUBJ, SUBK), compvect);
740 MAT2 &eigvect = const_cast<MAT2 &>(eigvect_); variable
745 if (compvect) gmm::copy(identity_matrix(), eigvect);
751 Householder_tridiagonalization(Tri, eigvect, compvect);
[all …]
/dports/science/code_saturne/code_saturne-7.1.0/src/turb/
H A Dcs_turbulence_model.c2078 cs_real_33_t eigvect = {{1.0,0.0,0.0},{0.0,1.0,0.0},{0.0,0.0,1.0}}; in cs_clip_turbulent_fluxes() local
2117 cs_math_33_eig_val_vec(rij,tol_jacobi,eigval,eigvect); in cs_clip_turbulent_fluxes()
2118 cs_math_33_3_product(eigvect,cvar_rit[cell_id],rot_rit); in cs_clip_turbulent_fluxes()
2139 cs_math_33t_3_product(eigvect,rot_rit,cvar_rit[cell_id]); in cs_clip_turbulent_fluxes()
/dports/math/maxima/maxima-5.43.2/share/diffequations/
H A Ddrawdf.mac245 for eigvect in eigvects do (
247 solution_plot(pt + 0.001 * normalize_vec(eigvect)),
248 solution_plot(pt - 0.001 * normalize_vec(eigvect)),