From ca3b8696702b307e3764fedd71d3e9687adf40ea Mon Sep 17 00:00:00 2001 From: Veit Gufler Date: Tue, 15 Dec 2020 21:29:33 +0100 Subject: [PATCH] Euler-Bernoulli Eigenmodes correct with scipy.linalg.eigh --- EasyBeam/EasyBeam.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/EasyBeam/EasyBeam.py b/EasyBeam/EasyBeam.py index 44054bd..f58f13e 100755 --- a/EasyBeam/EasyBeam.py +++ b/EasyBeam/EasyBeam.py @@ -185,10 +185,10 @@ def EigenvalueAnalysis(self, nEig=2, massMatrixType="consistent"): self.k = self.Assemble(self.StiffMatElem) self.m = self.Assemble(self.MassMatElem) try: - lambdaComplex, self.Phi = npla.eigh(self.k[self.DoF, :][:, self.DoF], + lambdaComplex, self.Phi = spla.eigh(self.k[self.DoF, :][:, self.DoF], self.m[self.DoF, :][:, self.DoF], eigvals=(0, nEig-1)) - self.EigenvalSolver = "numpy.linalg.eigh" + self.EigenvalSolver = "scipy.linalg.eigh" except: lambdaComplex, self.Phi = spla.eig(self.k[self.DoF, :][:, self.DoF], self.m[self.DoF, :][:, self.DoF])