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

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

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

revision 1.2 by ArceneL, Fri Sep 26 07:55:59 2003 UTC revision 1.3 by ArceneL, Fri Nov 21 10:30:15 2003 UTC
# Line 22  Line 22 
22  # A Laplacian on a spiky geometry.  # A Laplacian on a spiky geometry.
23    
24  require 'check.ph';  require 'check.ph';
25    CleanUp('laplacian-spike.bb');
26    
27  # The Povray scene  # The Povray scene
28  WriteFile('test.pov',<< 'EOF');  WriteFile("laplacian-spike.pov",<< 'EOF');
29  /* The sphere surrounds the box in such a way that only one corner  /* The sphere surrounds the box in such a way that only one corner
30  still appears */  still appears */
31    
# Line 47  vector n = (20,20,20); Line 48  vector n = (20,20,20);
48  vector a = (-4,-4,-4);  vector a = (-4,-4,-4);
49  vector b = (2,2,2);  vector b = (2,2,2);
50    
51  scene S = pov("test.pov");  scene S = pov("laplacian-spike.pov");
52  mesh M = structured(n,a,b);  mesh M = structured(n,a,b);
53  domain O = domain(S,outside(<1,0,0>));  domain O = domain(S,outside(<1,0,0>));
54    
# Line 59  solve(u) in O by M{ Line 60  solve(u) in O by M{
60      u = 0 on M;      u = 0 on M;
61  };  };
62    
63  save(medit,"test",M);  save(medit,"laplacian-spike",M);
64  save(medit,"test",u,M);  save(medit,"laplacian-spike",u,M);
65  EOF  EOF
66    
67  # Compares the resulting output with a reference  # Compares the resulting output with a reference
68  VectorCheck('test.bb','laplacian-spike-ref.bb',1e-6);  VectorCheck('laplacian-spike.bb','laplacian-spike-ref.bb',1e-6);
69    
70  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