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

Development branch for SQN research (DO NOT MERGE) #17

Draft
wants to merge 54 commits into
base: master
Choose a base branch
from

Conversation

robbietuk
Copy link
Owner

@robbietuk robbietuk commented Jun 22, 2021

DO NOT MERGE!

At the time of making this PR, the merged changes are not yet in master but are needed for my STIR build.
The changes not in master but are here are (tick indicates it is in master):

  • Objective function prior Hessian methods and tests
  • Fix to the loglikelihood Hessian method and test
  • Fix to the loglikelihood gradient to actually compute y/ybar -1

robbietuk and others added 30 commits May 26, 2021 09:57
* RDP accumulate_Hessian_times_input

* Added RDP Hessian documentation

* Add div by 0 safety

* Add epsilon to Hessian function
robbietuk and others added 18 commits June 9, 2021 15:30
* RDP Documententation

* QP use (off_)diagonal_second_derivative methods and documentation

* Restructure RDP and QP Hessian logic

* Improve RDP and QP documentation

* Implement log cosh accumulate_Hessian_times_input in terms of derivatives

* Documentation
…sian_concavity

This is a change to keep consistency with changes in UCL#901
* Add method get_is_convex(), which accesses prior is_convex parameter and only test convexity of priors if convex function

* compute_Hessian method added to generalised prior and errors by default but with different messages depending on is_convex 

* Test rename and create test_Hessian_convexity and test_Hessian_methods

* Add test for compute_Hessian

* RelativeDifferencePrior initialisation call set_defaults

* Modernise `compute_Hessian` for QP and LogCosh. Add `compute_Hessian` for RDP

* Correct RDP second derivative functions

* Major changes to test_Hessian_against_numerical, which loops over each voxel for perturbation response

* Add verbosity suppression to suppress gradient info calls
@robbietuk robbietuk marked this pull request as draft June 22, 2021 20:01
@robbietuk robbietuk changed the title Development branch for SQN research Development branch for SQN research (DO NOT MERGE) Jun 22, 2021
@robbietuk robbietuk force-pushed the robbietuk/Devel/SQN branch from ffd3ff8 to adb294e Compare June 25, 2021 14:30
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.

1 participant