Bim package: Difference between revisions

Bim package (edit)

Revision as of 22:38, 19 July 2012

48 bytes added , 19 July 2012

no edit summary

Carlo

43

edits

@@ Line 106: / Line 106: @@
 <pre>
-GammaD = bim2c_unknowns_on_side (mesh, [1 2]); 	          ## DIRICHLET NODES LIST
+GammaD = bim2c_unknowns_on_side (mesh, [1 2]); 	    ## DIRICHLET NODES LIST
-GammaN = bim2c_unknowns_on_side (mesh, [3 4]); 	          ## NEUMANN NODES LIST
+GammaN = bim2c_unknowns_on_side (mesh, [3 4]); 	    ## NEUMANN NODES LIST
 GammaN = setdiff (GammaN, GammaD);
-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:nnodes, union (GammaD, GammaN)); ## INTERIOR NODES LIST
 </pre>
@@ Line 125: / Line 125: @@
 Ani = A(GammaN, Omega);
-Aid = A(Ilist, GammaD);
+Aid = A(Omega, GammaD);
-Ain = A(Ilist, GammaN);
+Ain = A(Omega, GammaN);
-Aii = A(Ilist, Omega);
+Aii = A(Omega, Omega);
 bd = b(GammaD);
@@ Line 137: / Line 137: @@
 <pre>
-temp = [Ann Ani ; Ain Aii ] \ [ Fn+bn-And*ud(GammaD) ; bi-Aid*ud(GammaD)];
+temp = [Ann Ani ; Ain Aii ] \ [ jn+bn-And*ud(GammaD) ; bi-Aid*ud(GammaD)];
-u(Nlist)   = temp(1:numel (GammaN));
+u = ud;
-u(Omega)   = temp(length(un)+1:end);
+u(GammaN)  = temp(1:numel (GammaN));
+u(Omega)   = temp(length(GammaN)+1:end);
 </pre>
@@ Line 145: / Line 146: @@
 <pre>
-Fd = [Add Adi Adn] * u([GammaD; Omega; GammaN]) - bd;
+jd = [Add Adi Adn] * u([GammaD; Omega; GammaN]) - bd;
 </pre>
@@ Line 158: / Line 159: @@
 <pre>
-[jxglob, jyglob] = bim2c_global_flux (mesh, u, epsilon, 1, phi);
+[jxglob, jyglob] = bim2c_global_flux (mesh, u, epsilon*ones(nelems, 1), ones(nnodes, 1), ones(nnodes, 1), phi);
 </pre>
@@ Line 166: / Line 167: @@
 or [[https://wci.llnl.gov/codes/visit/|visit]]
 <pre>
-fpl_vtk_write_field ("vtkdata", mesh, {}, {[gx; gy]', "Gradient"}, 1);
+fpl_vtk_write_field ("vtkdata", mesh, {u, "Solution"}, {[gx; gy]', "Gradient"}, 1);
 </pre>
 [[Category:OctaveForge]]
 [[Category:Packages]]