← Back to branch summary

~ubuntu-branches/debian/squeeze/gmsh/squeeze

« back to all changes in this revision

Viewing changes to Numeric/GaussQuadratureLin.cpp

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme, Christophe Prud'homme
  • Date: 2009-07-13 15:49:21 UTC
  • mfrom: (7.1.4 sid)
  • Revision ID: james.westby@ubuntu.com-20090713154921-zer07j8wixwa07ig
Tags: 2.3.1.dfsg-4
http://bugs.debian.org/529972
[Christophe Prud'homme]
* Bug fix: "gmsh with cgns write support", thanks to Oliver Borm
  (Closes: #529972).
* debian/rules: make sure that Gmsh is built with occ support on all
  platforms thanks to Denis Barbier (#536435).

Show diffs side-by-side

added added

removed removed

Lines of Context:
Numeric/GaussQuadratureLin.cpp
 
1
// Gmsh - Copyright (C) 1997-2009 C. Geuzaine, J.-F. Remacle
 
2
//
 
3
// See the LICENSE.txt file for license information. Please report all
 
4
// bugs and problems to <gmsh@geuz.org>.
 
5
 
 
6
#include "Gauss.h"
 
7
#include "GaussLegendre1D.h"
 
8
 
 
9
IntPt * GQL[20] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
 
10
 
 
11
IntPt *getGQLPts(int order)
 
12
{ 
 
13
  int n = (order+1)/2;
 
14
  int index = n;
 
15
  if(!GQL[index]) {
 
16
    double *pt,*wt;
 
17
    gmshGaussLegendre1D(n,&pt,&wt);
 
18
    GQL[index] = new IntPt[n];
 
19
    for(int i=0; i < n; i++) {
 
20
      GQL[index][i].pt[0] = pt[i];
 
21
      GQL[index][i].pt[1] = 0.0;
 
22
      GQL[index][i].pt[2] = 0.0;
 
23
      GQL[index][i].weight = wt[i];
 
24
    }
 
25
  }
 
26
  return GQL[index];
 
27
}
 
28
 
 
29
int getNGQLPts(int order)
 
30
{ 
 
31
  return (order+1)/2;
 
32
}