2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Представление сумм в качестве Expression в FEniCS (dolfin)
Сообщение24.03.2016, 01:53 
Есть решение ДУсЧП, представленное в виде ряда. Хочу оценить погрешность решения, полученного МКЭ. Это можно сделать с помощью кода
Код:
errornorm(u, uh, norm_type='l2')

где uh -- МКЭ-решение, а u -- точное решение. В качестве u может быть подставлено либо Function, либо Expression. Не понятно, как создать Expression, отвечающий конечной сумме с произвольным количеством слагаемых.
Были такие идеи:
    1)Создать для каждого слагаемого по Expression'у и суммировать их. Однако, в результате получается объект класса Sum, который не принимается функцией errornorm.
    2)Создать строку, в которой будет записана искомая сумма в синтаксисе C++, и создать по ней Expression. Однако, такой способ оказывается неоправданно затратным по памяти.
Неужели в таком мощном инструменте нет возможности задавать функции, представимые суммой? :|

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group