← Back to branch summary

~ubuntu-branches/ubuntu/trusty/rheolef/trusty

« back to all changes in this revision

Viewing changes to nfem/form_element/2D_s_D_s.cc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito
  • Date: 2012-04-06 09:12:21 UTC
  • mfrom: (1.1.5)
  • Revision ID: package-import@ubuntu.com-20120406091221-m58me99p1nxqui49
Tags: 6.0-1
http://bugs.debian.org/667356
http://bugs.debian.org/661689
* New upstream release 6.0 (major changes):
  - massively distributed and parallel support
  - full FEM characteristic method (Lagrange-Gakerkin method) support
  - enhanced users documentation 
  - source code supports g++-4.7 (closes: #667356)
* debian/control: dependencies for MPI distributed solvers added
* debian/rules: build commands simplified
* debian/librheolef-dev.install: man1/* to man9/* added
* debian/changelog: package description rewritted (closes: #661689)

Show diffs side-by-side

added added

removed removed

Lines of Context:
nfem/form_element/2D_s_D_s.cc
 
1
///
 
2
/// This file is part of Rheolef.
 
3
///
 
4
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
 
5
///
 
6
/// Rheolef is free software; you can redistribute it and/or modify
 
7
/// it under the terms of the GNU General Public License as published by
 
8
/// the Free Software Foundation; either version 2 of the License, or
 
9
/// (at your option) any later version.
 
10
///
 
11
/// Rheolef is distributed in the hope that it will be useful,
 
12
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
/// GNU General Public License for more details.
 
15
///
 
16
/// You should have received a copy of the GNU General Public License
 
17
/// along with Rheolef; if not, write to the Free Software
 
18
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
19
///
 
20
/// =========================================================================
 
21
//
 
22
// 2D_s_D_s matrix on a surface, as defined by a level set
 
23
//
 
24
#include "2D_s_D_s.h"
 
25
#include "rheolef/ublas_matrix_range.h"
 
26
namespace rheolef {
 
27
extern void compute_matrix_2D_s_D_s (const geo_element& K, const ublas::vector<point>& x, const ublas::vector<Float>& f, ublas::matrix<Float>& mk);
 
28
} // namespace rheolef
 
29
using namespace rheolef;
 
30
using namespace std;
 
31
using namespace ublas;
 
32
 
 
33
void 
 
34
_2D_s_D_s::operator() (const geo_element& K, ublas::matrix<Float>& b) const
 
35
{
 
36
    size_type nloc = K.size();
 
37
    ublas::vector<point> x (nloc);
 
38
    ublas::vector<Float> f (nloc);
 
39
    tiny_vector<size_t> bgd_idof (nloc);
 
40
    size_type d = coordinate_dimension();
 
41
    b.resize (d*nloc, d*nloc);
 
42
    const space& Bh = get_first_space();
 
43
    Bh.set_global_dof (K, bgd_idof);
 
44
    const field& phi_h = _wh;
 
45
    const geo& bgd_lambda = phi_h.get_geo();
 
46
    for (size_t iloc = 0; iloc < nloc; iloc++) {
 
47
      x[iloc] = bgd_lambda.vertex(K[iloc]);
 
48
      f[iloc] = phi_h.at (bgd_idof[iloc]);
 
49
    }
 
50
    compute_matrix_2D_s_D_s (K, x, f, b);
 
51
}
 
52
void
 
53
_2D_s_D_s::check_after_initialize () const
 
54
{
 
55
  size_type d = coordinate_dimension();
 
56
  check_macro (
 
57
        d == get_second_space().n_component(),
 
58
        "unsupported non-vectorial second space for `2D_s_D_s' form");
 
59
  check_macro (
 
60
        d == get_first_space().n_component(),
 
61
        "unsupported non-vectorial first space for `2D_s_D_s' form");
 
62
  check_macro (get_first_space().get_approx()  == "P1", "`2D_s_D_s' form: P1 approx required");
 
63
  check_macro (get_second_space().get_approx() == "P1", "`2D_s_D_s' form: P1 approx required");
 
64
  check_macro (is_weighted(), "`2D_s_D_s' form may specify the level set function");
 
65
}