multifrontal_schur_solve

Source: tests/multifrontal/detail/multifrontal_schur_solve.cpp

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

// Containers.
#include <myramath/sparse/Pattern.h>
#include <myramath/sparse/PatternRange.h>
#include <myramath/sparse/PatternBuilder.h>
#include <myramath/multifrontal/symbolic/AssemblyTree.h>
#include <myramath/multifrontal/symbolic/SchurTree.h>

// Algorithms.
#include <myramath/sparse/laplacian3.h>

// JobGraph under test.
#include <myramath/jobgraph/verify.h>
#include <myramath/jobgraph/graphviz.h>
#include <myramath/multifrontal/detail/schursolve.h>

// Reporting.
#include <tests/myratest.h>

using namespace myra;
typedef multifrontal::Options Options;

ADD_TEST("multifrontal_schur_solve","[multifrontal][jobgraph]")
  {
  // Construct AssemblyTree of a 3D laplacian.
  int I = 24;
  int J = 24;
  int K = 24;  
  int N = I*J*K;
  Pattern A_pattern = stencil3(I,J,K);
  Permutation perm = bisect3(I,J,K);
  AssemblyTree atree(A_pattern,perm);
  // Construct SchurTree for k=0 plane.
  Natural3D n(I,J,K);
  PatternBuilder B_builder(N,I*J);
  int next = 0;
  for (int i = 0; i < I; ++i)
    for (int j = 0; j < J; ++j)
      B_builder.insert(n(i,j,0),next++);
  Pattern B_pattern = B_builder.make_Pattern();
  SchurTree stree(atree,B_pattern);
  // Verify the JobGraph.
  typedef ::multifrontal::detail::schursolve::JobGraphBase1 Graph;
  Graph graph(&stree);
  REQUIRE( verify(graph) );
  // Save .dot file to disk?
  graphviz(graph,"graph2.dot");  
  }

Results: [PASS]

Go back to Summary of /test programs.