1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
// Gmsh - Copyright (C) 1997-2020 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
#include "GmshMessage.h"
#include "GaussIntegration.h"
#include "GaussLegendre1D.h"
IntPt *GQL[20] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
IntPt *getGQLPts(int order)
{
// Number of Gauss Point:
// (order + 1) / 2 *ROUNDED UP*
int n = (order + 1) / (double)2 + 0.5;
int index = n;
if(index >= (int)(sizeof(GQL) / sizeof(IntPt *))) {
Msg::Error("Increase size of GQL in gauss quadrature line");
index = 0;
}
if(!GQL[index]) {
double *pt, *wt;
gmshGaussLegendre1D(n, &pt, &wt);
GQL[index] = new IntPt[n];
for(int i = 0; i < n; i++) {
GQL[index][i].pt[0] = pt[i];
GQL[index][i].pt[1] = 0.0;
GQL[index][i].pt[2] = 0.0;
GQL[index][i].weight = wt[i];
}
}
return GQL[index];
}
int getNGQLPts(int order) { return (order + 1) / (double)2 + 0.5; }
|
