Синтаксис

[cgxu,cgyu] = pdecgrad(p,t,c,u)
[cgxu,cgyu] = pdecgrad(p,t,c,u,time)
[cgxu,cgyu] = pdecgrad(p,t,c,u,time,sdl)

Описание

[cgxu,cgyu] = pdecgrad(p,t,c,u) возвращает поток, $$c\otimes \nabla u$$, оцененный в центре каждого треугольника.

Строка i cgxu содержит

Строка i cgyu содержит

Существует один столбец для каждого треугольника в t и в cgxu и в cgyu.

Градиент, вычисленный pdegrad, является на самом деле тем же самым везде в треугольной внутренней части, потому что pdegrad использует только линейные основные функции. Контуры треугольников являются особым случаем: здесь производные могут быть прерывистыми. Однако поток $$c\otimes \nabla u$$ может отличаться в треугольнике, потому что коэффициент c может отличаться.

Геометрия проблемы УЧП дана данными о mesh p и t. Детали о представлении данных mesh могут быть найдены в записи на initmesh.

Коэффициент c проблемы УЧП может быть дан во множестве путей. Смотрите Коэффициенты УЧП.

Скалярный time дополнительного аргумента используется для параболических и гиперболических проблем, если c зависит от t, время.

Дополнительный аргумент sdl ограничивает вычисление субдоменами в списке sdl.

Смотрите также

evaluateCGradient | evaluateGradient