Skip to content

Commit

Permalink
Another small fix
Browse files Browse the repository at this point in the history
  • Loading branch information
Rapo committed Aug 15, 2017
1 parent 32cf832 commit 090f649
Showing 1 changed file with 2 additions and 4 deletions.
6 changes: 2 additions & 4 deletions src/craig_bampton.jl
Expand Up @@ -10,15 +10,13 @@ r = retained DOF:s, l = internal DOF:s, n = the number of modes to keep.
"""
function craig_bampton(K, M, r, l, n)
eps_ = 0.000001
eps_ = 0.0000001
Krr = K[r,r]; Krl = K[r,l]
Klr = K[l,r]; Kll = K[l,l]
Mrr = M[r,r]; Mrl = M[r,l]
Mlr = M[l,r]; Mll = M[l,l]
if issparse(K) && issparse(M)
tmp = eigs(Kll, Mll)
w2 = tmp[1]
X1 = tmp[2]
w2, X1 = eigs(Kll, Mll)
else
w2 = eigvals(Kll,Mll)
X1 = eigvecs(Kll,Mll)
Expand Down

0 comments on commit 090f649

Please sign in to comment.