Lars Ruthotto

thanks for this, @stevenjgilmore ! Actually adding bound constraints makes a lot of sense here. I now add ``` lb = -4.*ones(2) ub = 4.*ones(2) lower_bounds!(opt,lb) upper_bounds!(opt,ub) ``` to my...

As @swufung said, increasing the flexibility of the package is a top priority of ours. There are a number of ways how to benefit more from the Julia ecosystem. Personally,...

Thanks for these references. I will read them with interest later. Again, I wouldn't be against trying adjoints, multiple shooting, and the other tools you've built in this package. Thanks...

Thanks. This sounds like a good thrust. Our Julia implementation here may be useful as a starting guess: https://github.com/EmoryMLIP/MFGnet.jl I won't have time to do much coding myself in the...

thanks for this excellent job. Can you please resolve the conflict and clean up the .m~ files?

Hey Hari, I resolved the conflicts and added three problems. The one concerning the documentation are critical. See if you also want to address the one regarding the performance.

Hi Hari, I have one comment before we can merge this (already did two minor changes to this). My goal is to finish this soon, but it's important to get...

Hi Hari! Thanks for your hard work and patience. It's getting close I think. One final thing that I see is to make the function getDiag in FEMPIREobjFctn work for...

Thanks, Hari. This works for me! Great job !

Do you want me to merge this right now, or should I wait for performance improvements? I've also noticed that the matrix free code does not run for non mass-preserving...