[PATCH] add extensive test for igc updates

[Jeffrey D. Oldham]

Richard Guenther wrote:

>With igc update optimizations it is important to check if they work
>correctly in all cases, so here is a testcase that (tires to)
>enumerate all possible cases.
>
>Ok?
>
>  
>
Testing is good.  Yes, please add it to the Pooma repository.

>Richard.
>
>
>2004Aug18  Richard Guenther <richard.guenther at uni-tuebingen.de>
>
>	* src/Array/tests/array_test30.cpp: new.
>  
>
>------------------------------------------------------------------------
>
>--- /dev/null	Tue May 18 17:20:27 2004
>+++ array_test30.cpp	Wed Aug 18 13:47:46 2004
>@@ -0,0 +1,143 @@
>+// -*- C++ -*-
>+// ACL:license
>+// ----------------------------------------------------------------------
>+// This software and ancillary information (herein called "SOFTWARE")
>+// called POOMA (Parallel Object-Oriented Methods and Applications) is
>+// made available under the terms described here.  The SOFTWARE has been
>+// approved for release with associated LA-CC Number LA-CC-98-65.
>+// 
>+// Unless otherwise indicated, this SOFTWARE has been authored by an
>+// employee or employees of the University of California, operator of the
>+// Los Alamos National Laboratory under Contract No. W-7405-ENG-36 with
>+// the U.S. Department of Energy.  The U.S. Government has rights to use,
>+// reproduce, and distribute this SOFTWARE. The public may copy, distribute,
>+// prepare derivative works and publicly display this SOFTWARE without 
>+// charge, provided that this Notice and any statement of authorship are 
>+// reproduced on all copies.  Neither the Government nor the University 
>+// makes any warranty, express or implied, or assumes any liability or 
>+// responsibility for the use of this SOFTWARE.
>+// 
>+// If SOFTWARE is modified to produce derivative works, such modified
>+// SOFTWARE should be clearly marked, so as not to confuse it with the
>+// version available from LANL.
>+// 
>+// For more information about POOMA, send e-mail to pooma at acl.lanl.gov,
>+// or visit the POOMA web page at http://www.acl.lanl.gov/pooma/.
>+// ----------------------------------------------------------------------
>+// ACL:license
>+
>+//-----------------------------------------------------------------------------
>+// array_test30: verify correctness of igc updates
>+//-----------------------------------------------------------------------------
>+
>+// Include files
>+
>+#include "Pooma/Arrays.h"
>+#include "Utilities/Tester.h"
>+#include <iostream>
>+
>+
>+template <class A1, class A2>
>+bool test(Pooma::Tester& tester,
>+	  const A1& a_mp, const A1& b_mp,
>+	  const A2& a_sp, const A2& b_sp,
>+	  const Loc<2>& delta1, const Loc<2>& delta2,
>+	  bool initial_f, const Loc<2>& initial)
>+{
>+  static int sequence = 0;
>+  Interval<2> I;
>+
>+  // initialize rhs arrays, ensure wrong igc values
>+  // via sequence number.
>+  I = b_sp.totalDomain();
>+  b_sp(I) = sequence + iota(I).comp(0) + I[0].size()*iota(I).comp(1);
>+  b_mp.engine().setGuards(0);
>+  b_mp(I) = b_sp(I);
>+
>+  // if requested, force initial update of a set of igcs
>+  if (initial_f) {
>+    b_sp(b_sp.physicalDomain()) = b_mp(b_sp.physicalDomain()+initial);
>+    b_sp(I) = sequence + iota(I).comp(0) + I[0].size()*iota(I).comp(1);
>+    Pooma::blockAndEvaluate();
>+  }
>+
>+  // do calculation both sp and mp
>+  I = a_sp.physicalDomain();
>+  a_sp(I) = b_sp(I+delta1) - b_sp(I+delta2);
>+  a_mp(I) = b_mp(I+delta1) - b_mp(I+delta2);
>+
>+  // check the results are the same everywhere
>+  bool res = all(a_sp(I) == a_mp(I));
>+  tester.out() << "For deltas " << delta1 << " and " << delta2 << " ";
>+  if (initial_f)
>+    tester.out() << "with initial " << initial << " ";
>+  tester.check("result is", res);
>+  if (!res) {
>+    int n = b_mp.layout().sizeGlobal();
>+    for (int i=0; i<n; ++i) {
>+      Array<2, int, Remote<Brick> > b(b_mp.engine().globalPatch(i));
>+      tester.out() << "Brick " << i << " " << intersect(b.domain(), b_mp.physicalDomain())
>+	      	   << " on context " << b.engine().owningContext()
>+		   << " is\n" << b(intersect(b.totalDomain(), b_mp.physicalDomain()))
>+		   << std::endl;
>+    }
>+    tester.out() << "Aborting." << std::endl;
>+    return false;
>+  }
>+
>+  sequence++;
>+
>+  return true;
>+}
>+
>+
>+int main(int argc, char *argv[])
>+{
>+  // Initialize POOMA and output stream, using Tester class
>+  Pooma::initialize(argc, argv);
>+  Pooma::Tester tester(argc, argv);
>+
>+  Interval<2> domain(12, 12);
>+  UniformGridLayout<2> layout_mp(domain, Loc<2>(3, 3),
>+				 GuardLayers<2>(2), DistributedTag());
>+  DomainLayout<2> layout_sp(domain, GuardLayers<2>(2));
>+
>+  Array<2, int, MultiPatch<UniformTag, Remote<Brick> > >
>+    a_mp(layout_mp), b_mp(layout_mp);
>+  Array<2, int, Brick>
>+    a_sp(layout_sp), b_sp(layout_sp);
>+
>+  // all 5^4 == 625 uninitialized cases
>+  for (int d1i = -2; d1i <= 2; ++d1i)
>+    for (int d1j = -2; d1j <= 2; ++d1j)
>+      for (int d2i = -2; d2i <= 2; ++d2i)
>+	for (int d2j = -2; d2j <= 2; ++d2j)
>+	  if (!test(tester, a_mp, b_mp, a_sp, b_sp,
>+		    Loc<2>(d1i, d1j), Loc<2>(d2i, d2j),
>+		    false, Loc<2>(0)))
>+	    goto out;
>+
>+  // all 5^4 == 625 initialized cases with simplified expression
>+  for (int ii = -2; ii <= 2; ++ii)
>+    for (int ij = -2; ij <= 2; ++ij)
>+      for (int d1i = -2; d1i <= 2; ++d1i)
>+        for (int d1j = -2; d1j <= 2; ++d1j)
>+	  if (!test(tester, a_mp, b_mp, a_sp, b_sp,
>+		    Loc<2>(d1i, d1j), Loc<2>(d1i, d1j),
>+		    true, Loc<2>(ii, ij)))
>+	    goto out;
>+
>+ out:
>+  tester.out() << "Best testing is done with all 1 to 9 processes" << std::endl;
>+
>+  int retval = tester.results("array_test30");
>+  Pooma::finalize();
>+  return retval;
>+}
>+
>+// ACL:rcsinfo
>+// ----------------------------------------------------------------------
>+// $RCSfile: array_test29.cpp,v $   $Author: pooma $
>+// $Revision: 1.1 $   $Date: 2004/07/20 18:41:00 $
>+// ----------------------------------------------------------------------
>+// ACL:rcsinfo
>  
>


-- 
Jeffrey D. Oldham
oldham at codesourcery.com