  a(1) = 4.0/5.0
  a(2) = -1.0/5.0
  a(3) = 4.0/105.0
  a(4) = -1.0/280.0

  q(1,1) = -1.695543604431898508749855654248370812054
  q(2,1) = 2.244525109969933844775276252674718713114
  q(3,1) = -0.002163206607021876385102025416828692597556
  q(4,1) = -0.8963677947307838167920190976769921451579
  q(5,1) = 0.2272938673777916084251861567415263328923
  q(6,1) = 0.1564962758424078514454013191006671546757
  q(7,1) = 0.003067629343740613524557473056482141287158
  q(8,1) = -0.03730827676416971624344442423120269215938
  q(9,1) = 0
  q(10,1) = 0
  q(11,1) = 0
  q(12,1) = 0

  q(1,2) = -0.4338209217401506177055540526837828066711
  q(2,2) = 0
  q(3,2) = 0.06631961192627293470305188523936353462439
  q(4,2) = 0.6021342012739069461791324177355805603656
  q(5,2) = -0.2213339967864188453859900127368926875520
  q(6,2) = -0.001759980003751516130708214237001705688807
  q(7,2) = -0.03586646217917856587352426166797879385310
  q(8,2) = 0.02432754750931966421359223835071189877503
  q(9,2) = 0
  q(10,2) = 0
  q(11,2) = 0
  q(12,2) = 0

  q(1,3) = 0.002477771724790106958560398469919805908522
  q(2,3) = -0.3930241559935857537756812928554282732813
  q(3,3) = 0
  q(4,3) = -0.6080263371305341035810612202824357742461
  q(5,3) = 2.372533124692864412954923458457782022478
  q(6,3) = -2.126112217612194556872952028715623246817
  q(7,3) = 0.9075309435559160631495435310307745215131
  q(8,3) = -0.1553791292372561688333328461049890555549
  q(9,3) = 0
  q(10,3) = 0
  q(11,3) = 0
  q(12,3) = 0

  q(1,4) = 0.1470043328582935680963279007411074464266
  q(2,4) = -0.5109179516689331400658365102850902594996
  q(3,4) = 0.08705685833207777685654755900085198684605
  q(4,4) = 0
  q(5,4) = -0.1306945333157231839535552435788086367322
  q(6,4) = 0.5981933016362239542598518794410035592720
  q(7,4) = -0.2031099149801233772210520666270517105660
  q(8,4) = 0.01246790713818440202771648130798761425319
  q(9,4) = 0
  q(10,4) = 0
  q(11,4) = 0
  q(12,4) = 0

  q(1,5) = -0.1624073990846984662267508265053107632238
  q(2,5) = 0.8182390368133059538132603839842499379861
  q(3,5) = -1.480018301574606037840376638130713134769
  q(4,5) = 0.5694185675497882973361521309100738568707
  q(5,5) = 0
  q(6,5) = 0.02929508293831017534661096502884570435894
  q(7,5) = 0.2825910005984090924748504462528740123097
  q(8,5) = -0.04846434332860791750685068920563907354793
  q(9,5) = -0.008653643911901097396895772334380539984709
  q(10,5) = 0
  q(11,5) = 0
  q(12,5) = 0

  q(1,6) = -0.03609686950604370828405582087760384548297
  q(2,6) = 0.002100328577309696555612805397001149641337
  q(3,6) = 0.4281425817419204360479874008765051769418
  q(4,6) = -0.8413238238875046736839550701643761633826
  q(5,6) = -0.009456755727630000971085851751473603919416
  q(6,6) = 0
  q(7,6) = 0.5172389789041733771234394082779034212042
  q(8,6) = -0.08760813565107717873901807278223424202774
  q(9,6) = 0.02979718129528502284256573906127239395909
  q(10,6) = -0.002793485746432970891490538036994286933664
  q(11,6) = 0
  q(12,6) = 0

  q(1,7) = -0.0009797678134978722516935350416937185004514
  q(2,7) = 0.05926834509975463070197111976550139351983
  q(3,7) = -0.2530570463899371286637621744353665227114
  q(4,7) = 0.3955555826583941989223787901885502942688
  q(5,7) = -0.1263166266018192756531792392597193430222
  q(6,7) = -0.7162192643577544899896537844517283192029
  q(7,7) = 0
  q(8,7) = 0.8209721963136350137518635528607308559584
  q(9,7) = -0.2166153552278720352907291696202456084323
  q(10,7) = 0.04126006766245181624585317516576106827282
  q(11,7) = -0.003868131343354857773048735171790100150577
  q(12,7) = 0

  q(1,8) = 0.01090012273307913185972171872891927027272
  q(2,8) = -0.03677379943661633440187210478144113487019
  q(3,8) = 0.03963287609450569641558898819403594845401
  q(4,8) = -0.02221139656912423686782339404035534965598
  q(5,8) = 0.01981665906734246925389947980969432794418
  q(6,8) = 0.1109698768905366566973791789571759073485
  q(7,8) = -0.7509903604688148129224205834189298636300
  q(8,8) = 0
  q(9,8) = 0.7926019635554773751160111698352821779736
  q(10,8) = -0.1981504908888693437790027924588205444934
  q(11,8) = 0.03774295064549892262457196046834677037969
  q(12,8) = -0.003538401623015523996053621293907509723096

  # from https://bitbucket.org/cactuscode/cactusnumerical/src/master/SummationByParts/src/GetScalProdDiag.F90
  CCTK_REAL, dimension(8), parameter :: bmask_8 = (/ 1498139.0_wp/5080320.0_wp,&
                                                     1107307.0_wp/725760.0_wp, &
                                                     20761.0_wp/80640.0_wp, &
                                                     1304999.0_wp/725760.0_wp, &
                                                     299527.0_wp/725760.0_wp, &
                                                     103097.0_wp/80640.0_wp, &
                                                     670091.0_wp/725760.0_wp, &
                                                     5127739.0_wp/5080320.0_wp/)
