← Back to branch summary

~ubuntu-branches/ubuntu/raring/rheolef/raring-proposed

« back to all changes in this revision

Viewing changes to nfem/pform_element/s_grad_grad.cc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito, Pierre Saramito, Sylvestre Ledru
  • Date: 2012-05-14 14:02:09 UTC
  • mfrom: (1.1.6)
  • Revision ID: package-import@ubuntu.com-20120514140209-dzbdlidkotyflf9e
Tags: 6.1-1
http://bugs.debian.org/671996
[ Pierre Saramito ]
* New upstream release 6.1 (minor changes):
  - support arbitrarily polynomial order Pk
  - source code supports g++-4.7 (closes: #671996)

[ Sylvestre Ledru ]
* update of the watch file

Show diffs side-by-side

added added

removed removed

Lines of Context:
nfem/pform_element/s_grad_grad.cc
29
29
void 
30
30
s_grad_grad<T,M>::operator() (const geo_element& K, ublas::matrix<T>& a) const
31
31
{
 
32
  check_macro (! base::is_on_band(), "unsupported banded level set extension");
 
33
 
32
34
  // 3d: curl(curl(psi)) = - Laplacian(psi) + grad(div(psi))
33
35
  //  => not s.d.p matrix but with look at pis with div(psi)=0 : stream function
34
36
  if (K.dimension() == 3) {
35
37
    size_type n_comp = base::get_first_space().size();
36
 
    matrix<Float> ai;
 
38
    matrix<T> ai;
37
39
    base::build_scalar_grad_grad (K, ai);
38
40
    size_type n1 = ai.size1();
39
41
    size_type n2 = ai.size2();