← Back to branch summary

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

« back to all changes in this revision

Viewing changes to nfem/lib/form_element_rep.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/lib/form_element_rep.cc
452
452
  reference_element hat_K = K;
453
453
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
454
454
  m.clear();
455
 
  quadrature::const_iterator first = _quad.begin(hat_K);
456
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
455
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
456
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
457
457
  vector<Float> F_bx_table (get_first_basis().size(hat_K));
458
458
  vector<Float> F_by_table (get_second_basis().size(hat_K));
459
459
  tensor DF;
598
598
  reference_element hat_K = K;
599
599
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
600
600
  m.clear();
601
 
  quadrature::const_iterator first = _quad.begin(hat_K);
602
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
601
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
602
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
603
603
  tensor DF;
604
604
  tensor invDF;
605
605
  for (size_type q = 0; first != last; first++, q++) {
641
641
  reference_element hat_K = K;
642
642
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
643
643
  m.clear();
644
 
  quadrature::const_iterator first = _quad.begin(hat_K);
645
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
644
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
645
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
646
646
  vector<point> F_bx_table_grad (get_first_basis().size(hat_K));
647
647
  vector<Float> F_by_table (get_second_basis().size(hat_K));
648
648
  tensor DF;
674
674
  reference_element hat_K = K;
675
675
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
676
676
  m.clear();
677
 
  quadrature::const_iterator first = _quad.begin(hat_K);
678
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
677
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
678
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
679
679
  vector<tensor> F_bx_table (get_first_basis().size(hat_K)); 
680
680
  vector<Float> F_by_table (get_second_basis().size(hat_K));
681
681
  tensor DF;
711
711
  reference_element hat_K = K;
712
712
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
713
713
  m.clear();
714
 
  quadrature::const_iterator first = _quad.begin(hat_K);
715
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
714
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
715
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
716
716
  vector<tensor > F_bx_table (get_first_basis().size(hat_K)); 
717
717
  vector<point> F_by_table (get_second_basis().size(hat_K));
718
718
  tensor DF;
748
748
  reference_element hat_K = K;
749
749
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
750
750
  m.clear();
751
 
  quadrature::const_iterator first = _quad.begin(hat_K);
752
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
751
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
752
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
753
753
  vector<tensor > F_bx_table (get_first_basis().size(hat_K));
754
754
  vector<tensor > F_by_table (get_second_basis().size(hat_K)); 
755
755
  tensor DF;
786
786
  reference_element hat_K = K;
787
787
  m.resize (get_second_basis().size(hat_K), get_first_basis().size(hat_K));
788
788
  m.clear();
789
 
  quadrature::const_iterator first = _quad.begin(hat_K);
790
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
789
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
790
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
791
791
  vector<point> F_bx_table (get_first_basis().size(hat_K));
792
792
  vector<point> F_by_table (get_second_basis().size(hat_K));
793
793
  tensor DF;
847
847
  size_type ny = get_second_basis().size(hat_K);
848
848
  m.resize (map_d*ny, map_d*nx);
849
849
  m.clear();
850
 
  quadrature::const_iterator first = _quad.begin(hat_K);
851
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
850
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
851
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
852
852
  tensor DF;
853
853
  tensor invDF;
854
854
  tensor invDFt;
901
901
  size_type ny = get_second_basis().size(hat_K);
902
902
  m.resize (map_d*ny, map_d*nx);
903
903
  m.clear();
904
 
  quadrature::const_iterator first = _quad.begin(hat_K);
905
 
  quadrature::const_iterator last  = _quad.end  (hat_K);
 
904
  quadrature<Float>::const_iterator first = _quad.begin(hat_K);
 
905
  quadrature<Float>::const_iterator last  = _quad.end  (hat_K);
906
906
  tensor DF;
907
907
  tensor invDF;
908
908
  tensor invDFt;