C     1D Gauss-Legendre weights for N = 4
      wg(  1) = 0.3478548/2.0
      wg(  2) = 0.6521452/2.0
      wg(  3) = 0.6521452/2.0
      wg(  4) = 0.3478548/2.0

C     1D Gauss-Legendre weights for N = 8
      wg(  5) = 0.1012285/2.0
      wg(  6) = 0.2223810/2.0
      wg(  7) = 0.3137066/2.0
      wg(  8) = 0.3626838/2.0
      wg(  9) = 0.3626838/2.0
      wg( 10) = 0.3137066/2.0
      wg( 11) = 0.2223810/2.0
      wg( 12) = 0.1012285/2.0

C     1D Gauss-Legendre weights for N = 16
      wg( 13) = 0.0271524/2.0
      wg( 14) = 0.0622535/2.0
      wg( 15) = 0.0951585/2.0
      wg( 16) = 0.1246290/2.0
      wg( 17) = 0.1495960/2.0
      wg( 18) = 0.1691565/2.0
      wg( 19) = 0.1826034/2.0
      wg( 20) = 0.1894506/2.0
      wg( 21) = 0.1894506/2.0
      wg( 22) = 0.1826034/2.0
      wg( 23) = 0.1691565/2.0
      wg( 24) = 0.1495960/2.0
      wg( 25) = 0.1246290/2.0
      wg( 26) = 0.0951585/2.0
      wg( 27) = 0.0622535/2.0
      wg( 28) = 0.0271524/2.0

C     1D Gauss-Legendre weights for N = 32
      wg( 29) = 0.00701861/2.0
      wg( 30) = 0.01627439/2.0
      wg( 31) = 0.02539207/2.0
      wg( 32) = 0.03427386/2.0
      wg( 33) = 0.04283590/2.0
      wg( 34) = 0.05099806/2.0
      wg( 35) = 0.05868409/2.0
      wg( 36) = 0.06582222/2.0
      wg( 37) = 0.07234579/2.0
      wg( 38) = 0.07819390/2.0
      wg( 39) = 0.08331192/2.0
      wg( 40) = 0.08765209/2.0
      wg( 41) = 0.09117388/2.0
      wg( 42) = 0.09384440/2.0
      wg( 43) = 0.09563872/2.0
      wg( 44) = 0.09654009/2.0
      wg( 45) = 0.09654009/2.0
      wg( 46) = 0.09563872/2.0
      wg( 47) = 0.09384440/2.0
      wg( 48) = 0.09117388/2.0
      wg( 49) = 0.08765209/2.0
      wg( 50) = 0.08331192/2.0
      wg( 51) = 0.07819390/2.0
      wg( 52) = 0.07234579/2.0
      wg( 53) = 0.06582222/2.0
      wg( 54) = 0.05868409/2.0
      wg( 55) = 0.05099806/2.0
      wg( 56) = 0.04283590/2.0
      wg( 57) = 0.03427386/2.0
      wg( 58) = 0.02539207/2.0
      wg( 59) = 0.01627439/2.0
      wg( 60) = 0.00701861/2.0

C     1D Gauss-Legendre weights for N = 48
      wg( 61) = 0.00315335/2.0
      wg( 62) = 0.00732755/2.0
      wg( 63) = 0.01147723/2.0
      wg( 64) = 0.01557932/2.0
      wg( 65) = 0.01961616/2.0
      wg( 66) = 0.02357076/2.0
      wg( 67) = 0.02742651/2.0
      wg( 68) = 0.03116723/2.0
      wg( 69) = 0.03477722/2.0
      wg( 70) = 0.03824135/2.0
      wg( 71) = 0.04154508/2.0
      wg( 72) = 0.04467456/2.0
      wg( 73) = 0.04761666/2.0
      wg( 74) = 0.05035904/2.0
      wg( 75) = 0.05289019/2.0
      wg( 76) = 0.05519950/2.0
      wg( 77) = 0.05727729/2.0
      wg( 78) = 0.05911484/2.0
      wg( 79) = 0.06070444/2.0
      wg( 80) = 0.06203942/2.0
      wg( 81) = 0.06311419/2.0
      wg( 82) = 0.06392424/2.0
      wg( 83) = 0.06446616/2.0
      wg( 84) = 0.06473770/2.0
      wg( 85) = 0.06473770/2.0
      wg( 86) = 0.06446616/2.0
      wg( 87) = 0.06392424/2.0
      wg( 88) = 0.06311419/2.0
      wg( 89) = 0.06203942/2.0
      wg( 90) = 0.06070444/2.0
      wg( 91) = 0.05911484/2.0
      wg( 92) = 0.05727729/2.0
      wg( 93) = 0.05519950/2.0
      wg( 94) = 0.05289019/2.0
      wg( 95) = 0.05035904/2.0
      wg( 96) = 0.04761666/2.0
      wg( 97) = 0.04467456/2.0
      wg( 98) = 0.04154508/2.0
      wg( 99) = 0.03824135/2.0
      wg(100) = 0.03477722/2.0
      wg(101) = 0.03116723/2.0
      wg(102) = 0.02742651/2.0
      wg(103) = 0.02357076/2.0
      wg(104) = 0.01961616/2.0
      wg(105) = 0.01557932/2.0
      wg(106) = 0.01147723/2.0
      wg(107) = 0.00732755/2.0
      wg(108) = 0.00315335/2.0

C     1D Gauss-Legendre abscissas for N = 4
      xg(  1) = -0.8611363/2.0 + 0.5
      xg(  2) = -0.3399810/2.0 + 0.5
      xg(  3) =  0.3399810/2.0 + 0.5
      xg(  4) =  0.8611363/2.0 + 0.5

C     1D Gauss-Legendre abscissas for N = 8
      xg(  5) = -0.9602899/2.0 + 0.5
      xg(  6) = -0.7966665/2.0 + 0.5
      xg(  7) = -0.5255324/2.0 + 0.5
      xg(  8) = -0.1834346/2.0 + 0.5
      xg(  9) =  0.1834346/2.0 + 0.5
      xg( 10) =  0.5255324/2.0 + 0.5
      xg( 11) =  0.7966665/2.0 + 0.5
      xg( 12) =  0.9602899/2.0 + 0.5

C     1D Gauss-Legendre abscissas for N = 16
      xg( 13) = -0.9894009/2.0 + 0.5
      xg( 14) = -0.9445750/2.0 + 0.5
      xg( 15) = -0.8656312/2.0 + 0.5
      xg( 16) = -0.7554044/2.0 + 0.5
      xg( 17) = -0.6178762/2.0 + 0.5
      xg( 18) = -0.4580168/2.0 + 0.5
      xg( 19) = -0.2816036/2.0 + 0.5
      xg( 20) = -0.0950125/2.0 + 0.5
      xg( 21) =  0.0950125/2.0 + 0.5
      xg( 22) =  0.2816036/2.0 + 0.5
      xg( 23) =  0.4580168/2.0 + 0.5
      xg( 24) =  0.6178762/2.0 + 0.5
      xg( 25) =  0.7554044/2.0 + 0.5
      xg( 26) =  0.8656312/2.0 + 0.5
      xg( 27) =  0.9445750/2.0 + 0.5
      xg( 28) =  0.9894009/2.0 + 0.5

C     1D Gauss-Legendre abscissas for N = 32
      xg( 29) = -0.9972639/2.0 + 0.5
      xg( 30) = -0.9856115/2.0 + 0.5
      xg( 31) = -0.9647623/2.0 + 0.5
      xg( 32) = -0.9349061/2.0 + 0.5
      xg( 33) = -0.8963212/2.0 + 0.5
      xg( 34) = -0.8493676/2.0 + 0.5
      xg( 35) = -0.7944838/2.0 + 0.5
      xg( 36) = -0.7321821/2.0 + 0.5
      xg( 37) = -0.6630443/2.0 + 0.5
      xg( 38) = -0.5877158/2.0 + 0.5
      xg( 39) = -0.5068999/2.0 + 0.5
      xg( 40) = -0.4213513/2.0 + 0.5
      xg( 41) = -0.3318686/2.0 + 0.5
      xg( 42) = -0.2392874/2.0 + 0.5
      xg( 43) = -0.1444720/2.0 + 0.5
      xg( 44) = -0.04830767/2.0 + 0.5
      xg( 45) =  0.04830767/2.0 + 0.5
      xg( 46) =  0.1444720/2.0 + 0.5
      xg( 47) =  0.2392874/2.0 + 0.5
      xg( 48) =  0.3318686/2.0 + 0.5
      xg( 49) =  0.4213513/2.0 + 0.5
      xg( 50) =  0.5068999/2.0 + 0.5
      xg( 51) =  0.5877158/2.0 + 0.5
      xg( 52) =  0.6630443/2.0 + 0.5
      xg( 53) =  0.7321821/2.0 + 0.5
      xg( 54) =  0.7944838/2.0 + 0.5
      xg( 55) =  0.8493676/2.0 + 0.5
      xg( 56) =  0.8963212/2.0 + 0.5
      xg( 57) =  0.9349061/2.0 + 0.5
      xg( 58) =  0.9647623/2.0 + 0.5
      xg( 59) =  0.9856115/2.0 + 0.5
      xg( 60) =  0.9972639/2.0 + 0.5

C     1D Gauss-Legendre abscissas for N = 48
      xg( 61) = -0.9987710/2.0 + 0.5
      xg( 62) = -0.9935302/2.0 + 0.5
      xg( 63) = -0.9841246/2.0 + 0.5
      xg( 64) = -0.9705916/2.0 + 0.5
      xg( 65) = -0.9529877/2.0 + 0.5
      xg( 66) = -0.9313867/2.0 + 0.5
      xg( 67) = -0.9058791/2.0 + 0.5
      xg( 68) = -0.8765720/2.0 + 0.5
      xg( 69) = -0.8435883/2.0 + 0.5
      xg( 70) = -0.8070662/2.0 + 0.5
      xg( 71) = -0.7671590/2.0 + 0.5
      xg( 72) = -0.7240341/2.0 + 0.5
      xg( 73) = -0.6778724/2.0 + 0.5
      xg( 74) = -0.6288674/2.0 + 0.5
      xg( 75) = -0.5772247/2.0 + 0.5
      xg( 76) = -0.5231610/2.0 + 0.5
      xg( 77) = -0.4669029/2.0 + 0.5
      xg( 78) = -0.4086865/2.0 + 0.5
      xg( 79) = -0.3487559/2.0 + 0.5
      xg( 80) = -0.2873625/2.0 + 0.5
      xg( 81) = -0.2247638/2.0 + 0.5
      xg( 82) = -0.1612224/2.0 + 0.5
      xg( 83) = -0.09700470/2.0 + 0.5
      xg( 84) = -0.03238017/2.0 + 0.5
      xg( 85) =  0.03238017/2.0 + 0.5
      xg( 86) =  0.09700470/2.0 + 0.5
      xg( 87) =  0.1612224/2.0 + 0.5
      xg( 88) =  0.2247638/2.0 + 0.5
      xg( 89) =  0.2873625/2.0 + 0.5
      xg( 90) =  0.3487559/2.0 + 0.5
      xg( 91) =  0.4086865/2.0 + 0.5
      xg( 92) =  0.4669029/2.0 + 0.5
      xg( 93) =  0.5231610/2.0 + 0.5
      xg( 94) =  0.5772247/2.0 + 0.5
      xg( 95) =  0.6288674/2.0 + 0.5
      xg( 96) =  0.6778724/2.0 + 0.5
      xg( 97) =  0.7240341/2.0 + 0.5
      xg( 98) =  0.7671590/2.0 + 0.5
      xg( 99) =  0.8070662/2.0 + 0.5
      xg(100) =  0.8435883/2.0 + 0.5
      xg(101) =  0.8765720/2.0 + 0.5
      xg(102) =  0.9058791/2.0 + 0.5
      xg(103) =  0.9313867/2.0 + 0.5
      xg(104) =  0.9529877/2.0 + 0.5
      xg(105) =  0.9705916/2.0 + 0.5
      xg(106) =  0.9841246/2.0 + 0.5
      xg(107) =  0.9935302/2.0 + 0.5
      xg(108) =  0.9987710/2.0 + 0.5