  a(1) = 3.0/4.0
  a(2) = -3.0/20.0
  a(3) = 1.0/60.0

  q(1,1) = -1.582533518939116418785258993332844897062
  q(2,1) = 2.033426786468126253898161347360808173712
  q(3,1) = -0.1417052898146741610733887894481170575600
  q(4,1) = -0.4501096599735708523162117824920488989702
  q(5,1) = 0.1042956382142412661862395105494407610836
  q(6,1) = 0.03662604404499391209045870736276191879693
  q(7,1) = 0
  q(8,1) = 0
  q(9,1) = 0

  q(1,2) = -0.4620701275035953590186631853846278325646
  q(2,2) = 0
  q(3,2) = 0.2873679417026202568532985205129449923126
  q(4,2) = 0.2585974499280928196267362923074433487080
  q(5,2) = -0.06894808744606961472005221923058251153103
  q(6,2) = -0.01494717668104810274131940820517799692506
  q(7,2) = 0
  q(8,2) = 0
  q(9,2) = 0

  q(1,3) = 0.07134398748360337973038301686379010397038
  q(2,3) = -0.6366933020423417826592908754928085932593
  q(3,3) = 0
  q(4,3) = 0.6067199374180168986519150843189505198519
  q(5,3) = -0.02338660408468356531858175098561718651857
  q(6,3) = -0.01798401877459493040442547470431484404443
  q(7,3) = 0
  q(8,3) = 0
  q(9,3) = 0

  q(1,4) = 0.1146397975178068401430112823144985150596
  q(2,4) = -0.2898424301162697370942324201800071793273
  q(3,4) = -0.3069262456316931913128086944558079603132
  q(4,4) = 0
  q(5,4) = 0.5203848121857539166740071338174418292578
  q(6,4) = -0.05169127637022742348368508279860701098408
  q(7,4) = 0.01343534241462959507370778130248180630715
  q(8,4) = 0
  q(9,4) = 0

  q(1,5) = -0.03614399304268576976452921364705641609825
  q(2,5) = 0.1051508663818248421520867474440761344449
  q(3,5) = 0.01609777419666805778308369351834662756172
  q(4,5) = -0.7080721616106272031118456849378369336023
  q(5,5) = 0
  q(6,5) = 0.7692160858661111736140494493705980473867
  q(7,5) = -0.1645296432652024882569506157166433921544
  q(8,5) = 0.01828107147391138758410562396851593246160
  q(9,5) = 0

  q(1,6) = -0.01141318406360863692889821914555232596651
  q(2,6) = 0.02049729840293952857599941220163960606616
  q(3,6) = 0.01113095018331244864875173213474522093204
  q(4,6) = 0.06324365883611076515355091406993789453750
  q(5,6) = -0.6916640154753724474963890679085181638850
  q(6,6) = 0
  q(7,6) = 0.7397091390607520376247117645715851236273
  q(8,6) = -0.1479418278121504075249423529143170247255
  q(9,6) = 0.01643798086801671194721581699047966941394

  # from https://bitbucket.org/cactuscode/cactusnumerical/src/master/SummationByParts/src/GetScalProdDiag.F90
  CCTK_REAL, dimension(6), parameter :: bmask_6 = (/ 13649.0_wp/43200._wp, &
                                                     12013.0_wp/8640._wp, &
                                                     2711.0_wp/4320.0_wp, &
                                                     5359.0_wp/4320.0_wp, &
                                                     7877.0_wp/8640.0_wp, &
                                                     43801.0_wp/43200.0_wp /)
