Skip to content

Commit

Permalink
remove som annoying outputs
Browse files Browse the repository at this point in the history
  • Loading branch information
lijas committed Oct 27, 2023
1 parent 2bafc4f commit dd46c8b
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 5 deletions.
1 change: 0 additions & 1 deletion src/assembling.jl
Original file line number Diff line number Diff line change
Expand Up @@ -59,7 +59,6 @@ end

function assemble_stiffnessmatrix_and_forcevector!(dh, state::StateVariables, globaldata)
for (partid, part) in enumerate(globaldata.parts)
@info "Assembling for partid $(partid), $(typeof(part))"
assemble_stiffnessmatrix_and_forcevector!(dh, part, state)
end
end
Expand Down
4 changes: 0 additions & 4 deletions src/solvers/newton_solver.jl
Original file line number Diff line number Diff line change
Expand Up @@ -35,30 +35,26 @@ function step!(solver::NewtonSolver, state::StateVariables, globaldata, ntries=0
#Apply boundary conditions
update!(ch, state.t)
apply!(state.d, ch)

state.newton_itr = 0
while true
state.newton_itr += 1

fill!(state.system_arrays, 0.0)

@info "[NEWTONSOLVER] Assembling"
@timeit "Assembling" assemble_stiffnessmatrix_and_forcevector!(dh, state, globaldata)
@timeit "ExternalForces" apply_external_forces!(dh, globaldata.efh, state, globaldata)
@timeit "Apply constraint" apply_constraints!(dh, globaldata.constraints, state, globaldata)

r = state.system_arrays.fⁱ - state.system_arrays.fᵉ
K = state.system_arrays.Kⁱ #- state.system_arrays.Kᵉ
#Solve
@info "[NEWTONSOLVER] Applying"
apply_zero!(K, r, ch)

#Kt = ThreadedSparseMatrixCSC(K)' #Note the transpose
prob = LinearSolve.LinearProblem(K, -r)
precon = solver.preconditioner(K)
linsolve = LinearSolve.init(prob, solver.linearsolver, Pl=precon)

@info "[NEWTONSOLVER] Solving"
@timeit "Solve" sol = LinearSolve.solve(linsolve)
ΔΔd = sol.u
apply_zero!(ΔΔd, ch)
Expand Down

0 comments on commit dd46c8b

Please sign in to comment.