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

Tracking pseudotime and interface for PseudoSteppingSchemes #240

Merged
merged 31 commits into from
Apr 11, 2022

Conversation

glwagner
Copy link
Member

@glwagner glwagner commented Apr 8, 2022

@sandreza you might be interested in this.

@glwagner
Copy link
Member Author

glwagner commented Apr 9, 2022

@adelinehillier we're getting

[ Info:   Testing EnsembleKalmanInversion construcation and basic iteration (Nbatch = 2)...
EnsembleKalmanInversions construction and iteration tests (Nbatch = 2): Error During Test at /Users/gregorywagner/Projects/ParameterEstimocean.jl/test/test_ensemble_kalman_inversion.jl:71
  Got exception outside of a @test
  PosDefException: matrix is not Hermitian; Cholesky factorization failed.
  Stacktrace:
    [1] checkpositivedefinite(info::Int64)
      @ LinearAlgebra /Users/julia/buildbot/worker/package_macos64/build/usr/share/julia/stdlib/v1.6/LinearAlgebra/src/factorization.jl:18
    [2] cholesky!(A::Matrix{Float64}, ::Val{false}; check::Bool)
      @ LinearAlgebra /Users/julia/buildbot/worker/package_macos64/build/usr/share/julia/stdlib/v1.6/LinearAlgebra/src/cholesky.jl:282
    [3] #cholesky#134
      @ /Users/julia/buildbot/worker/package_macos64/build/usr/share/julia/stdlib/v1.6/LinearAlgebra/src/cholesky.jl:378 [inlined]
    [4] cholesky (repeats 2 times)
      @ /Users/julia/buildbot/worker/package_macos64/build/usr/share/julia/stdlib/v1.6/LinearAlgebra/src/cholesky.jl:378 [inlined]

But this PR is supposed to be a refactor which doesn't affect functionality, so I'm not sure what happened.

@adelinehillier
Copy link
Collaborator

But this PR is supposed to be a refactor which doesn't affect functionality, so I'm not sure what happened.

Very strange! I'll look into it.

@@ -139,7 +139,6 @@ function EnsembleKalmanInversion(inverse_problem;
Xᵢ = unconstrained_parameters
iteration = 0
pseudotime = 0.0
pseudo_Δt = 0.0
Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thank you!!

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should we use default 1.0?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think 1.0 does make the most sense

Copy link
Contributor

@sandreza sandreza left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

lgtm

@glwagner glwagner merged commit df2d4d5 into main Apr 11, 2022
@glwagner glwagner deleted the glw/pseudotime branch April 11, 2022 18:24
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

Successfully merging this pull request may close these issues.

3 participants