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