— «Ты вот так просто всё понял?!»
— «Неужели пропасть между людьми и правда настолько огромна?!»
— Сунь Юй ведь был самым настоящим аспирантом прикладной математики. Он потратил целых два дня, изучая проект и сопутствующие материалы, но едва-едва уловил общую логику и ограничения существующих методов вычислений.
— «И этот „СЛ« на том конце провода... Неужели он действительно во всём разобрался, просто взглянув на файл и выслушав моё единственное объяснение?!»»
— «Быть того не может. Я выложил за этого эксперта целое состояние — две тысячи юаней. Если он пойдёт по ложному пути, предложенное им решение окажется совсем мимо кассы. Мои деньги просто вылетят в трубу!»
— При этой мысли Сунь Юй торопливо заговорил, изо всех сил стараясь звучать тактично: — Э-э, в этом проекте есть несколько довольно сложных мест. Пожалуй, мне стоит...
— Не успел он договорить, как голос в наушниках просто проигнорировал его и продолжил.
Типичный анализ устойчивости конструкций. Вы использовали стандартное для этой отрасли конечно-элементное программное обеспечение для расчета на устойчивость. Суть задачи — в решении обобщенной проблемы собственных значений: ([K]-λ[КГ]){φ}= 0.
— Голос «СЛ» звучал неторопливо.
— Сунь Юй, уже было открывший рот, так и замер с приоткрытыми губами.
— Потому что сказанное им оказалось совсем верным.
— Голос продолжил: — Судя по вашему описанию и этим из файла, трудность, с которой вы столкнулись, кроется в критическом коэффициенте запаса устойчивости λ. Именно его вам и нужно найти, и в теории он должен лежать в относительно четко оговоренном диапазоне. В файле указано значение 4,0 ± 0,3.
Но результаты многократных расчетов в вашей программе порой падают до 3,2, а иногда беспричинно подпрыгивают до 4,8. Результаты расходятся и никак не могут стабильно сойтись в ожидаемом теоретическом интервале.
— Сунь Юй остолбенел.
— «Он понял, что в этой роли полного клоуна выступал он сам».
— «Я, наверное, израсходовал всю удачу первой половины жизни, раз мне посчастливилось наткнуться на гения уровня „СЛ«. И что я учудил? Сначала проверял его навыки какой-то задачкой для бакалавров, а потом еще и сомневался, поймет ли он мой проект после одного-единственного объяснения!»»
— Сунь Юй поднял руку и отвесил себе легкую пощечину.
— Затем, словно воришка, он торопливо поднес мышку к экрану и нажал кнопку отключения микрофона.
— Было правдой то, что голос «СЛ» звучал молодо, напрочь лишенный той усталости, что неизменно приносит многолетняя работа в профессиональной сфере.
— Но теперь для Сунь Юя это спокойствие казалось почти пугающим.
— Он говорил о материале с полным безразличием, словно механически зачитывая текст с бумажки.
Итак, суть проблемы ясна. Ваш текущий решатель показывает чисшую неустойчивость при работе с кластеризацией собственных значений или плохо обусловленными матричными задачами. От меня вам требуется совсем новый алгоритм, который обладает лучшей численной устойчивостью и способен точно вычислить искомое наименьшее собственное значение.
— Стоило ему умолкнуть, как Сунь Юй заметил, что серая аватарка «СЛ» в окне конференции вдруг загорелась.
— Тот включил камеру.
— В груди Сунь Юя шевельнулось любопытство, и его взгляд метнулся к маленькому окошку видеосвязи.
— Однако на видео не оказалось лица гения, каким он его себе воображал.
— Камера была наклонена вниз и направлена прямо на чистый лист бумаги, лежавший на столе.
— Мгновение спустя в кадре появилась рука, сжимающая ручку.
— «СЛ» пояснял Сунь Юю свои действия, одновременно выводя формулы.
Я разработаю гибридный подход. Он состоит из нескольких шагов.
Во-первых, для вашей обобщенной проблемы собственных значений [K]x = λ[КГ]x нам необходимо построить умелый предобуславливатель [M], чтобы улучшить спектральные свойства матрицы коэффициентов.
— Кончик ручки заскользил по бумаге.
Здесь я предлагаю использовать технику предобуславливания на основе неполного разложения Холецкого или алгебраической многосеточной метода, применяемую напрямую к вашей матрице жесткости [K]. Поскольку, судя по вашим этим, матрица [K] обычно является хорошо обусловленной положительно определенной матрицей...
Во-вторых, мы задействуем алгоритм решателя, который более устойчив к плохо обусловленным задачам и кластеризации собственных значений.
Я рекомендую метод Якоби — Дэвидсона. По сравнению с тем, что вы используете сейчас, или более традиционными методами Ланцоша либо итераций подпространства, он в меньшей степени зависит от вектора начального приближения...
В-третьих, чтобы еще больше повысить точность и эффективность поиска минимального собственного значения, нам необходимо внедрить метод спектрального преобразования.
В рамках метода Якоби — Дэвидсона примените спектральное сдвиговое преобразование. Выберите сдвиг σ, близкий к вашему ожидаемому минимальному собственному значению λ_cr — в вашем случае это означает выбор σ = 4,0. Это превращает исходную задачу в стандартную проблему собственных значений, с которой гораздо проще справиться...
Наконец, адаптивная корректировка. Для решения упомянутой вами проблемы чувствительности сетки ключевые внутренние параметры алгоритма — такие как допуск сходимости и размерность подпространства — должны адаптироваться на основе оценок локального числа обусловленности матрицы, а не оставаться фиксированными.
В то же время по ходу работы алгоритма вам необходимо в режиме реального времени отслеживать скорость изменения нормы невязки и оценку собственного значения. Как только обнаруживается стагнация сходимости, предобуславливатель следует динамически обновлять...
— К тому моменту, когда он закончил, лист формата А4 уже был сплошь исписан густой вязью переменных и формул.
— Глаза Сунь Юя прилипли к экрану, но мозг уже давно перестал воспринимать информацию.
— Вдруг Сунь Юй вздрогнул.
— «Кто я? Где я?»
— Ему казалось, что многие из терминов и понятий, о которых говорил этот парень, он уже видел или слышал в аспирантуре, но стоило сложить их воедино, как они превращались для него в бессвязный бред.
— Впрочем, всё это уже не имело никакого значения.
— Теперь он был уверен: этот «СЛ» обладал настолько заоблачным уровнем математических познаний, который полностью превосходил его собственное понимание. Ему оставалось лишь благоговейно смотреть на него снизу вверх.
— И этого было достаточно!
— «Мне нужно лишь уцепиться за этого гиганта, и моя жизнь удалась!»
— В этот момент видеоокно погасло. «СЛ» отключил камеру.
— Он продолжил: — Основные шаги алгоритма и ключевые формулы расписаны на том листе бумаги. Что вам теперь нужно сделать, так это не торопясь досконально во всем разобраться, а затем перевести эти шаги в псевдокод. После этого найдите в своей компании программиста, который поможет воплотить его в виде скрипта, и наконец внедрите этот скрипт в свое программное обеспечение для конечных элементов.
— В голове у Сунь Юя зашумело. «Чего? Я? Да я против него — как монах Тан против целой толпы демонов!»
— Сунь Юй поспешно сказал: — Нет-нет-нет, пожалуйста, не могли бы вы сами написать для меня псевдокод? Как видите, наше время еще не вышло.
— Боясь получить отказ, он тут же добавил: — Послушайте, только взгляните на часы! Наша консультация началась около десяти, а еще даже нет одиннадцати. Трехчасовой тариф... У нас в запасе еще уйма времени!
— Теперь Сунь Юй считал, что единственный недостаток дорогих вещей — это их цена. Во всем остальном они были сплошным достоинством.
— На другом конце наушников воцарилось короткое молчание.
Ладно. Дайте мне минуту.
— Несколько минут спустя в интерфейсе конференции всплыло уведомление об отправке файла.
— «СЛ» поделился двумя изображениями.
— Сунь Юй скачал их, чтобы взглянуть. На первом красовался тот самый лист с набросками алгоритма и ключевыми формулами, а на втором был записан свежесоставленный псевдокод.
— С души Сунь Юя словно свалилась тяжелая плита.
— Он взглянул на часы. Было уже половина двенадцатого.
— Он быстро проговорил в голосовой чат: — Ладно, учитель, я закругляюсь и берусь за дело. Поздно уже, вам бы отдохнуть.
— Молодой голос на том ответил: — Хорошо, я отключаюсь. Время я записал; у нас осталось еще восемьдесят три минуты. Для следующей консультации свяжитесь со мной за три дня до встречи.
— Затем Сунь Юй увидел, как «СЛ» вышел из конференции.

Комментарии

Загрузка...