Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-10-02 22:31:31 +00:00
Code Issues Packages Projects Releases Wiki Activity

better heat example

Browse Source
This commit is contained in:
Jukka Aho
2016-06-09 23:17:57 +03:00
parent 06a4623acc
commit 74ed32a96c
2 changed files with 49 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/solvers.jl
+8
View File
@@ -312,14 +312,22 @@ function call(solver::Solver)
info("Starting nonlinear iteration #$(solver.iteration)")
# 2.1 update linearized assemblies (if needed)
info("Assembling problems ...")
tic()
for problem in solver.problems
problem.assembly.changed = true # force reassembly
assemble!(problem, solver.time)
end
t1 = round(toq(), 2)
info("Assembled in $t1 seconds.")
# 2.2 call solver for linearized system (default: direct lu factorization)
info("Solve linear system ...")
tic()
u, la = solve_linear_system(solver, Val{solver.linear_system_solver})
push!(solver.norms, (norm(u), norm(la)))
t1 = round(toq(), 2)
info("Solved Ax = b in $t1 seconds.")
# 2.3 update solution back to elements
for problem in solver.problems