← Back to branch summary

~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to sandbox/bdm/main.cpp

  • Committer: Niclas Jansson
  • Date: 2011-06-10 14:33:43 UTC
  • Revision ID: njansson@csc.kth.se-20110610143343-d21p4am8rghiojfm
Added rudimentary header to binary files

Show diffs side-by-side

added added

removed removed

Lines of Context:
sandbox/bdm/main.cpp
1
 
// Copyright (C) 2007 Anders Logg and Marie Rognes
2
 
// Licensed under the GNU LGPL Version 2.1.
3
 
//
4
 
// First added:  2007-04-20
5
 
// Last changed: 2007-04-20
6
 
//
7
 
// This demo program solves the mixed formulation of
8
 
// Poisson's equation:
9
 
//
10
 
//     q + grad(u) = 0
11
 
//          div(q) = f
12
 
//
13
 
// The corresponding weak (variational problem)
14
 
//
15
 
//     <v, q> - <div(v), u> = 0       for all v
16
 
//               <w, div q> = <w, f>  for all w
17
 
//
18
 
// is solved using BDM (Brezzi-Douglas-Marini) elements
19
 
// of degree 1 (v, q) and DG (discontinuous Galerkin)
20
 
// elements of degree 0 for (w, u).
21
 
 
22
 
#include <dolfin.h>
23
 
#include "MixedPoisson.h"
24
 
 
25
 
using namespace dolfin;
26
 
 
27
 
int main()
28
 
{
29
 
  UnitSquare mesh(1, 1);
30
 
  Matrix A;
31
 
  MixedPoissonBilinearForm a;
32
 
 
33
 
  // Assemble matrix
34
 
  assemble(A, a, mesh);
35
 
  
36
 
  // Display matrix values
37
 
  A.disp();
38
 
 
39
 
  return 0;
40
 
}