"Fossies" - the Fresh Open Source Software Archive

Member "getdp-3.1.0-source/Kernel/Gauss_Quadrangle.h" (29 Dec 2018, 2206 Bytes) of package /linux/privat/getdp-3.1.0-source.tgz:


As a special service "Fossies" has tried to format the requested source page into HTML format using (guessed) C and C++ source code syntax highlighting (style: standard) with prefixed line numbers and code folding option. Alternatively you can here view or download the uninterpreted source code file. For more information about "Gauss_Quadrangle.h" see the Fossies "Dox" file reference documentation and the latest Fossies "Diffs" side-by-side code changes report: 3.0.4_vs_3.1.0.

    1 // GetDP - Copyright (C) 1997-2019 P. Dular and C. Geuzaine, University of Liege
    2 //
    3 // See the LICENSE.txt file for license information. Please report all
    4 // issues on https://gitlab.onelab.info/getdp/getdp/issues.
    5 
    6 /* 1 integration point */
    7 
    8 static double xq1[1] = {0.};
    9 static double yq1[1] = {0.};
   10 static double pq1[1] = {4.};
   11 
   12 /* 3 integration points */
   13 
   14 static double xq3[3] = {0.816496580928,-0.408248290464,-0.408248290464};
   15 static double yq3[3] = {0.0,0.840896415255,-0.840896415255};
   16 static double pq3[3] = {1.33333333333,1.33333333333,1.33333333333};
   17 
   18 /* 4 integration points */
   19 
   20 static double xq4[4] = {0.577350269189,-0.577350269189,0.577350269189,-0.577350269189};
   21 static double yq4[4] = {0.577350269189,0.577350269189,-0.577350269189,-0.577350269189};
   22 static double pq4[4] = {1.,1.,1.,1.};
   23 
   24 /* 7 integration points */
   25 
   26 // WARNING: It seems that there is a problem with the values with 7 integration points ! Geuzaine (16/06/2016)
   27 
   28 static double xq7[7] = {0.0,-0.683130051064,0.683130051064,0.890654421782,
   29             -0.890654421782,0.374256642286,-0.374256642286};
   30 static double yq7[7] = {0.0,-0.683130051064,0.683130051064,-0.374256642286,0.374256642286,
   31             -0.890654421782,0.890654421782};
   32 static double pq7[7] = {1.142857142857,0.595238095238,0.595238095238,
   33             0.416666666666,0.416666666666,0.416666666666,0.416666666666};
   34 
   35 /*
   36    GAUSS QUADRANGLE WITH 1/R SINGULARITY OVER NODE (-1,-1,0)
   37 
   38    ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia,
   39          Some numerical integration formulae over triangles and
   40      squares with a 1/R singularity,
   41      Appl. Math. Modelling, Vol 5, June 1981, pp 209--211
   42 */
   43 
   44 /* 1 integration point */
   45 
   46 static double xqs1[1] = {-0.26501817};
   47 static double yqs1[1] = {-0.26501817};
   48 static double pqs1[1] = {3.52549435};
   49 
   50 /* 3 integration points */
   51 
   52 static double xqs3[3] = {-0.58105530,1.0,-0.21877566};
   53 static double yqs3[3] = {-0.58105530,-0.21877566,1.0};
   54 static double pqs3[3] = {2.37881900,0.57333767,0.57333767};
   55 
   56 /* 4 integration points */
   57 
   58 static double xqs4[4] = {-0.37512304,0.69629093,-0.92928746,-0.15602536};
   59 static double yqs4[4] = {-0.92928746,-0.15602536,-0.37512304,0.69629093};
   60 static double pqs4[4] = {1.02276580,0.73998134,1.02276580,0.73998134};