Тестирование
Работоспособность и эффективность построенного алгоритма переменного шага исследовалась на решении линейных задач из [62]
где A и B- квадратные матрицы размерности N ? Nсоответственно с элементами
- N-мерный вектор с постоянными
коэффициентами, v(t) - скалярная функция.
Во всех приведенных ниже примерах
; матрица Aдиагональная, а матрица B трехдиагональная с элементами
Функция v (t) определяется по формуле
Тестовая задача, аналитически приведенная к разрешенному виду, представлена в [78]. Решение разрешенной задачи выполнялось алгоритмом на основе одностадийной формулы типа Розенброка, описанной в [58]. Таким образом, тестовые задачи решены в классической и неявной форме.
105
Из результатов расчетов следует, что вычислительные затраты для алгоритмов решения разрешенных и не разрешенных относительно производной задач практически совпадают. Однако в специальной численной формуле нет необходимости в предварительном разрешении системы (1.5), что предпочтительнее в случае анализа нелинейных задач.
4.2