From: <cd...@us...> - 2010-03-26 17:24:54
|
Revision: 7131 http://octave.svn.sourceforge.net/octave/?rev=7131&view=rev Author: cdf Date: 2010-03-26 17:24:48 +0000 (Fri, 26 Mar 2010) Log Message: ----------- tutorial Modified Paths: -------------- trunk/octave-forge/extra/bim/DESCRIPTION trunk/octave-forge/extra/bim/doc/tutorial.html Modified: trunk/octave-forge/extra/bim/DESCRIPTION =================================================================== --- trunk/octave-forge/extra/bim/DESCRIPTION 2010-03-26 16:34:55 UTC (rev 7130) +++ trunk/octave-forge/extra/bim/DESCRIPTION 2010-03-26 17:24:48 UTC (rev 7131) @@ -5,7 +5,7 @@ Maintainer: Culpo Massimiliano Title: PDE Solver using a Finite Element/Finite Volume approach Description: Package for solving Diffusion Advection Reaction (DAR) Partial Differential Equations based on the Finite Volume Scharfetter-Gummel (FVSG) method a.k.a Box Integration Method (BIM) -Depends: octave (>= 3.0.0), fpl, msh +Depends: octave (>= 3.2.0), fpl, msh Autoload: no License: GPL version 2 or later SVNRelease: 6866 \ No newline at end of file Modified: trunk/octave-forge/extra/bim/doc/tutorial.html =================================================================== --- trunk/octave-forge/extra/bim/doc/tutorial.html 2010-03-26 16:34:55 UTC (rev 7130) +++ trunk/octave-forge/extra/bim/doc/tutorial.html 2010-03-26 17:24:48 UTC (rev 7131) @@ -23,8 +23,8 @@ For a linear problem only the values at boundary nodes are actually relevant<br> <p><pre> -xu = mesh.p(1,:); -yu = mesh.p(2,:); +xu = mesh.p(1,:).'; +yu = mesh.p(2,:).'; nelems = columns(mesh.t); nnodes = columns(mesh.p); uin = 3*xu; @@ -93,6 +93,9 @@ temp = [Ann Ani ; Ain Aii ] \ [ Fn+bn-And*ud ; bi-Aid*ud]; un = temp(1:length(un)); ui = temp(length(un)+1:end); +u(Dlist) = ud; +u(Ilist) = ui; +u(Nlist) = un; </pre></p> <b> Compute the fluxes through Dirichlet sides</b><br> @@ -113,7 +116,8 @@ <B> Save data for later visualization</B><BR> <p><pre> -fpl_dx_write_field("dxdata",mesh,[gx gy],"Gradient",1,2,1); +fpl_dx_write_field("dxdata",mesh,[gx; gy]',"Gradient",1,2,1); +fpl_vtk_write_field ("vtkdata", mesh, {}, {[gx; gy]', "Gradient"}, 1); </pre></p> </body> This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |