boorlakov / fengine Goto Github PK

Описание проблемы
При запуске простой наиболее простой и возможной задачи, ответ не сходится с ожидаемым.

Чтобы воспроизвести
Шаги, чтобы воспроизвести ошибки:

1. Откройте тестовый проект FiniteElementsMethodTests
2. Откройте файл FemSolverTest
3. Наведитесь на соответствующий тест FemSolverWithSimpleIterationTest_WhenPassSimpleFuncAndNonUniformGrid_ShouldReturnCorrectResult
4. Запустите в режиме Debug
5. Смотрите на полученный результат и сравните его с ожидаемым

Ожидаемое поведение
Простой массив значений {2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0}

Screenshots

Desktop:

• OS: [MacOS]
• IDE: [Rider]
• Version: [2021.3.3]
• .NET Version: 6.0.1

Дополнительная информация
Возможно, ответ не сходится из-за неправильного формирования функции uK, которая впоследствии интегрируется. Возможно требуется передавать кусочно заданную функцию, а не всю, как это есть сейчас. Задача тяжелая, требуется расследование

Сделать это можно с помощью этого, например.

X axis - ось итераций
Y axis - ось значений невязки

Для того, чтобы код был более расширяемым и меньше зависил от конкретных реализаций интегрирования/сборок СЛАУ можно применить DI

https://github.com/boorlakov/umf/blob/a926d97b0af855ebf712a30e4bb4b0fec2e69710/umf-1/Fdm.cs#L4

Необходимо добавить коэфициент релаксации для ускорения сходимости простой итерации.

Коэф. релаксации должен уметь вводиться напрямую пользователем и уметь вычисляться автоматически как минимизация функционала.

• Учет коэф. релаксации в методе (backend)
• Отображение на стороне UI как дополительное поле с параметрами (автоматическое вычисление, или задается пользователем) (frontend)

Реализовать и протестировать метод Ньютона как на бекенде, так и учесть опцию на фронтенде