/dports/math/R-cran-expm/expm/src/ |
H A D | logm-eigen.c | 46 …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 D | expm-eigen.c | 42 …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 D | benchEigenSolver.cpp | 96 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 D | benchEigenSolver.cpp | 96 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 D | benchEigenSolver.cpp | 96 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 D | benchEigenSolver.cpp | 96 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 D | benchEigenSolver.cpp | 96 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 D | benchEigenSolver.cpp | 96 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 D | NKM_SurfMat.hpp | 94 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 D | NKM_SurfMat.cpp | 765 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 D | tolEllipse.R | 61 eigvect <- eA$vectors[, order(ev)[2]] functionVar 65 alfa <- atan(eigvect[2]/eigvect[1])
|
/dports/math/gmm++/gmm-5.4/include/gmm/ |
H A D | gmm_range_basis.h | 58 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 D | gmm_dense_qr.h | 393 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 D | gmm_dense_qr.h | 392 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 D | cs_turbulence_model.c | 2078 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 D | drawdf.mac | 245 for eigvect in eigvects do ( 247 solution_plot(pt + 0.001 * normalize_vec(eigvect)), 248 solution_plot(pt - 0.001 * normalize_vec(eigvect)),
|