Skip to content

Commit

Permalink
remove offsets in newton
Browse files Browse the repository at this point in the history
see #232
  • Loading branch information
adtzlr committed Aug 6, 2022
1 parent 12d6b4f commit 665f72f
Showing 1 changed file with 3 additions and 8 deletions.
11 changes: 3 additions & 8 deletions felupe/tools/_newton.py
Original file line number Diff line number Diff line change
Expand Up @@ -123,7 +123,7 @@ def solve(A, b, x, dof1, dof0, offsets=None, ext0=None, solver=spsolve):
system = fesolve.partition(x, A, dof1, dof0, -b)
dx = fesolve.solve(*system, ext0, solver=solver)

return np.split(dx, offsets)
return dx


def check(dx, x, f, tol):
Expand Down Expand Up @@ -221,14 +221,9 @@ def newtonrhapson(

# solve linear system and update solution
sig = inspect.signature(solve)

try:
offsets = x.offsets
except:
offsets = []

keys = ["x", "dof1", "dof0", "offsets", "ext0", "solver"]
values = [x, dof1, dof0, offsets, ext0, solver]
keys = ["x", "dof1", "dof0", "ext0", "solver"]
values = [x, dof1, dof0, ext0, solver]

for key, value in zip(keys, values):

Expand Down

0 comments on commit 665f72f

Please sign in to comment.