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

  q(1,1) = -21600.0/13649.0
  q(2,1) = 81763.0/40947.0
  q(3,1) = 131.0/27298.0
  q(4,1) = -9143.0/13649.0
  q(5,1) = 20539.0/81894.0
  q(6,1) = 0
  q(7,1) = 0
  q(8,1) = 0
  q(9,1) = 0

  q(1,2) = -81763.0/180195.0
  q(2,2) = 0
  q(3,2) = 7357.0/36039.0
  q(4,2) = 30637.0/72078.0
  q(5,2) = -2328.0/12013.0
  q(6,2) = 6611.0/360390.0
  q(7,2) = 0
  q(8,2) = 0
  q(9,2) = 0

  q(1,3) = -131.0/54220.0
  q(2,3) = -7357.0/16266.0
  q(3,3) = 0
  q(4,3) = 645.0/2711.0
  q(5,3) = 11237.0/32532.0
  q(6,3) = -3487.0/27110.0
  q(7,3) = 0
  q(8,3) = 0
  q(9,3) = 0

  q(1,4) = 9143.0/53590.0
  q(2,4) = -30637.0/64308.0
  q(3,4) =  -645.0/5359.0
  q(4,4) = 0
  q(5,4) = 13733.0/32154.0
  q(6,4) = -67.0/4660.0
  q(7,4) = 72.0/5359.0
  q(8,4) = 0
  q(9,4) = 0

  q(1,5) = -20539.0/236310.0
  q(2,5) = 2328.0/7877.0
  q(3,5) = -11237.0/47262.0
  q(4,5) = -13733.0/23631.0
  q(5,5) = 0
  q(6,5) = 89387.0/118155.0
  q(7,5) = -1296.0/7877.0
  q(8,5) = 144.0/7877.0
  q(9,5) = 0

  q(1,6) = 0
  q(2,6) = -6611.0/262806.0
  q(3,6) = 3487.0/43801.0
  q(4,6) = 1541.0/87602.0
  q(5,6) = -89387.0/131403.0
  q(6,6) = 0
  q(7,6) = 32400.0/43801.0
  q(8,6) = -6480.0/43801.0
  q(9,6) = 720.0/43801.0

  # 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 /)
