Mercurial > hg > kwantix
view results/plotresults_vtk.m @ 24:ad9e3d28ce9b
Implement op<< for bvp::interpolator
author | Jordi Gutiérrez Hermoso <jordigh@gmail.com> |
---|---|
date | Fri, 21 Aug 2009 16:19:58 -0500 |
parents | 4c5bac1f2612 |
children |
line wrap: on
line source
skip = 1; for i = skip:skip:500 name = "u"; if(i < 10000) name = strcat(name, "0"); endif if(i < 1000) name = strcat(name, "0"); endif if(i < 100) name = strcat(name, "0"); endif if(i < 10) name = strcat(name, "0"); endif name = strcat(name,num2str(i)); #vtk_xlim([-1,1]); #vtk_ylim([-1,1]); #vtk_zlim([0,0.2]); filename = strcat(name,".map"); load(filename); h = eval(name); T = delaunay(h(:,1),h(:,2)); vtk_clear; vtk_trisurf(T,h(:,1),h(:,2),h(:,3)); ##foo = vtk_trimesh(T,h(:,1),h(:,2),h(:,3)); vtk_axis('off'); vtk_print(strcat(name,".png"),"-dpng"); clear eval(name); endfor