[Getfem-commits] r4948 - /trunk/getfem/interface/src/scilab/demos/demo_slices.sce

[tomas . ligursky]

Author: ligut2am
Date: Fri Apr 10 17:08:25 2015
New Revision: 4948

URL: http://svn.gna.org/viewcvs/getfem?rev=4948&view=rev
Log:
removing an old demo

Removed:
    trunk/getfem/interface/src/scilab/demos/demo_slices.sce

Removed: trunk/getfem/interface/src/scilab/demos/demo_slices.sce
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/src/scilab/demos/demo_slices.sce?rev=4947&view=auto
==============================================================================
--- trunk/getfem/interface/src/scilab/demos/demo_slices.sce     (original)
+++ trunk/getfem/interface/src/scilab/demos/demo_slices.sce     (removed)
@@ -1,57 +0,0 @@
-// not working, not part of the getfem-interface distrib
-
-if getos()=='Windows' then
-  // Under Windows, all the trace messages are available in the dos console
-  // Under Linuxs, all the trace messages are redirected to the Scilab console
-  consolebox('on');
-end
-gf_util('trace level',3);
-gf_util('warning level',3);
-
-[mf] = gfMeshFem('load','signorini_cou.mesh_fem'); m = mf.linked_mesh;
-load signorini_cou.data; U=signorini_cou';
-
-mfdu = gf_mesh_fem(m,1);
-// the P2 fem is not derivable across elements, hence we use a discontinuous
-// fem for the derivative of U.
-gf_mesh_fem_set(mfdu,'fem',gf_fem('FEM_PRODUCT(FEM_PRODUCT(FEM_PK_DISCONTINUOUS(1,1),FEM_PK_DISCONTINUOUS(1,1)),FEM_PK_DISCONTINUOUS(1,1))'));
-
-// on output size(DU)=[3,3,nbdof(mfdu)]
-DU = gf_compute(mf,U,'gradient',mfdu);
-
-// from the derivative, we compute the von mises stress
-VM = zeros(1,gf_mesh_fem_get(mfdu,'nbdof'));
-N  = gf_mesh_get(m,'dim');
-
-for i=1:size(DU,3),
-  t = DU(:,:,i);
-  E = (t+t')/2;
-  VM(i) = sum(E(:).^2) - (1./N)*sum(diag(E))^2;
-end
-
-lambda = 1;
-VM = 4*lambda^2*VM;
-
-nrefine = 6;
-sl1 = gf_slice(list('boundary',list('none')),m,nrefine);
-c   = [0.1;0.1;20.1];
-x = [1;0;0];
-y = [0;1;0];
-z = [0;0;1];
-sl2 = 
gf_slice(list('boundary',list('union',list('planar',+1,c,x),list('planar',+1,c,y),list('planar',+1,c,z))),m,nrefine);
-//sl2 = 
gf_slice(list('boundary',list('union',list('planar',+1,c,x),list('planar',+1,c,y))),m,nrefine);
-
-P  = gf_slice_get(sl2,'pts'); 
-dP = gf_compute(mf,U,'interpolate on',sl2);
-gf_slice_set(sl2, 'pts', P+dP);
-
-VMsl = gf_compute(mfdu,VM,'interpolate on',sl2);
-scf(1);
-h = gf_plot_slice(sl2,'mesh','on','data',VMsl);
-//view(-80,-15); axis off; camlight;
-
-scf(2);
-h = gf_plot_slice(sl1,'mesh_faces','on','mesh','on');
-//view(-85,-15); axis off; camlight; 
-//set(h,'facecolor',[.8 0 0]);
-