TODO in future release

[Buyun-Liang]

1. For avoiding unnecessary AD gradients, the GRANSO line search could be modified so that the AD gradient is only computed on demand, i.e., only compute the gradient when sufficient decrease holds

2. Adding the MPI parallel full BFGS Python code Tim has

[Buyun-Liang]

3. gradually switch to the vector_norm and matrix_norm functions

[Buyun-Liang]

https://pytorch.org/docs/stable/generated/torch.norm.html