Changeset 3466

Timestamp:

04/08/10 10:27:37 (15 years ago)

Author:

seroussi

Message:

minor: added gauss points for Segments

File:

• issm/trunk/src/c/shared/Numerics/GaussPoints.cpp (modified) (1 diff)

issm/trunk/src/c/shared/Numerics/GaussPoints.cpp

@@ -2041 +2041 @@
                 
         }
+        else if (num_gauss==4){
+                //order=3, num_gauss=4. Can integrate polynomials of degree 0 to 5
+
+                gauss_weights=(double*)xmalloc(num_gauss*sizeof(double));
+                gauss_weights[0]=0.347854845137454;
+                gauss_weights[1]=0.652145154862546;
+                gauss_weights[2]=0.652145154862546;
+                gauss_weights[3]=0.347854845137454;
+                
+                segment_coord=(double*)xmalloc(num_gauss*sizeof(double));
+                segment_coord[0]=-0.861136311594053;
+                segment_coord[1]=-0.339981043584856;
+                segment_coord[2]=0.339981043584856;
+                segment_coord[3]=0.861136311594053;
+                
+        }
+        else if (num_gauss==5){
+                //order=3, num_gauss=4. Can integrate polynomials of degree 0 to 5
+
+                gauss_weights=(double*)xmalloc(num_gauss*sizeof(double));
+                gauss_weights[0]=0.236926885056189;
+                gauss_weights[1]=0.478628670499366;
+                gauss_weights[2]=0.568888888888889;
+                gauss_weights[3]=0.478628670499366;
+                gauss_weights[4]=0.236926885056189;
+                
+                segment_coord=(double*)xmalloc(num_gauss*sizeof(double));
+                segment_coord[0]=-0.906179845938664;
+                segment_coord[1]=-0.538469310105683;
+                segment_coord[2]=0.0;
+                segment_coord[3]=0.538469310105683;
+                segment_coord[4]=0.906179845938664;
+                
+        }
         else{
                 _printf_("GaussSegment error message: order not supported yet");