Bim package: Difference between revisions

Jump to navigation Jump to search
23 bytes removed ,  19 July 2012
no edit summary
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,delta,zeta);
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:length (uin), union (GammaD, GammaN)); ## INTERIOR NODES LIST
Omega = setdiff (1:nnodes, union (GammaD, GammaN)); ## INTERIOR NODES LIST
</pre>
</pre>


Line 161: Line 161:
</pre>
</pre>


<B> Save data for later visualization</B><BR>
<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]]  
43

edits

Navigation menu