/dports/science/jmol/jmol-14.32.7/src/org/jmol/util/ |
H A D | Tensor.java | 198 return eigenValues; in getInfo() 277 return (eigenValues[0] + eigenValues[1] + eigenValues[2]) / 3; in isotropy() 295 return Math.abs(eigenValues[2] - eigenValues[0]); in span() 333 return eigenValues[2] - (eigenValues[0] + eigenValues[1]) / 2; in anisotropy() 378 t.eigenValues = eigenValues; in copyTensor() 592 eigenValues = values; in newTensorType() 600 + (eigenValues[1] >= 0 ? 2 : 0) + (eigenValues[2] >= 0 ? 4 : 0); in newTensorType() 681 && eigenValues[2] - eigenValues[1] < eigenValues[1] - eigenValues[0]) { in sortAndNormalize() 686 eigenValues[0] = eigenValues[2]; in sortAndNormalize() 696 float f = Math.abs(eigenValues[0] + eigenValues[1] + eigenValues[2]); in isEquiv() [all …]
|
/dports/math/dune-common/dune-common-df65b1282ea89ad40d2cb6565983f7e633ccce31/dune/common/ |
H A D | fmatrixev.hh | 296 eigenValues[0] = matrix[0][0]; in eigenValuesVectorsImpl() 322 temp[0][0] -= eigenValues[0]; in eigenValuesVectorsImpl() 371 if (eigenValues[0] > eigenValues[1]) in eigenValuesVectorsImpl() 373 std::swap(eigenValues[0], eigenValues[1]); in eigenValuesVectorsImpl() 376 if (eigenValues[1] > eigenValues[2]) in eigenValuesVectorsImpl() 378 std::swap(eigenValues[1], eigenValues[2]); in eigenValuesVectorsImpl() 381 if (eigenValues[0] > eigenValues[1]) in eigenValuesVectorsImpl() 383 std::swap(eigenValues[0], eigenValues[1]); in eigenValuesVectorsImpl() 393 Vector eval(eigenValues); in eigenValuesVectorsImpl() 418 eigenValues *= maxAbsElement; in eigenValuesVectorsImpl() [all …]
|
/dports/science/berkeleygw/BGW-2.0.0/MeanField/EPM/ |
H A D | epm.f90 | 213 eigenStates%eigenValues(n,k) = eigenStates%eigenValues(n,k) * hartree2eV 232 homo=eigenStates%eigenValues(n,k) 349 …dv1=-(eigenStates%eigenValues(n, 2)+eigenStates%eigenValues(n, 3)-2.0d0*eigenStates%eigenValues(n,… 351 …dv2=-(eigenStates%eigenValues(n, 4)+eigenStates%eigenValues(n, 5)-2.0d0*eigenStates%eigenValues(n,… 353 …dv3=-(eigenStates%eigenValues(n, 6)+eigenStates%eigenValues(n, 7)-2.0d0*eigenStates%eigenValues(n,… 363 …dv1=-(eigenStates%eigenValues(n, 2)+eigenStates%eigenValues(n, 3)-2.0d0*eigenStates%eigenValues(n,… 365 …dv2=-(eigenStates%eigenValues(n, 4)+eigenStates%eigenValues(n, 5)-2.0d0*eigenStates%eigenValues(n,… 367 …dv3=-(eigenStates%eigenValues(n, 6)+eigenStates%eigenValues(n, 7)-2.0d0*eigenStates%eigenValues(n,… 377 …dc1=(eigenStates%eigenValues(n, 9)+eigenStates%eigenValues(n,10)-2.0d0*eigenStates%eigenValues(n, … 379 …dc2=(eigenStates%eigenValues(n,11)+eigenStates%eigenValues(n,12)-2.0d0*eigenStates%eigenValues(n, … [all …]
|
H A D | eigenStatesMod.f90 | 25 real(double), pointer :: eigenValues(:,:) 54 nullify( eigenStates%eigenValues ) 64 allocate( eigenStates%eigenValues( eigenStates%numStates, & 101 deallocate( eigenStates%eigenVectors, eigenStates%eigenValues, & 130 write(*,'(i5,f23.16)') kpt, eigenStates%eigenValues(j, kpt)
|
/dports/misc/elki/elki-release0.7.1-1166-gfb1fffdf3/elki-core-math/src/main/java/de/lmu/ifi/dbs/elki/math/linearalgebra/pca/filter/ |
H A D | DropEigenPairFilter.java | 68 public int filter(double[] eigenValues) { in filter() argument 69 int contrastMaximum = eigenValues.length; in filter() 74 for(int i = 0; i < eigenValues.length; i++) { in filter() 75 eigenValueSum += eigenValues[i]; in filter() 78 final double weakEigenvalue = walpha * eigenValueSum / eigenValues.length; in filter() 80 double prev_sum = eigenValues[eigenValues.length - 1]; in filter() 82 for(int i = eigenValues.length - 2; i >= 0; i++) { in filter() 83 double curr_sum = prev_sum + eigenValues[i]; in filter() 84 double curr_rel = eigenValues[i] / curr_sum * i; in filter() 86 if(eigenValues[i] >= weakEigenvalue) { in filter()
|
H A D | ProgressiveEigenPairFilter.java | 115 public int filter(double[] eigenValues) { in filter() argument 118 for(int i = 0; i < eigenValues.length; i++) { in filter() 119 totalSum += eigenValues[i]; in filter() 121 double expectedVariance = totalSum / eigenValues.length * walpha; in filter() 125 for(int i = 0; i < eigenValues.length - 1; i++) { in filter() 127 if(eigenValues[i] < expectedVariance) { in filter() 130 currSum += eigenValues[i]; in filter() 132 double alpha = 1.0 - (1.0 - palpha) * (1.0 - (i + 1) / (double) eigenValues.length); in filter() 140 return eigenValues.length; in filter()
|
H A D | RelativeEigenPairFilter.java | 74 public int filter(double[] eigenValues) { in filter() argument 77 double eigenValueSum = eigenValues[eigenValues.length - 1]; in filter() 78 for(int i = eigenValues.length - 2; i >= 0; i--) { in filter() 79 eigenValueSum += eigenValues[i]; in filter() 80 double needEigenvalue = eigenValueSum / (eigenValues.length - i) * ralpha; in filter() 81 if(eigenValues[i] >= needEigenvalue) { in filter() 85 return eigenValues.length; in filter()
|
H A D | WeakEigenPairFilter.java | 63 public int filter(double[] eigenValues) { in filter() argument 66 for(int i = 0; i < eigenValues.length; i++) { in filter() 67 totalSum += eigenValues[i]; in filter() 69 double expectEigenvalue = totalSum / eigenValues.length * walpha; in filter() 72 for(int i = 0; i < eigenValues.length; i++) { in filter() 73 if(eigenValues[i] <= expectEigenvalue) { in filter() 80 return eigenValues.length; in filter()
|
H A D | PercentageEigenPairFilter.java | 65 public int filter(double[] eigenValues) { in filter() argument 68 for(int i = 0; i < eigenValues.length; i++) { in filter() 69 totalSum += eigenValues[i]; in filter() 74 for(int i = 0; i < eigenValues.length; i++) { in filter() 75 currSum += eigenValues[i]; in filter() 80 return eigenValues.length; in filter()
|
H A D | LimitEigenPairFilter.java | 71 public int filter(double[] eigenValues) { in filter() argument 73 double limit = absolute ? delta : eigenValues[0] * delta; in filter() 76 for(int i = 0; i < eigenValues.length; i++) { in filter() 77 double eigenValue = Math.abs(eigenValues[i]); in filter() 83 return eigenValues.length; in filter()
|
/dports/math/mlpack/mlpack-3.4.2/src/mlpack/core/data/scaler_methods/ |
H A D | pca_whitening.hpp | 75 eig_sym(eigenValues, eigenVectors, mlpack::math::ColumnCovariance( in Fit() 77 eigenValues += epsilon; in Fit() 89 if (eigenValues.is_empty() || eigenVectors.is_empty()) in Transform() 96 output = arma::diagmat(1.0 / (arma::sqrt(eigenValues))) * eigenVectors.t() in Transform() 109 output = arma::diagmat(arma::sqrt(eigenValues)) * inv(eigenVectors.t()) in InverseTransform() 117 const arma::vec& EigenValues() const { return eigenValues; } in EigenValues() 126 ar & BOOST_SERIALIZATION_NVP(eigenValues); in serialize() 140 arma::vec eigenValues; member in mlpack::data::PCAWhitening
|
/dports/science/dftbplus/dftbplus-19.1/prog/modes/ |
H A D | modes.F90 | 28 real(dp), allocatable :: eigenValues(:) variable 37 allocate(eigenValues(3 * nMovedAtom)) 60 call heev(dynMatrix,eigenValues,'U','V') 62 call heev(dynMatrix,eigenValues,'U','N') 66 eigenValues = sign(sqrt(abs(eigenValues)),eigenValues) 69 write(stdout, '(i5,f8.2)')ii,eigenValues(ii)*Hartree__cm 75 call taggedWriter%write(12, "frequencies", eigenValues) 125 write(123,*)'Eigenmode',iMode,eigenValues(iMode)*Hartree__cm,'cm-1' 142 write(123,*)'Eigenmode',iMode,eigenValues(iMode)*Hartree__cm,'cm-1'
|
/dports/science/InsightToolkit/ITK-5.0.1/Modules/Core/Common/include/ |
H A D | itkSymmetricEigenAnalysis.h | 78 { return std::abs(eigenValues[a]) < std::abs(eigenValues[b]); } in sortEigenValuesByMagnitude() 80 auto tmpCopy = eigenValues; in sortEigenValuesByMagnitude() 488 auto copyEigenValues = eigenValues; 506 EigenValues[row] = eigenValues[row]; 552 auto copyEigenValues = eigenValues; 622 EigenValues[i] = eigenValues[i]; 662 EigenValues[i] = eigenValues[i]; 824 auto copyEigenValues = eigenValues; 886 auto copyEigenValues = eigenValues; 947 EigenValues[i] = eigenValues[i]; [all …]
|
/dports/emulators/mess/mame-mame0226/3rdparty/bimg/3rdparty/nvtt/nvmath/ |
H A D | fitting.cpp | 458 if (eigenValues[2] > eigenValues[0] && eigenValues[2] > eigenValues[1]) in eigenSolveSymmetric3() 460 swap(eigenValues[0], eigenValues[2]); in eigenSolveSymmetric3() 463 if (eigenValues[1] > eigenValues[0]) in eigenSolveSymmetric3() 465 swap(eigenValues[0], eigenValues[1]); in eigenSolveSymmetric3() 468 if (eigenValues[2] > eigenValues[1]) in eigenSolveSymmetric3() 470 swap(eigenValues[1], eigenValues[2]); in eigenSolveSymmetric3() 474 nvDebugCheck(eigenValues[0] >= eigenValues[1] && eigenValues[0] >= eigenValues[2]); in eigenSolveSymmetric3() 653 if (eigenValues[j] > eigenValues[i]) in eigenSolveSymmetric4() 655 swap(eigenValues[i], eigenValues[j]); in eigenSolveSymmetric4() 661 …nvDebugCheck(eigenValues[0] >= eigenValues[1] && eigenValues[0] >= eigenValues[2] && eigenValues[0… in eigenSolveSymmetric4() [all …]
|
/dports/emulators/mame/mame-mame0226/3rdparty/bimg/3rdparty/nvtt/nvmath/ |
H A D | fitting.cpp | 458 if (eigenValues[2] > eigenValues[0] && eigenValues[2] > eigenValues[1]) in eigenSolveSymmetric3() 460 swap(eigenValues[0], eigenValues[2]); in eigenSolveSymmetric3() 463 if (eigenValues[1] > eigenValues[0]) in eigenSolveSymmetric3() 465 swap(eigenValues[0], eigenValues[1]); in eigenSolveSymmetric3() 468 if (eigenValues[2] > eigenValues[1]) in eigenSolveSymmetric3() 470 swap(eigenValues[1], eigenValues[2]); in eigenSolveSymmetric3() 474 nvDebugCheck(eigenValues[0] >= eigenValues[1] && eigenValues[0] >= eigenValues[2]); in eigenSolveSymmetric3() 653 if (eigenValues[j] > eigenValues[i]) in eigenSolveSymmetric4() 655 swap(eigenValues[i], eigenValues[j]); in eigenSolveSymmetric4() 661 …nvDebugCheck(eigenValues[0] >= eigenValues[1] && eigenValues[0] >= eigenValues[2] && eigenValues[0… in eigenSolveSymmetric4() [all …]
|
/dports/audio/faust/faust-2.37.3/tools/physicalModeling/mesh2faust/vega/libraries/minivector/ |
H A D | mat3d.cpp | 98 Vec3d eigenValues; in SVD() local 101 eigen_sym(normalEq, eigenValues, eigenVectors); in SVD() 124 Sigma[0] = (eigenValues[0] > 0.0) ? sqrt(eigenValues[0]) : 0.0; in SVD() 125 Sigma[1] = (eigenValues[1] > 0.0) ? sqrt(eigenValues[1]) : 0.0; in SVD() 126 Sigma[2] = (eigenValues[2] > 0.0) ? sqrt(eigenValues[2]) : 0.0; in SVD()
|
/dports/math/openturns/openturns-1.18/lib/src/Base/Algo/ |
H A D | KarhunenLoeveP1Algorithm.cxx | 337 Point eigenValues(nev); in run() local 377 eigenValues, in run() 380 LOGDEBUG(OSS(false) << "eigenValues=" << eigenValues); in run() 385 computedVariance += eigenValues[k]; in run() 389 …Variance = computedVariance + (augmentedDimension - eigenValues.getSize()) * eigenValues[eigenValu… in run() 390 LOGDEBUG(OSS(false) << "eigenValues=" << eigenValues); in run() 429 eigenValues[i] = -eigenPairs[i].first; in run() 430 cumulatedVariance += eigenValues[i]; in run() 446 selectedEV.add(eigenValues[K]); in run() 447 selectedVariance += eigenValues[K]; in run() [all …]
|
/dports/science/jstrack/jstrack/extlib/tcllib1.19/math/ |
H A D | pca.tcl | 32 variable eigenValues {} 83 set eigenValues {} 85 lappend eigenValues [expr {$c**2 * $factor}] 92 set retainedValues $eigenValues 113 variable eigenValues 117 return $eigenValues 271 set retainedValues [lrange $eigenValues 0 [expr {$numberUsed-1}]] 277 set retainedValues $eigenValues
|
/dports/devel/tcllib/tcllib-1.20/modules/math/ |
H A D | pca.tcl | 32 variable eigenValues {} 83 set eigenValues {} 85 lappend eigenValues [expr {$c**2 * $factor}] 92 set retainedValues $eigenValues 113 variable eigenValues 117 return $eigenValues 271 set retainedValues [lrange $eigenValues 0 [expr {$numberUsed-1}]] 277 set retainedValues $eigenValues
|
/dports/devel/tcllibc/tcllib-1.20/modules/math/ |
H A D | pca.tcl | 32 variable eigenValues {} 83 set eigenValues {} 85 lappend eigenValues [expr {$c**2 * $factor}] 92 set retainedValues $eigenValues 113 variable eigenValues 117 return $eigenValues 271 set retainedValues [lrange $eigenValues 0 [expr {$numberUsed-1}]] 277 set retainedValues $eigenValues
|
/dports/math/apache-commons-math/commons-math3-3.6.1-src/src/test/java/org/apache/commons/math3/linear/ |
H A D | EigenDecompositionTest.java | 153 double[] eigenValues = decomposition.getRealEigenvalues(); in testMath308() local 196 double[] eigenValues = decomposition.getRealEigenvalues(); in testMathpbx02() local 241 double[] eigenValues = decomposition.getRealEigenvalues(); in testMathpbx03() local 270 double[] eigenValues = ed.getRealEigenvalues(); in testTridiagonal() local 271 Assert.assertEquals(ref.length, eigenValues.length); in testTridiagonal() 297 double[] eigenValues = ed.getRealEigenvalues(); in testEigenvalues() local 298 Assert.assertEquals(refValues.length, eigenValues.length); in testEigenvalues() 300 Assert.assertEquals(refValues[i], eigenValues[i], 3.0e-15); in testEigenvalues() 315 double[] eigenValues = ed.getRealEigenvalues(); in testBigMatrix() local 316 Assert.assertEquals(bigValues.length, eigenValues.length); in testBigMatrix() [all …]
|
/dports/science/InsightToolkit/ITK-5.0.1/Modules/Core/Common/test/ |
H A D | itkSymmetricSecondRankTensorTest.cxx | 214 Double3DTensorType::EigenValuesArrayType eigenValues; in itkSymmetricSecondRankTensorTest() local 218 tensor3D.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkSymmetricSecondRankTensorTest() 222 std::cout << eigenValues << std::endl; in itkSymmetricSecondRankTensorTest() 237 if( std::fabs( expectedValues[i] - eigenValues[i] ) > tolerance ) in itkSymmetricSecondRankTensorTest() 241 std::cerr << "eigenValues = " << eigenValues << std::endl; in itkSymmetricSecondRankTensorTest() 273 tensor3D.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkSymmetricSecondRankTensorTest() 277 std::cout << eigenValues << std::endl; in itkSymmetricSecondRankTensorTest() 294 std::cerr << "eigenValues = " << eigenValues << std::endl; in itkSymmetricSecondRankTensorTest() 327 tensor3D.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkSymmetricSecondRankTensorTest() 331 std::cout << eigenValues << std::endl; in itkSymmetricSecondRankTensorTest() [all …]
|
/dports/science/InsightToolkit/ITK-5.0.1/Modules/Filtering/DiffusionTensorImage/test/ |
H A D | itkDiffusionTensor3DTest.cxx | 231 Double3DTensorType::EigenValuesArrayType eigenValues; in itkDiffusionTensor3DTest() local 234 tensor2.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkDiffusionTensor3DTest() 237 std::cout << eigenValues << std::endl; in itkDiffusionTensor3DTest() 252 if( std::fabs( expectedValues[i] - eigenValues[i] ) > tolerance ) in itkDiffusionTensor3DTest() 256 std::cerr << "eigenValues = " << eigenValues << std::endl; in itkDiffusionTensor3DTest() 276 tensor2.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkDiffusionTensor3DTest() 279 std::cout << eigenValues << std::endl; in itkDiffusionTensor3DTest() 296 std::cerr << "eigenValues = " << eigenValues << std::endl; in itkDiffusionTensor3DTest() 316 tensor2.ComputeEigenAnalysis( eigenValues, eigenVectors ); in itkDiffusionTensor3DTest() 319 std::cout << eigenValues << std::endl; in itkDiffusionTensor3DTest() [all …]
|
/dports/science/jmol/jmol-14.32.7/src/org/jmol/quantum/ |
H A D | NMRNoeMatrix.java | 233 double[][] eigenValues; field in NMRNoeMatrix 545 eigenValues[i][z] = relaxMatrix[i][z]; in Diagonalise() 592 double d = (eigenValues[p][p] - eigenValues[q][q]) in rotate() 593 / (2.0 * eigenValues[p][q]); in rotate() 597 eigenValues[p][p] += t * eigenValues[p][q]; in rotate() 598 eigenValues[q][q] -= t * eigenValues[p][q]; in rotate() 599 eigenValues[p][q] = eigenValues[q][p] = 0.0; in rotate() 602 double akp = c * eigenValues[k][p] + s * eigenValues[k][q]; in rotate() 603 double akq = -s * eigenValues[k][p] + c * eigenValues[k][q]; in rotate() 604 eigenValues[k][p] = eigenValues[p][k] = akp; in rotate() [all …]
|
/dports/misc/openvdb/openvdb-9.0.0/openvdb/openvdb/math/ |
H A D | Maps.cc | 106 Vec3d eigenValues; in createSymmetricMap() local 109 bool converged = math::diagonalizeSymmetricMatrix(m, Umatrix, eigenValues); in createSymmetricMap() 115 ScaleMap diagonal(eigenValues); in createSymmetricMap() 232 Vec3d eigenValues; in approxInverse() local 234 diagonalizeSymmetricMatrix(mat3, eigenVectors, eigenValues); in approxInverse() 238 if (std::abs(eigenValues[i]) < 10.0 * math::Tolerance<double>::value()) { in approxInverse() 241 d[i][i] = 1.f/eigenValues[i]; in approxInverse()
|