43
edits
No edit summary |
No edit summary |
||
Line 95: | Line 95: | ||
<pre> | <pre> | ||
AdvDiff = bim2a_advection_diffusion (mesh, epsilon, 1, 1, phi); | AdvDiff = bim2a_advection_diffusion (mesh, epsilon, 1, 1, phi); | ||
Mass = bim2a_reaction (mesh, | Mass = bim2a_reaction (mesh, 1, 1); | ||
b = bim2a_rhs (mesh,f,g); | b = bim2a_rhs (mesh,f,g); | ||
A = AdvDiff + Mass; | A = AdvDiff + Mass; | ||
Line 112: | Line 112: | ||
jn = zeros (length (GammaN),1); ## PRESCRIBED NEUMANN FLUXES | jn = zeros (length (GammaN),1); ## PRESCRIBED NEUMANN FLUXES | ||
ud = 3*xu; ## DIRICHLET DATUM | ud = 3*xu; ## DIRICHLET DATUM | ||
Omega = setdiff (1: | Omega = setdiff (1:nnodes, union (GammaD, GammaN)); ## INTERIOR NODES LIST | ||
</pre> | </pre> | ||
Line 161: | Line 161: | ||
</pre> | </pre> | ||
<B> | <B> Export data to VTK format</B> | ||
The resut can be exported to vtk format to visualize with [[http://www.paraview.org|paraview]] | The resut can be exported to vtk format to visualize with [[http://www.paraview.org|paraview]] |
edits