/[ff3d]/ff3d/testsuite/laplacian-lens.pl
ViewVC logotype

Diff of /ff3d/testsuite/laplacian-lens.pl

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.2 by ArceneL, Fri Oct 24 12:30:31 2003 UTC revision 1.3 by ArceneL, Fri Nov 21 10:30:14 2003 UTC
# Line 23  Line 23 
23  # that it defines boundary conditions on pieces of objects.  # that it defines boundary conditions on pieces of objects.
24    
25  require 'check.ph';  require 'check.ph';
26    CleanUp('laplacian-lens.bb');
27    
28  # The Povray scene  # The Povray scene
29  WriteFile('test.pov',<< 'EOF');  WriteFile("laplacian-lens.pov",<< 'EOF');
30    
31  /* We define two spheres to be able to apply different boundary  /* We define two spheres to be able to apply different boundary
32  conditions on each */  conditions on each */
# Line 52  vector n = (16,16,16); Line 53  vector n = (16,16,16);
53  vector a = (-2,-2,-2);  vector a = (-2,-2,-2);
54  vector b = (3,2,2);  vector b = (3,2,2);
55    
56  scene S = pov("test.pov");  scene S = pov("laplacian-lens.pov");
57  mesh M = structured(n,a,b);  mesh M = structured(n,a,b);
58  domain O = domain(S,outside(<1,0,0>) or outside(<0,1,0>));  domain O = domain(S,outside(<1,0,0>) or outside(<0,1,0>));
59    
# Line 65  solve(u) in O by M{ Line 66  solve(u) in O by M{
66      u = 0 on M;      u = 0 on M;
67  };  };
68    
69  save(medit,"test",u,M);  save(medit,"laplacian-lens",u,M);
70    
71  // Also saves mesh, just for debugging purposes  // Also saves mesh, just for debugging purposes
72  save(medit,"test",M);  save(medit,"laplacian-lens",M);
73  EOF  EOF
74    
75  # Compares the resulting output with a reference  # Compares the resulting output with a reference
76  VectorCheck('test.bb','laplacian-lens-ref.bb',1e-6);  VectorCheck('laplacian-lens.bb','laplacian-lens-ref.bb',1e-6);
77    
78  Report();  Report();

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26