Какова формула вычисления остатка в fipy?

Question

Я знаю, что мы можем получить остаток, если использовать «развертку» для решения уравнения, но какова формула для вычисления остатка в fipy? В чем разница между двумя сторонами уравнения? И почему остатки моего уравнения уменьшаются очень медленно? Может кто-нибудь помочь мне? Спасибо!

Answer 1

Невязка, возвращаемая sweep(), является L2-нормой L*x - b (где L - матрица, x - вектор решения, а b - вектор правой части для система линейных уравнений, полученная в результате дискретизации вашего уравнения (а)). Я подал проблему № 734 , чтобы добавить документацию по этому поводу.

Что касается того, почему невязка уменьшается очень медленно, невозможно сказать, не зная, какие уравнения вы решаете и как вы решаете им.