/*-------- Third Party --------*/
#include <deal.II/numerics/vector_tools.h>
-#include "braid.h"
-#include "braid_test.h"
+#include <braid.h>
+#include <braid_test.h>
/*-------- Project --------*/
#include "HeatEquation.hh"
#ifndef _HEATEQUATION_H_
#define _HEATEQUATION_H_
-#include <fstream>
-
#include <deal.II/base/utilities.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/function.h>
#include <deal.II/lac/sparse_matrix.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/precondition.h>
-#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/lac/affine_constraints.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_refinement.h>
#include <deal.II/numerics/matrix_tools.h>
#include <deal.II/base/convergence_table.h>
+#include <fstream>
+
using namespace dealii;
// The HeatEquation class is describes the finite element
FE_Q<dim> fe;
DoFHandler<dim> dof_handler;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
SparsityPattern sparsity_pattern;
SparseMatrix<double> mass_matrix;
Vector<double> system_rhs;
- std::ofstream myfile;
+ std::ofstream myfile;
const double theta;
* Author: Joshua Christopher, Colorado State University, 2018
*/
-#include <fstream>
-#include <iostream>
-
#include "BraidFuncs.hh"
#include "HeatEquation.hh"
#include "Utilities.hh"
+#include <fstream>
+#include <iostream>
+
int main(int argc, char *argv[])
{
try
braid_Core core;
double tstart = 0.0;
double tstop = 0.002;
- double ntime = 10;
+ int ntime = 10;
my_App *app = new(my_App);
braid_Init(MPI_COMM_WORLD, comm, tstart, tstop, ntime, app,