Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

propagate verbose keyword to linsolve #1571

Open
vpuri3 opened this issue Jan 18, 2022 · 2 comments
Open

propagate verbose keyword to linsolve #1571

vpuri3 opened this issue Jan 18, 2022 · 2 comments

Comments

@vpuri3
Copy link
Member

vpuri3 commented Jan 18, 2022

if islinear(f)
du1 = rate_prototype
uf = nothing
jac_config = nothing
linprob = LinearProblem(W,_vec(k); u0=_vec(dz))
Pl,Pr = wrapprecs(alg.precs(W,nothing,u,p,t,nothing,nothing,nothing,nothing)...,weight)
linsolve = init(linprob,alg.linsolve,alias_A=true,alias_b=true,
Pl = Pl, Pr = Pr)
else
du1 = zero(rate_prototype)
if isdae
uf = DAEResidualJacobianWrapper(alg,f,p,α,inv*dt),tmp,uprev,t)
else
uf = build_uf(alg,nf,t,p,Val(true))
end
jac_config = build_jac_config(alg,nf,uf,du1,uprev,u,ztmp,dz)
linprob = LinearProblem(W,_vec(k); u0=_vec(dz))
Pl,Pr = wrapprecs(alg.precs(W,nothing,u,p,t,nothing,nothing,nothing,nothing)...,weight)
linsolve = init(linprob,alg.linsolve,alias_A=true,alias_b=true,
Pl = Pl, Pr = Pr)
end

@ChrisRackauckas
Copy link
Member

Hmm, that might be a bit much because verbose is true by default and meant only for warnings. We need different levels of verbosity.

@vpuri3
Copy link
Member Author

vpuri3 commented Jan 27, 2022

yeah it's a linear solve embedded in a nonlinear solve embedded in an integrator. we need verbosity keyword for each embedded solve i guess

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants