Few days ago I worked on gradinet descent method for optimized stitching. I wrote the dirty implementation on python + numpy. I've done some primitive tests and verified in C# little program, which I made for viewing points before and after optimization. I've discovered that algorithm converges, and even it requiers sometimes lots of time (but still admissible), it seems to me that I'm on right way.