-
Notifications
You must be signed in to change notification settings - Fork 120
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
CRHMC with Cooling Non-SphericalGaussians #330
CRHMC with Cooling Non-SphericalGaussians #330
Conversation
typedef typename NonSphericalGaussianFunctor::parameters<NT, Point> func_params; | ||
|
||
typedef crhmc_input<MT, Point, Func, Grad, ZeroFunctor<Point>> Input; | ||
//typedef crhmc_input<MT, Point, Func, Grad, Hess> Input; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we need this line?
|
||
Input input = convert2crhmc_input<Input, Polytope, Func, Grad, ZeroFunctor<Point>>(P, f, g, zerof); | ||
|
||
//Input input = convert2crhmc_input<Input, Polytope, Func, Grad, Hess>(P, f, g, h); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we need this line?
typedef typename NonSphericalGaussianFunctor::parameters<NT, Point> func_params; | ||
|
||
typedef crhmc_input<MT, Point, Func, Grad, ZeroFunctor<Point>> Input; | ||
//typedef crhmc_input<MT, Point, Func, Grad, Hess> Input; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we need this line?
func_params initial_f_params = func_params(Point(dimension), a_vals[0], -1, inv_covariance_matrix); | ||
Func initial_f(initial_f_params); | ||
Grad initial_g(initial_f_params); | ||
Hess initial_h(initial_f_params); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Since we use the zero functor for the hessian, do we need initial_h
?
|
||
Func f(f_params); | ||
Grad g(f_params); | ||
Hess h(f_params); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we use Hess h
?
|
||
Input input = convert2crhmc_input<Input, Polytope, Func, Grad, ZeroFunctor<Point>>(Pin_copy, f, g, zerof); | ||
|
||
//Input input = convert2crhmc_input<Input, Polytope, Func, Grad, Hess>(P, f, g, h); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we need this line?
|
||
//restore the new eta and start point, by looking at the walk after its operations | ||
initial_eta = walk.get_current_eta(); | ||
//start_point = walk.getPoint(); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Do we need this line?
typedef typename NonSphericalGaussianFunctor::parameters<NT, Point> func_params; | ||
|
||
typedef crhmc_input<MT, Point, Func, Grad, ZeroFunctor<Point>> Input; | ||
//typedef crhmc_input<MT, Point, Func, Grad, Hess> Input; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'd remove everything related to the Hessian, we should use the zero functor only.
|
||
//create the crhmc problem | ||
Input input = convert2crhmc_input<Input, Polytope, Func, Grad, ZeroFunctor<Point>>(newPin, f, g, zerof); | ||
//Input input = convert2crhmc_input<Input, Polytope, Func, Grad, Hess>(P, f, g, h); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Shall we delete this line?
7fc7895
to
3c3e904
Compare
Thanks @vgnecula! That's a crucial step towards a new implementation for volume estimation. |
No description provided.