From d7510e9159cb7b7fd7daa5b239e2ee6513eb3d7a Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <nicolas.marsic@gmail.com>
Date: Tue, 13 Nov 2012 12:51:23 +0000
Subject: [PATCH] Line Quadrature: bug fix -- Round Up
---
Numeric/GaussQuadratureLin.cpp | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/Numeric/GaussQuadratureLin.cpp b/Numeric/GaussQuadratureLin.cpp
index 931294f5f0..e340e57c6e 100644
--- a/Numeric/GaussQuadratureLin.cpp
+++ b/Numeric/GaussQuadratureLin.cpp
@@ -9,8 +9,11 @@
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)
-{
- int n = (order+1)/2;
+{
+ // Number of Gauss Point:
+ // (order + 1) / 2 *ROUNDED UP*
+
+ int n = (order + 1) / (double)2 + 0.5;
int index = n;
if(!GQL[index]) {
double *pt,*wt;
@@ -28,5 +31,5 @@ IntPt *getGQLPts(int order)
int getNGQLPts(int order)
{
- return (order+1)/2;
+ return (order + 1) / (double)2 + 0.5;
}
--
GitLab