GHAで使う非線形最適化真面目にプログラムの精度検証を始めたら、ニュートンガウス法(仮称)のロバスト性の低さが目についてきました。恐らく初期値が悪いのも大きな要因です。計算量が多いとは言え、何だかんだで平田のアルゴリズムは頑健なのでしょう。