Skip to content

Commit

Permalink
Euler-Bernoulli Eigenmodes correct with scipy.linalg.eigh
Browse files Browse the repository at this point in the history
  • Loading branch information
Veit Gufler committed Dec 15, 2020
1 parent afe3e40 commit ca3b869
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions EasyBeam/EasyBeam.py
Original file line number Diff line number Diff line change
Expand Up @@ -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])
Expand Down

0 comments on commit ca3b869

Please sign in to comment.