/dports/science/octopus/octopus-10.5/src/math/ |
H A D | lalg_adv.F90 | 208 SAFE_ALLOCATE(evectors(1:nn, 1:nn)) 214 evectors(1:nn, 1:nn) = aa(1:nn, 1:nn) 231 evectors(1:nn, 1:nn) = aa(1:nn, 1:nn) 237 ex(1:nn, 1:nn) = evectors(1:nn, 1:nn) 239 call lalg_inverter(nn, evectors) 242 evectors(1:nn, ii) = zevalues(1:nn)*evectors(1:nn, ii) 250 SAFE_DEALLOCATE_A(evectors) 286 SAFE_ALLOCATE(evectors(1:nn, 1:nn)) 292 evectors = aa 321 call lalg_inverter(nn, evectors) [all …]
|
/dports/science/lammps/lammps-stable_29Sep2021/src/BODY/ |
H A D | body_nparticle.cpp | 140 double evectors[3][3]; in data_body() local 141 int ierror = MathEigen::jacobi3(tensor,inertia,evectors); in data_body() 159 ex_space[0] = evectors[0][0]; in data_body() 160 ex_space[1] = evectors[1][0]; in data_body() 161 ex_space[2] = evectors[2][0]; in data_body() 162 ey_space[0] = evectors[0][1]; in data_body() 163 ey_space[1] = evectors[1][1]; in data_body() 164 ey_space[2] = evectors[2][1]; in data_body() 165 ez_space[0] = evectors[0][2]; in data_body() 166 ez_space[1] = evectors[1][2]; in data_body() [all …]
|
H A D | body_rounded_polygon.cpp | 202 double evectors[3][3]; in data_body() local 203 int ierror = MathEigen::jacobi3(tensor,inertia,evectors); in data_body() 221 ex_space[0] = evectors[0][0]; in data_body() 222 ex_space[1] = evectors[1][0]; in data_body() 223 ex_space[2] = evectors[2][0]; in data_body() 224 ey_space[0] = evectors[0][1]; in data_body() 225 ey_space[1] = evectors[1][1]; in data_body() 226 ey_space[2] = evectors[2][1]; in data_body() 227 ez_space[0] = evectors[0][2]; in data_body() 228 ez_space[1] = evectors[1][2]; in data_body() [all …]
|
H A D | body_rounded_polyhedron.cpp | 246 double evectors[3][3]; in data_body() local 247 int ierror = MathEigen::jacobi3(tensor,inertia,evectors); in data_body() 265 ex_space[0] = evectors[0][0]; in data_body() 266 ex_space[1] = evectors[1][0]; in data_body() 267 ex_space[2] = evectors[2][0]; in data_body() 268 ey_space[0] = evectors[0][1]; in data_body() 269 ey_space[1] = evectors[1][1]; in data_body() 270 ey_space[2] = evectors[2][1]; in data_body() 271 ez_space[0] = evectors[0][2]; in data_body() 272 ez_space[1] = evectors[1][2]; in data_body() [all …]
|
/dports/math/cgal/CGAL-5.3/include/CGAL/Classification/ |
H A D | Local_eigen_analysis.h | 550 std::array<float, 9> evectors = make_array( 0.f, 0.f, 0.f, in compute() local 555 (covariance, evalues, evectors); in compute() 566 sum = evectors[0] + evectors[1] + evectors[2]; in compute() 569 evectors[i] = evectors[i] / sum; in compute() 570 m_content->smallest_eigenvectors[index] = make_array( float(evectors[0]), float(evectors[1]) ); in compute() 627 std::array<float, 9> evectors = {{ 0.f, 0.f, 0.f, in compute_triangles() local 632 (covariance, evalues, evectors); in compute_triangles() 644 sum = evectors[0] + evectors[1] + evectors[2]; in compute_triangles() 647 evectors[i] = evectors[i] / sum; in compute_triangles() 648 …_eigenvectors[get(get(CGAL::face_index,g), query)] = {{ float(evectors[0]), float(evectors[1]), }}; in compute_triangles()
|
/dports/science/lammps/lammps-stable_29Sep2021/src/ |
H A D | compute_omega_chunk.cpp | 203 double ione[3][3],inverse[3][3],evectors[3][3]; in compute_array() local 256 int ierror = MathEigen::jacobi3(ione,idiag,evectors); in compute_array() 260 ex[0] = evectors[0][0]; in compute_array() 261 ex[1] = evectors[1][0]; in compute_array() 262 ex[2] = evectors[2][0]; in compute_array() 263 ey[0] = evectors[0][1]; in compute_array() 264 ey[1] = evectors[1][1]; in compute_array() 265 ey[2] = evectors[2][1]; in compute_array() 266 ez[0] = evectors[0][2]; in compute_array() 267 ez[1] = evectors[1][2]; in compute_array() [all …]
|
H A D | atom_vec_tri.cpp | 552 double tensor[3][3],evectors[3][3]; in data_atom_bonus() local 560 int ierror = MathEigen::jacobi3(tensor,bonus[nlocal_bonus].inertia,evectors); in data_atom_bonus() 564 ex_space[0] = evectors[0][0]; in data_atom_bonus() 565 ex_space[1] = evectors[1][0]; in data_atom_bonus() 566 ex_space[2] = evectors[2][0]; in data_atom_bonus() 567 ey_space[0] = evectors[0][1]; in data_atom_bonus() 568 ey_space[1] = evectors[1][1]; in data_atom_bonus() 569 ey_space[2] = evectors[2][1]; in data_atom_bonus() 570 ez_space[0] = evectors[0][2]; in data_atom_bonus() 571 ez_space[1] = evectors[1][2]; in data_atom_bonus() [all …]
|
H A D | group.cpp | 1690 double evectors[3][3],inverse[3][3]; in omega() local 1740 int ierror = MathEigen::jacobi3(inertia, idiag, evectors); in omega() 1744 ex[0] = evectors[0][0]; in omega() 1745 ex[1] = evectors[1][0]; in omega() 1746 ex[2] = evectors[2][0]; in omega() 1747 ey[0] = evectors[0][1]; in omega() 1748 ey[1] = evectors[1][1]; in omega() 1749 ey[2] = evectors[2][1]; in omega() 1750 ez[0] = evectors[0][2]; in omega() 1751 ez[1] = evectors[1][2]; in omega() [all …]
|
/dports/math/reduce/Reduce-svn5758-src/psl/dist/util/ |
H A D | f-evectors.sl | 39 %%% evectors that runs fast. To use this facility, LOAD (don't 40 %%% IMPORT) it at compiletime. It does an (on fast-evectors), 59 (put 'egetv 'fast-flag '*fast-evectors) 60 (put 'eputv 'fast-flag '*fast-evectors) 61 (put 'esizev 'fast-flag '*fast-evectors) 62 (put 'eupbv 'fast-flag '*fast-evectors) 79 (loadtime (on fast-evectors))
|
H A D | Makefile | 67 $(PL)/f-evectors.b $(PL)/fast-int.b $(PL)/pipes.b \ 159 $(PL)/f-evectors.b: $(PU)/f-evectors.sl $(PL)/if-system.b 160 $(PDIST)/make-lap $(PU) f-evectors 259 $(PL)/f-vectors.b $(PL)/f-evectors.b $(PL)/strings.b \ 264 $(PL)/f-vectors.b $(PL)/f-evectors.b \ 269 $(PL)/f-vectors.b $(PL)/f-evectors.b \ 274 $(PL)/f-vectors.b $(PL)/f-evectors.b \ 279 $(PL)/f-vectors.b $(PL)/f-evectors.b \ 284 $(PL)/f-vectors.b $(PL)/f-evectors.b \ 298 $(PL)/numeric-ops.b $(PL)/f-vectors.b $(PL)/f-evectors.b \
|
/dports/math/lis/lis-2.0.30/graphics/ |
H A D | gnuplot5.sh | 19 $srcdir/etest5 $srcdir/testmat.mtx /dev/null evectors.mtx /dev/null /dev/null -e si -ie ii -ss $sub… 24 filename=evectors.mtx 35 rm evectors.mtx
|
H A D | gnuplot5.bat | 20 @%bindir%\etest5.exe %srcdir%\testmat.mtx nul evectors.mtx nul nul -e si -ie ii -ss %subspace% 25 @set filename=evectors.mtx 34 @del evectors.mtx
|
/dports/science/liggghts/LIGGGHTS-PUBLIC-3.8.0-26-g6e873439/src/ |
H A D | fix_template_multisphere.cpp | 444 double evectors[3][3]; in calc_eigensystem() local 446 int ierror = MathExtra::jacobi(moi_,static_cast<double*>(inertia_),evectors); in calc_eigensystem() 449 ex_space_[0] = evectors[0][0]; in calc_eigensystem() 450 ex_space_[1] = evectors[1][0]; in calc_eigensystem() 451 ex_space_[2] = evectors[2][0]; in calc_eigensystem() 452 ey_space_[0] = evectors[0][1]; in calc_eigensystem() 453 ey_space_[1] = evectors[1][1]; in calc_eigensystem() 454 ey_space_[2] = evectors[2][1]; in calc_eigensystem() 455 ez_space_[0] = evectors[0][2]; in calc_eigensystem() 456 ez_space_[1] = evectors[1][2]; in calc_eigensystem() [all …]
|
H A D | math_extra.cpp | 133 int jacobi(double matrix[3][3], double *evalues, double evectors[3][3]) in jacobi() 139 for (j = 0; j < 3; j++) evectors[i][j] = 0.0; in jacobi() 140 evectors[i][i] = 1.0; in jacobi() 183 for (k = 0; k < 3; k++) rotate(evectors,k,i,k,j,s,tau); in jacobi() 196 int jacobi(double **matrix, double *evalues, double **evectors) in jacobi() argument 202 for (j = 0; j < 3; j++) evectors[i][j] = 0.0; in jacobi() 203 evectors[i][i] = 1.0; in jacobi() 246 for (k = 0; k < 3; k++) rotate(evectors,k,i,k,j,s,tau); in jacobi()
|
/dports/science/octopus/octopus-10.5/src/states/ |
H A D | modelmb_density_matrix_inc.F90 | 164 R_TYPE, allocatable :: densmatr(:, :), evectors(:, :) 200 SAFE_ALLOCATE(evectors(1:mb_1part%npt, 1:mb_1part%npt)) 218 evectors = densmatr_tmp 219 call lalg_eigensolve(mb_1part%npt, evectors, evalues, bof, err_code) 224 evectors = evectors/sqrt(mb_1part%vol_elem_1part) 250 write(iunit,'(es11.3,es11.3)') evectors(ll,jj) 307 SAFE_DEALLOCATE_A(evectors)
|
/dports/science/hypre/hypre-2.23.0/src/FEI_mv/femli/ |
H A D | driver_util.c | 45 double *matrix, *evalues, *evectors, *daux1, *daux2; in testEigen() local 56 evectors = hypre_TAlloc(double, mDim * mDim , HYPRE_MEMORY_HOST); in testEigen() 60 mli_computespectrum_(&mDim, &mDim, matrix, evalues, &matz, evectors, in testEigen() 64 hypre_TFree(evectors, HYPRE_MEMORY_HOST); in testEigen()
|
H A D | mli_amgsa_elem.cxx | 53 double *evalues, *evectors, *dAux1, *dAux2; in setupUsingFEData() local 211 evectors = NULL; in setupUsingFEData() 254 evectors = new double[macroMatDim*macroMatDim]; in setupUsingFEData() 290 &matz, evectors, dAux1, dAux2, &ierr); in setupUsingFEData() 302 evectors[j*blockSize+k+i*macroMatDim]; in setupUsingFEData() 309 delete [] evectors; in setupUsingFEData()
|
/dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/nwpw/ |
H A D | nwpw_vectors_input.F | 207 $ ('nwpw_vectors_output:evectors input filename',0,INPUT_ERR) 215 $ ('nwpw_vectors_output:evectors output filename',0,INPUT_ERR) 289 > 'input evectors' 298 > 'output evectors'
|
/dports/graphics/opencv/opencv-4.5.3/modules/core/test/ |
H A D | test_eigen.cpp | 86 …bool check_pair_count(const cv::Mat& src, const cv::Mat& evalues, const cv::Mat& evectors, int low… 181 …k_pair_count(const cv::Mat& src, const cv::Mat& evalues, const cv::Mat& evectors, int low_index, i… in check_pair_count() argument 186 if (!(evectors.rows == right_eigen_pair_count && evectors.cols == right_eigen_pair_count)) in check_pair_count() 188 …std::cout << endl; std::cout << "Checking sizes of eigen vectors matrix " << evectors << "..." << … in check_pair_count() 189 …std::cout << "Number of rows: " << evectors.rows << " Number of cols: " << evectors.cols << endl; in check_pair_count()
|
/dports/science/octopus/octopus-10.5/src/system/ |
H A D | subspace_inc.F90 | 168 R_TYPE, allocatable :: hs(:, :), hpsi(:, :, :), evectors(:, :) 248 SAFE_ALLOCATE(evectors(1:nrow, 1:ncol)) 277 call elpa%eigenvectors(hs, eigenval, evectors, info) 295 w = eigenval(1), z = evectors(1, 1), iz = 1, jz = 1, descz = hs_desc(1), & 310 w = eigenval(1), z = evectors(1, 1), iz = 1, jz = 1, descz = hs_desc(1), & 324 …w = eigenval(1), z = evectors(1, 1), iz = 1, jz = 1, descz = hs_desc(1), work = rttmp, lwork = -1,… 335 …w = eigenval(1), z = evectors(1, 1), iz = 1, jz = 1, descz = hs_desc(1), work = work(1), lwork = l… 367 evectors(1, 1), 1, 1, hs_desc(1), &
|
/dports/science/lammps/lammps-stable_29Sep2021/src/EXTRA-COMPUTE/ |
H A D | compute_gyration_shape.cpp | 89 double ione[3][3], evalues[3], evectors[3][3]; in compute_vector() local 98 int ierror = MathEigen::jacobi3(ione,evalues,evectors); in compute_vector()
|
H A D | compute_gyration_shape_chunk.cpp | 117 double ione[3][3], evalues[3], evectors[3][3]; in compute_array() local 128 int ierror = MathEigen::jacobi3(ione,evalues,evectors); in compute_array()
|
/dports/science/liggghts/LIGGGHTS-PUBLIC-3.8.0-26-g6e873439/src/POEMS/ |
H A D | fix_poems.cpp | 628 double **tensor,**evectors; in init() local 646 ex_space[isegment][0] = evectors[0][0]; in init() 647 ex_space[isegment][1] = evectors[1][0]; in init() 648 ex_space[isegment][2] = evectors[2][0]; in init() 650 ey_space[isegment][0] = evectors[0][1]; in init() 651 ey_space[isegment][1] = evectors[1][1]; in init() 652 ey_space[isegment][2] = evectors[2][1]; in init() 654 ez_space[isegment][0] = evectors[0][2]; in init() 655 ez_space[isegment][1] = evectors[1][2]; in init() 696 memory->destroy(evectors); in init() [all …]
|
/dports/science/lammps/lammps-stable_29Sep2021/src/POEMS/ |
H A D | fix_poems.cpp | 480 double **tensor,**evectors; in init() local 498 ex_space[ibody][0] = evectors[0][0]; in init() 499 ex_space[ibody][1] = evectors[1][0]; in init() 500 ex_space[ibody][2] = evectors[2][0]; in init() 502 ey_space[ibody][0] = evectors[0][1]; in init() 503 ey_space[ibody][1] = evectors[1][1]; in init() 504 ey_space[ibody][2] = evectors[2][1]; in init() 506 ez_space[ibody][0] = evectors[0][2]; in init() 507 ez_space[ibody][1] = evectors[1][2]; in init() 508 ez_space[ibody][2] = evectors[2][2]; in init() [all …]
|
/dports/science/jmol/jmol-14.32.7/src/org/jmol/adapter/readers/xtal/ |
H A D | CrystalReader.java | 714 P3[] evectors = new P3[3]; in processNextCriticalPoint() local 715 evectors[0] = P3.new3((float) ev[0][0], (float) ev[0][1], in processNextCriticalPoint() 717 evectors[1] = P3.new3((float) ev[1][0], (float) ev[1][1], in processNextCriticalPoint() 719 evectors[2] = P3.new3((float) ev[2][0], (float) ev[2][1], in processNextCriticalPoint() 721 System.out.println("evpts " + evectors[0] + " " + evectors[1] + " " in processNextCriticalPoint() 722 + evectors[2]); in processNextCriticalPoint() 723 m.put("eigenvectors", evectors); in processNextCriticalPoint() 724 Tensor t = new Tensor().setFromEigenVectors(evectors, evalues, "cp", in processNextCriticalPoint()
|