!=========================================================================== ! subroutine der0_2D_cf_vc_pre( ierror, error, & xbegin,xend,ybegin,yend,zbegin,zend, & bnd4AA, bnd4AB, bnd4I, & bnd6AA, bnd6AB, bnd6AC, bnd6BB, bnd6BC, bnd6I, & bnd8AA, bnd8AB, bnd8AC, bnd8AD, bnd8BB, bnd8BC, & bnd8BD, bnd8CC, bnd8CD, bnd8I, bnd10I, & iperx,ipery,iperz, orderi, orderb, & biL, biR, bjL, bjR, bkL, bkR, & fiL, fiR, fjL, fjR, fkL, fkR, & iiL, iiR, ijL, ijR, ikL, ikR, & piL, piR, pjL, pjR, pkL, pkR, & dfiL,dfiR,dfjL,dfjR,dfkL,dfkR, & nproLx,nproRx,nproLy,nproRy,nproLz,nproRz, & proL,proR,promax,width) ! include 'der0_2D_cf_vc_include.h' ! !====================================================================== ! ! Initialize IERROR to zero implying no errors ! ierror = 0 ! ! Boundaries of usable f patch piL = iiL - biL piR = iiR + biR pjL = ijL - bjL pjR = ijR + bjR pkL = ikL - bkL pkR = ikR + bkR ! ! Is the usable patch width larger than the width of f? if( min0(piL-fiL,fiR-piR,pjL-fjL,fjR-pjR,pkL-fkL,fkR-pkR) & .lt. 0) then ierror = -10 error = 'DER_0_2D: Interior + Green width > function width' return endif ! ! Is the interior width larger than the width of df? if( min0(iiL-dfiL,dfiR-iiR,ijL-dfjL,dfjR-ijR,ikL-dfkL,dfkR-ikR) & .lt. 0) then ierror = -20 error = 'DER_0_2D: Interior width > derivative width' return endif ! ! Have the number of green cells gone negative anywhere? if( min0(biL,biR,bjL,bjR,bkL,bkR) .lt. 0 ) then ierror = -40 error = 'DER_0_2D: Green cell width < zero' return endif ! ! Has the interior width gone negative? if( min0(iiR-iiL,ijR-ijL,ikR-ikL) .lt. 0 ) then ierror = -50 error = 'DER_0_2D: Interior Width < zero' return endif ! ! Has the boundary order been chosen properly?? ! if( orderi .eq. 4) then if(orderb.ne.4) then ierror = -60 error = 'DER_0_2D: No such boundary order: orderb ' return endif elseif( orderi .eq. 6) then if(orderb.ne.6) then ierror = -61 error = 'DER_0_2D: No such boundary order: orderb ' return endif elseif( orderi .eq. 8) then if(orderb.ne.8) then ierror = -62 error = 'DER_0_2D: No such boundary order: orderb ' return endif endif ! ! Stencil Widths ! ! Centered stencils ! proL = -1 proR = -1 ! if( orderi .eq. 2 ) then proL = 0 proR = 0 elseif( orderi .eq. 4 ) then proL = 1 proR = 1 elseif( orderi .eq. 6 ) then proL = 2 proR = 2 elseif( orderi .eq. 8 ) then proL = 3 proR = 3 elseif( orderi .eq. 10 ) then proL = 4 proR = 4 endif width = proL + proR + 2 ! if( ( proL .lt. 0) .or. ( proR .lt. 0) ) then ierror = -80 error = 'DER_0_2D: No such scheme' return endif ! if( (iperx .lt. 0) .or. ( iperx .gt. 1) .or. & (ipery .lt. 0) .or. ( ipery .gt. 1) .or. & (iperz .lt. 0) .or. ( iperz .gt. 1) ) then ierror = -90 error = 'DER_0_2D: Periodicity flags are set incorrectly' return endif ! ! Find maximum protrusion. ! promax = max0( proL, proR ) ! ! Determine where to begin Interior DO loops and how many ! boundary points need to be closed. Periodic domains ! will always (better) have enough green cells. ! ! X-direction if(( biL .ge. promax ) .or. (iperx .eq. 1)) then nproLx = 0 xbegin = iiL else nproLx = promax - biL xbegin = iiL + nproLx endif ! if(( biR .ge. promax ) .or. (iperx .eq. 1)) then nproRx = 0 xend = iiR - 1 else nproRx = promax - biR xend = iiR - nproRx - 1 endif ! ! Y-direction ! if(( bjL .ge. promax ) .or. (ipery .eq. 1)) then nproLy = 0 ybegin = ijL else nproLy = promax - bjL ybegin = ijL + nproLy endif ! if(( bjR .ge. promax ) .or. (ipery .eq. 1)) then nproRy = 0 yend = ijR - 1 else nproRy = promax - bjR yend = ijR - nproRy - 1 endif ! ! Z-direction ! if(( bkL .ge. promax ) .or. (iperz .eq. 1)) then nproLz = 0 zbegin = ikL else nproLz = promax - bkL zbegin = ikL + nproLz endif ! if(( bkR .ge. promax ) .or. (iperz .eq. 1)) then nproRz = 0 zend = ikR - 1 else nproRz = promax - bkR zend = ikR - nproRz - 1 endif ! ! if( ( orderi .eq. 10) .and. ! & (min0(biL,biR,bjL,bjR,bkL,bkR) .lt. 4) ) then ! ierror = -100 ! error = 'DER_0_2D: 10th-order has no boundary closures' ! return ! endif ! ! Set up Boundary stencil coefficients ! ! 2nd-order ! bnd2AA(1,1) = 1.d0/4.d0 bnd2AA(1,2) = 1.d0/4.d0 bnd2AA(2,1) = 1.d0/4.d0 bnd2AA(2,2) = 1.d0/4.d0 ! ! 4th-order ! bnd4AA(1,1) = 25.d0/256.d0 bnd4AA(1,2) = 75.d0/256.d0 bnd4AA(1,3) = -25.d0/256.d0 bnd4AA(1,4) = 5.d0/256.d0 bnd4AA(2,1) = 75.d0/256.d0 bnd4AA(2,2) = 225.d0/256.d0 bnd4AA(2,3) = -75.d0/256.d0 bnd4AA(2,4) = 15.d0/256.d0 bnd4AA(3,1) = -25.d0/256.d0 bnd4AA(3,2) = -75.d0/256.d0 bnd4AA(3,3) = 25.d0/256.d0 bnd4AA(3,4) = -5.d0/256.d0 bnd4AA(4,1) = 5.d0/256.d0 bnd4AA(4,2) = 15.d0/256.d0 bnd4AA(4,3) = -5.d0/256.d0 bnd4AA(4,4) = 1.d0/256.d0 ! bnd4AB(1,1) = -5.d0/256.d0 bnd4AB(1,2) = 45.d0/256.d0 bnd4AB(1,3) = 45.d0/256.d0 bnd4AB(1,4) = -5.d0/256.d0 bnd4AB(2,1) = -15.d0/256.d0 bnd4AB(2,2) = 135.d0/256.d0 bnd4AB(2,3) = 135.d0/256.d0 bnd4AB(2,4) = -15.d0/256.d0 bnd4AB(3,1) = 5.d0/256.d0 bnd4AB(3,2) = -45.d0/256.d0 bnd4AB(3,3) = -45.d0/256.d0 bnd4AB(3,4) = 5.d0/256.d0 bnd4AB(4,1) = -1.d0/256.d0 bnd4AB(4,2) = 9.d0/256.d0 bnd4AB(4,3) = 9.d0/256.d0 bnd4AB(4,4) = -1.d0/256.d0 ! bnd4I(1,1) = 1.d0/256.d0 bnd4I(1,2) = -9.d0/256.d0 bnd4I(1,3) = -9.d0/256.d0 bnd4I(1,4) = 1.d0/256.d0 bnd4I(2,1) = -9.d0/256.d0 bnd4I(2,2) = 81.d0/256.d0 bnd4I(2,3) = 81.d0/256.d0 bnd4I(2,4) = -9.d0/256.d0 bnd4I(3,1) = -9.d0/256.d0 bnd4I(3,2) = 81.d0/256.d0 bnd4I(3,3) = 81.d0/256.d0 bnd4I(3,4) = -9.d0/256.d0 bnd4I(4,1) = 1.d0/256.d0 bnd4I(4,2) = -9.d0/256.d0 bnd4I(4,3) = -9.d0/256.d0 bnd4I(4,4) = 1.d0/256.d0 ! ! 6th-order ! bnd6AA(1,1) = 3969.d0/65536.d0 bnd6AA(1,2) = 19845.d0/65536.d0 bnd6AA(1,3) = -13230.d0/65536.d0 bnd6AA(1,4) = 7938.d0/65536.d0 bnd6AA(1,5) = -2835.d0/65536.d0 bnd6AA(1,6) = 441.d0/65536.d0 bnd6AA(2,1) = 19845.d0/65536.d0 bnd6AA(2,2) = 99225.d0/65536.d0 bnd6AA(2,3) = -66150.d0/65536.d0 bnd6AA(2,4) = 39690.d0/65536.d0 bnd6AA(2,5) = -14175.d0/65536.d0 bnd6AA(2,6) = 2205.d0/65536.d0 bnd6AA(3,1) = -13230.d0/65536.d0 bnd6AA(3,2) = -66150.d0/65536.d0 bnd6AA(3,3) = 44100.d0/65536.d0 bnd6AA(3,4) = -26460.d0/65536.d0 bnd6AA(3,5) = 9450.d0/65536.d0 bnd6AA(3,6) = -1470.d0/65536.d0 bnd6AA(4,1) = 7938.d0/65536.d0 bnd6AA(4,2) = 39690.d0/65536.d0 bnd6AA(4,3) = -26460.d0/65536.d0 bnd6AA(4,4) = 15876.d0/65536.d0 bnd6AA(4,5) = -5670.d0/65536.d0 bnd6AA(4,6) = 882.d0/65536.d0 bnd6AA(5,1) = -2835.d0/65536.d0 bnd6AA(5,2) = -14175.d0/65536.d0 bnd6AA(5,3) = 9450.d0/65536.d0 bnd6AA(5,4) = -5670.d0/65536.d0 bnd6AA(5,5) = 2025.d0/65536.d0 bnd6AA(5,6) = -315.d0/65536.d0 bnd6AA(6,1) = 441.d0/65536.d0 bnd6AA(6,2) = 2205.d0/65536.d0 bnd6AA(6,3) = -1470.d0/65536.d0 bnd6AA(6,4) = 882.d0/65536.d0 bnd6AA(6,5) = -315.d0/65536.d0 bnd6AA(6,6) = 49.d0/65536.d0 ! bnd6AB(1,1) = -441.d0/65536.d0 bnd6AB(1,2) = 6615.d0/65536.d0 bnd6AB(1,3) = 13230.d0/65536.d0 bnd6AB(1,4) = -4410.d0/65536.d0 bnd6AB(1,5) = 1323.d0/65536.d0 bnd6AB(1,6) = -189.d0/65536.d0 bnd6AB(2,1) = -2205.d0/65536.d0 bnd6AB(2,2) = 33075.d0/65536.d0 bnd6AB(2,3) = 66150.d0/65536.d0 bnd6AB(2,4) = -22050.d0/65536.d0 bnd6AB(2,5) = 6615.d0/65536.d0 bnd6AB(2,6) = -945.d0/65536.d0 bnd6AB(3,1) = 1470.d0/65536.d0 bnd6AB(3,2) = -22050.d0/65536.d0 bnd6AB(3,3) = -44100.d0/65536.d0 bnd6AB(3,4) = 14700.d0/65536.d0 bnd6AB(3,5) = -4410.d0/65536.d0 bnd6AB(3,6) = 630.d0/65536.d0 bnd6AB(4,1) = -882.d0/65536.d0 bnd6AB(4,2) = 13230.d0/65536.d0 bnd6AB(4,3) = 26460.d0/65536.d0 bnd6AB(4,4) = -8820.d0/65536.d0 bnd6AB(4,5) = 2646.d0/65536.d0 bnd6AB(4,6) = -378.d0/65536.d0 bnd6AB(5,1) = 315.d0/65536.d0 bnd6AB(5,2) = -4725.d0/65536.d0 bnd6AB(5,3) = -9450.d0/65536.d0 bnd6AB(5,4) = 3150.d0/65536.d0 bnd6AB(5,5) = -945.d0/65536.d0 bnd6AB(5,6) = 135.d0/65536.d0 bnd6AB(6,1) = -49.d0/65536.d0 bnd6AB(6,2) = 735.d0/65536.d0 bnd6AB(6,3) = 1470.d0/65536.d0 bnd6AB(6,4) = -490.d0/65536.d0 bnd6AB(6,5) = 147.d0/65536.d0 bnd6AB(6,6) = -21.d0/65536.d0 ! bnd6AC(1,1) = 189.d0/65536.d0 bnd6AC(1,2) = -1575.d0/65536.d0 bnd6AC(1,3) = 9450.d0/65536.d0 bnd6AC(1,4) = 9450.d0/65536.d0 bnd6AC(1,5) = -1575.d0/65536.d0 bnd6AC(1,6) = 189.d0/65536.d0 bnd6AC(2,1) = 945.d0/65536.d0 bnd6AC(2,2) = -7875.d0/65536.d0 bnd6AC(2,3) = 47250.d0/65536.d0 bnd6AC(2,4) = 47250.d0/65536.d0 bnd6AC(2,5) = -7875.d0/65536.d0 bnd6AC(2,6) = 945.d0/65536.d0 bnd6AC(3,1) = -630.d0/65536.d0 bnd6AC(3,2) = 5250.d0/65536.d0 bnd6AC(3,3) = -31500.d0/65536.d0 bnd6AC(3,4) = -31500.d0/65536.d0 bnd6AC(3,5) = 5250.d0/65536.d0 bnd6AC(3,6) = -630.d0/65536.d0 bnd6AC(4,1) = 378.d0/65536.d0 bnd6AC(4,2) = -3150.d0/65536.d0 bnd6AC(4,3) = 18900.d0/65536.d0 bnd6AC(4,4) = 18900.d0/65536.d0 bnd6AC(4,5) = -3150.d0/65536.d0 bnd6AC(4,6) = 378.d0/65536.d0 bnd6AC(5,1) = -135.d0/65536.d0 bnd6AC(5,2) = 1125.d0/65536.d0 bnd6AC(5,3) = -6750.d0/65536.d0 bnd6AC(5,4) = -6750.d0/65536.d0 bnd6AC(5,5) = 1125.d0/65536.d0 bnd6AC(5,6) = -135.d0/65536.d0 bnd6AC(6,1) = 21.d0/65536.d0 bnd6AC(6,2) = -175.d0/65536.d0 bnd6AC(6,3) = 1050.d0/65536.d0 bnd6AC(6,4) = 1050.d0/65536.d0 bnd6AC(6,5) = -175.d0/65536.d0 bnd6AC(6,6) = 21.d0/65536.d0 ! bnd6BB(1,1) = 49.d0/65536.d0 bnd6BB(1,2) = -735.d0/65536.d0 bnd6BB(1,3) = -1470.d0/65536.d0 bnd6BB(1,4) = 490.d0/65536.d0 bnd6BB(1,5) = -147.d0/65536.d0 bnd6BB(1,6) = 21.d0/65536.d0 bnd6BB(2,1) = -735.d0/65536.d0 bnd6BB(2,2) = 11025.d0/65536.d0 bnd6BB(2,3) = 22050.d0/65536.d0 bnd6BB(2,4) = -7350.d0/65536.d0 bnd6BB(2,5) = 2205.d0/65536.d0 bnd6BB(2,6) = -315.d0/65536.d0 bnd6BB(3,1) = -1470.d0/65536.d0 bnd6BB(3,2) = 22050.d0/65536.d0 bnd6BB(3,3) = 44100.d0/65536.d0 bnd6BB(3,4) = -14700.d0/65536.d0 bnd6BB(3,5) = 4410.d0/65536.d0 bnd6BB(3,6) = -630.d0/65536.d0 bnd6BB(4,1) = 490.d0/65536.d0 bnd6BB(4,2) = -7350.d0/65536.d0 bnd6BB(4,3) = -14700.d0/65536.d0 bnd6BB(4,4) = 4900.d0/65536.d0 bnd6BB(4,5) = -1470.d0/65536.d0 bnd6BB(4,6) = 210.d0/65536.d0 bnd6BB(5,1) = -147.d0/65536.d0 bnd6BB(5,2) = 2205.d0/65536.d0 bnd6BB(5,3) = 4410.d0/65536.d0 bnd6BB(5,4) = -1470.d0/65536.d0 bnd6BB(5,5) = 441.d0/65536.d0 bnd6BB(5,6) = -63.d0/65536.d0 bnd6BB(6,1) = 21.d0/65536.d0 bnd6BB(6,2) = -315.d0/65536.d0 bnd6BB(6,3) = -630.d0/65536.d0 bnd6BB(6,4) = 210.d0/65536.d0 bnd6BB(6,5) = -63.d0/65536.d0 bnd6BB(6,6) = 9.d0/65536.d0 ! bnd6BC(1,1) = -21.d0/65536.d0 bnd6BC(1,2) = 175.d0/65536.d0 bnd6BC(1,3) = -1050.d0/65536.d0 bnd6BC(1,4) = -1050.d0/65536.d0 bnd6BC(1,5) = 175.d0/65536.d0 bnd6BC(1,6) = -21.d0/65536.d0 bnd6BC(2,1) = 315.d0/65536.d0 bnd6BC(2,2) = -2625.d0/65536.d0 bnd6BC(2,3) = 15750.d0/65536.d0 bnd6BC(2,4) = 15750.d0/65536.d0 bnd6BC(2,5) = -2625.d0/65536.d0 bnd6BC(2,6) = 315.d0/65536.d0 bnd6BC(3,1) = 630.d0/65536.d0 bnd6BC(3,2) = -5250.d0/65536.d0 bnd6BC(3,3) = 31500.d0/65536.d0 bnd6BC(3,4) = 31500.d0/65536.d0 bnd6BC(3,5) = -5250.d0/65536.d0 bnd6BC(3,6) = 630.d0/65536.d0 bnd6BC(4,1) = -210.d0/65536.d0 bnd6BC(4,2) = 1750.d0/65536.d0 bnd6BC(4,3) = -10500.d0/65536.d0 bnd6BC(4,4) = -10500.d0/65536.d0 bnd6BC(4,5) = 1750.d0/65536.d0 bnd6BC(4,6) = -210.d0/65536.d0 bnd6BC(5,1) = 63.d0/65536.d0 bnd6BC(5,2) = -525.d0/65536.d0 bnd6BC(5,3) = 3150.d0/65536.d0 bnd6BC(5,4) = 3150.d0/65536.d0 bnd6BC(5,5) = -525.d0/65536.d0 bnd6BC(5,6) = 63.d0/65536.d0 bnd6BC(6,1) = -9.d0/65536.d0 bnd6BC(6,2) = 75.d0/65536.d0 bnd6BC(6,3) = -450.d0/65536.d0 bnd6BC(6,4) = -450.d0/65536.d0 bnd6BC(6,5) = 75.d0/65536.d0 bnd6BC(6,6) = -9.d0/65536.d0 ! bnd6I(1,1) = 9.d0/65536.d0 bnd6I(1,2) = -75.d0/65536.d0 bnd6I(1,3) = 450.d0/65536.d0 bnd6I(1,4) = 450.d0/65536.d0 bnd6I(1,5) = -75.d0/65536.d0 bnd6I(1,6) = 9.d0/65536.d0 bnd6I(2,1) = -75.d0/65536.d0 bnd6I(2,2) = 625.d0/65536.d0 bnd6I(2,3) = -3750.d0/65536.d0 bnd6I(2,4) = -3750.d0/65536.d0 bnd6I(2,5) = 625.d0/65536.d0 bnd6I(2,6) = -75.d0/65536.d0 bnd6I(3,1) = 450.d0/65536.d0 bnd6I(3,2) = -3750.d0/65536.d0 bnd6I(3,3) = 22500.d0/65536.d0 bnd6I(3,4) = 22500.d0/65536.d0 bnd6I(3,5) = -3750.d0/65536.d0 bnd6I(3,6) = 450.d0/65536.d0 bnd6I(4,1) = 450.d0/65536.d0 bnd6I(4,2) = -3750.d0/65536.d0 bnd6I(4,3) = 22500.d0/65536.d0 bnd6I(4,4) = 22500.d0/65536.d0 bnd6I(4,5) = -3750.d0/65536.d0 bnd6I(4,6) = 450.d0/65536.d0 bnd6I(5,1) = -75.d0/65536.d0 bnd6I(5,2) = 625.d0/65536.d0 bnd6I(5,3) = -3750.d0/65536.d0 bnd6I(5,4) = -3750.d0/65536.d0 bnd6I(5,5) = 625.d0/65536.d0 bnd6I(5,6) = -75.d0/65536.d0 bnd6I(6,1) = 9.d0/65536.d0 bnd6I(6,2) = -75.d0/65536.d0 bnd6I(6,3) = 450.d0/65536.d0 bnd6I(6,4) = 450.d0/65536.d0 bnd6I(6,5) = -75.d0/65536.d0 bnd6I(6,6) = 9.d0/65536.d0 ! ! 8th-order ! bnd8AA(1,1) = 184041.d0/4194304.d0 bnd8AA(1,2) = 1288287.d0/4194304.d0 bnd8AA(1,3) = -1288287.d0/4194304.d0 bnd8AA(1,4) = 1288287.d0/4194304.d0 bnd8AA(1,5) = -920205.d0/4194304.d0 bnd8AA(1,6) = 429429.d0/4194304.d0 bnd8AA(1,7) = -117117.d0/4194304.d0 bnd8AA(1,8) = 14157.d0/4194304.d0 bnd8AA(2,1) = 1288287.d0/4194304.d0 bnd8AA(2,2) = 9018009.d0/4194304.d0 bnd8AA(2,3) = -9018009.d0/4194304.d0 bnd8AA(2,4) = 9018009.d0/4194304.d0 bnd8AA(2,5) = -6441435.d0/4194304.d0 bnd8AA(2,6) = 3006003.d0/4194304.d0 bnd8AA(2,7) = -819819.d0/4194304.d0 bnd8AA(2,8) = 99099.d0/4194304.d0 bnd8AA(3,1) = -1288287.d0/4194304.d0 bnd8AA(3,2) = -9018009.d0/4194304.d0 bnd8AA(3,3) = 9018009.d0/4194304.d0 bnd8AA(3,4) = -9018009.d0/4194304.d0 bnd8AA(3,5) = 6441435.d0/4194304.d0 bnd8AA(3,6) = -3006003.d0/4194304.d0 bnd8AA(3,7) = 819819.d0/4194304.d0 bnd8AA(3,8) = -99099.d0/4194304.d0 bnd8AA(4,1) = 1288287.d0/4194304.d0 bnd8AA(4,2) = 9018009.d0/4194304.d0 bnd8AA(4,3) = -9018009.d0/4194304.d0 bnd8AA(4,4) = 9018009.d0/4194304.d0 bnd8AA(4,5) = -6441435.d0/4194304.d0 bnd8AA(4,6) = 3006003.d0/4194304.d0 bnd8AA(4,7) = -819819.d0/4194304.d0 bnd8AA(4,8) = 99099.d0/4194304.d0 bnd8AA(5,1) = -920205.d0/4194304.d0 bnd8AA(5,2) = -6441435.d0/4194304.d0 bnd8AA(5,3) = 6441435.d0/4194304.d0 bnd8AA(5,4) = -6441435.d0/4194304.d0 bnd8AA(5,5) = 4601025.d0/4194304.d0 bnd8AA(5,6) = -2147145.d0/4194304.d0 bnd8AA(5,7) = 585585.d0/4194304.d0 bnd8AA(5,8) = -70785.d0/4194304.d0 bnd8AA(6,1) = 429429.d0/4194304.d0 bnd8AA(6,2) = 3006003.d0/4194304.d0 bnd8AA(6,3) = -3006003.d0/4194304.d0 bnd8AA(6,4) = 3006003.d0/4194304.d0 bnd8AA(6,5) = -2147145.d0/4194304.d0 bnd8AA(6,6) = 1002001.d0/4194304.d0 bnd8AA(6,7) = -273273.d0/4194304.d0 bnd8AA(6,8) = 33033.d0/4194304.d0 bnd8AA(7,1) = -117117.d0/4194304.d0 bnd8AA(7,2) = -819819.d0/4194304.d0 bnd8AA(7,3) = 819819.d0/4194304.d0 bnd8AA(7,4) = -819819.d0/4194304.d0 bnd8AA(7,5) = 585585.d0/4194304.d0 bnd8AA(7,6) = -273273.d0/4194304.d0 bnd8AA(7,7) = 74529.d0/4194304.d0 bnd8AA(7,8) = -9009.d0/4194304.d0 bnd8AA(8,1) = 14157.d0/4194304.d0 bnd8AA(8,2) = 99099.d0/4194304.d0 bnd8AA(8,3) = -99099.d0/4194304.d0 bnd8AA(8,4) = 99099.d0/4194304.d0 bnd8AA(8,5) = -70785.d0/4194304.d0 bnd8AA(8,6) = 33033.d0/4194304.d0 bnd8AA(8,7) = -9009.d0/4194304.d0 bnd8AA(8,8) = 1089.d0/4194304.d0 ! bnd8AB(1,1) = -14157.d0/4194304.d0 bnd8AB(1,2) = 297297.d0/4194304.d0 bnd8AB(1,3) = 891891.d0/4194304.d0 bnd8AB(1,4) = -495495.d0/4194304.d0 bnd8AB(1,5) = 297297.d0/4194304.d0 bnd8AB(1,6) = -127413.d0/4194304.d0 bnd8AB(1,7) = 33033.d0/4194304.d0 bnd8AB(1,8) = -3861.d0/4194304.d0 bnd8AB(2,1) = -99099.d0/4194304.d0 bnd8AB(2,2) = 2081079.d0/4194304.d0 bnd8AB(2,3) = 6243237.d0/4194304.d0 bnd8AB(2,4) = -3468465.d0/4194304.d0 bnd8AB(2,5) = 2081079.d0/4194304.d0 bnd8AB(2,6) = -891891.d0/4194304.d0 bnd8AB(2,7) = 231231.d0/4194304.d0 bnd8AB(2,8) = -27027.d0/4194304.d0 bnd8AB(3,1) = 99099.d0/4194304.d0 bnd8AB(3,2) = -2081079.d0/4194304.d0 bnd8AB(3,3) = -6243237.d0/4194304.d0 bnd8AB(3,4) = 3468465.d0/4194304.d0 bnd8AB(3,5) = -2081079.d0/4194304.d0 bnd8AB(3,6) = 891891.d0/4194304.d0 bnd8AB(3,7) = -231231.d0/4194304.d0 bnd8AB(3,8) = 27027.d0/4194304.d0 bnd8AB(4,1) = -99099.d0/4194304.d0 bnd8AB(4,2) = 2081079.d0/4194304.d0 bnd8AB(4,3) = 6243237.d0/4194304.d0 bnd8AB(4,4) = -3468465.d0/4194304.d0 bnd8AB(4,5) = 2081079.d0/4194304.d0 bnd8AB(4,6) = -891891.d0/4194304.d0 bnd8AB(4,7) = 231231.d0/4194304.d0 bnd8AB(4,8) = -27027.d0/4194304.d0 bnd8AB(5,1) = 70785.d0/4194304.d0 bnd8AB(5,2) = -1486485.d0/4194304.d0 bnd8AB(5,3) = -4459455.d0/4194304.d0 bnd8AB(5,4) = 2477475.d0/4194304.d0 bnd8AB(5,5) = -1486485.d0/4194304.d0 bnd8AB(5,6) = 637065.d0/4194304.d0 bnd8AB(5,7) = -165165.d0/4194304.d0 bnd8AB(5,8) = 19305.d0/4194304.d0 bnd8AB(6,1) = -33033.d0/4194304.d0 bnd8AB(6,2) = 693693.d0/4194304.d0 bnd8AB(6,3) = 2081079.d0/4194304.d0 bnd8AB(6,4) = -1156155.d0/4194304.d0 bnd8AB(6,5) = 693693.d0/4194304.d0 bnd8AB(6,6) = -297297.d0/4194304.d0 bnd8AB(6,7) = 77077.d0/4194304.d0 bnd8AB(6,8) = -9009.d0/4194304.d0 bnd8AB(7,1) = 9009.d0/4194304.d0 bnd8AB(7,2) = -189189.d0/4194304.d0 bnd8AB(7,3) = -567567.d0/4194304.d0 bnd8AB(7,4) = 315315.d0/4194304.d0 bnd8AB(7,5) = -189189.d0/4194304.d0 bnd8AB(7,6) = 81081.d0/4194304.d0 bnd8AB(7,7) = -21021.d0/4194304.d0 bnd8AB(7,8) = 2457.d0/4194304.d0 bnd8AB(8,1) = -1089.d0/4194304.d0 bnd8AB(8,2) = 22869.d0/4194304.d0 bnd8AB(8,3) = 68607.d0/4194304.d0 bnd8AB(8,4) = -38115.d0/4194304.d0 bnd8AB(8,5) = 22869.d0/4194304.d0 bnd8AB(8,6) = -9801.d0/4194304.d0 bnd8AB(8,7) = 2541.d0/4194304.d0 bnd8AB(8,8) = -297.d0/4194304.d0 ! bnd8AC(1,1) = 3861.d0/4194304.d0 bnd8AC(1,2) = -45045.d0/4194304.d0 bnd8AC(1,3) = 405405.d0/4194304.d0 bnd8AC(1,4) = 675675.d0/4194304.d0 bnd8AC(1,5) = -225225.d0/4194304.d0 bnd8AC(1,6) = 81081.d0/4194304.d0 bnd8AC(1,7) = -19305.d0/4194304.d0 bnd8AC(1,8) = 2145.d0/4194304.d0 bnd8AC(2,1) = 27027.d0/4194304.d0 bnd8AC(2,2) = -315315.d0/4194304.d0 bnd8AC(2,3) = 2837835.d0/4194304.d0 bnd8AC(2,4) = 4729725.d0/4194304.d0 bnd8AC(2,5) = -1576575.d0/4194304.d0 bnd8AC(2,6) = 567567.d0/4194304.d0 bnd8AC(2,7) = -135135.d0/4194304.d0 bnd8AC(2,8) = 15015.d0/4194304.d0 bnd8AC(3,1) = -27027.d0/4194304.d0 bnd8AC(3,2) = 315315.d0/4194304.d0 bnd8AC(3,3) = -2837835.d0/4194304.d0 bnd8AC(3,4) = -4729725.d0/4194304.d0 bnd8AC(3,5) = 1576575.d0/4194304.d0 bnd8AC(3,6) = -567567.d0/4194304.d0 bnd8AC(3,7) = 135135.d0/4194304.d0 bnd8AC(3,8) = -15015.d0/4194304.d0 bnd8AC(4,1) = 27027.d0/4194304.d0 bnd8AC(4,2) = -315315.d0/4194304.d0 bnd8AC(4,3) = 2837835.d0/4194304.d0 bnd8AC(4,4) = 4729725.d0/4194304.d0 bnd8AC(4,5) = -1576575.d0/4194304.d0 bnd8AC(4,6) = 567567.d0/4194304.d0 bnd8AC(4,7) = -135135.d0/4194304.d0 bnd8AC(4,8) = 15015.d0/4194304.d0 bnd8AC(5,1) = -19305.d0/4194304.d0 bnd8AC(5,2) = 225225.d0/4194304.d0 bnd8AC(5,3) = -2027025.d0/4194304.d0 bnd8AC(5,4) = -3378375.d0/4194304.d0 bnd8AC(5,5) = 1126125.d0/4194304.d0 bnd8AC(5,6) = -405405.d0/4194304.d0 bnd8AC(5,7) = 96525.d0/4194304.d0 bnd8AC(5,8) = -10725.d0/4194304.d0 bnd8AC(6,1) = 9009.d0/4194304.d0 bnd8AC(6,2) = -105105.d0/4194304.d0 bnd8AC(6,3) = 945945.d0/4194304.d0 bnd8AC(6,4) = 1576575.d0/4194304.d0 bnd8AC(6,5) = -525525.d0/4194304.d0 bnd8AC(6,6) = 189189.d0/4194304.d0 bnd8AC(6,7) = -45045.d0/4194304.d0 bnd8AC(6,8) = 5005.d0/4194304.d0 bnd8AC(7,1) = -2457.d0/4194304.d0 bnd8AC(7,2) = 28665.d0/4194304.d0 bnd8AC(7,3) = -257985.d0/4194304.d0 bnd8AC(7,4) = -429975.d0/4194304.d0 bnd8AC(7,5) = 143325.d0/4194304.d0 bnd8AC(7,6) = -51597.d0/4194304.d0 bnd8AC(7,7) = 12285.d0/4194304.d0 bnd8AC(7,8) = -1365.d0/4194304.d0 bnd8AC(8,1) = 297.d0/4194304.d0 bnd8AC(8,2) = -3465.d0/4194304.d0 bnd8AC(8,3) = 31185.d0/4194304.d0 bnd8AC(8,4) = 51975.d0/4194304.d0 bnd8AC(8,5) = -17325.d0/4194304.d0 bnd8AC(8,6) = 6237.d0/4194304.d0 bnd8AC(8,7) = -1485.d0/4194304.d0 bnd8AC(8,8) = 165.d0/4194304.d0 ! bnd8AD(1,1) = -2145.d0/4194304.d0 bnd8AD(1,2) = 21021.d0/4194304.d0 bnd8AD(1,3) = -105105.d0/4194304.d0 bnd8AD(1,4) = 525525.d0/4194304.d0 bnd8AD(1,5) = 525525.d0/4194304.d0 bnd8AD(1,6) = -105105.d0/4194304.d0 bnd8AD(1,7) = 21021.d0/4194304.d0 bnd8AD(1,8) = -2145.d0/4194304.d0 bnd8AD(2,1) = -15015.d0/4194304.d0 bnd8AD(2,2) = 147147.d0/4194304.d0 bnd8AD(2,3) = -735735.d0/4194304.d0 bnd8AD(2,4) = 3678675.d0/4194304.d0 bnd8AD(2,5) = 3678675.d0/4194304.d0 bnd8AD(2,6) = -735735.d0/4194304.d0 bnd8AD(2,7) = 147147.d0/4194304.d0 bnd8AD(2,8) = -15015.d0/4194304.d0 bnd8AD(3,1) = 15015.d0/4194304.d0 bnd8AD(3,2) = -147147.d0/4194304.d0 bnd8AD(3,3) = 735735.d0/4194304.d0 bnd8AD(3,4) = -3678675.d0/4194304.d0 bnd8AD(3,5) = -3678675.d0/4194304.d0 bnd8AD(3,6) = 735735.d0/4194304.d0 bnd8AD(3,7) = -147147.d0/4194304.d0 bnd8AD(3,8) = 15015.d0/4194304.d0 bnd8AD(4,1) = -15015.d0/4194304.d0 bnd8AD(4,2) = 147147.d0/4194304.d0 bnd8AD(4,3) = -735735.d0/4194304.d0 bnd8AD(4,4) = 3678675.d0/4194304.d0 bnd8AD(4,5) = 3678675.d0/4194304.d0 bnd8AD(4,6) = -735735.d0/4194304.d0 bnd8AD(4,7) = 147147.d0/4194304.d0 bnd8AD(4,8) = -15015.d0/4194304.d0 bnd8AD(5,1) = 10725.d0/4194304.d0 bnd8AD(5,2) = -105105.d0/4194304.d0 bnd8AD(5,3) = 525525.d0/4194304.d0 bnd8AD(5,4) = -2627625.d0/4194304.d0 bnd8AD(5,5) = -2627625.d0/4194304.d0 bnd8AD(5,6) = 525525.d0/4194304.d0 bnd8AD(5,7) = -105105.d0/4194304.d0 bnd8AD(5,8) = 10725.d0/4194304.d0 bnd8AD(6,1) = -5005.d0/4194304.d0 bnd8AD(6,2) = 49049.d0/4194304.d0 bnd8AD(6,3) = -245245.d0/4194304.d0 bnd8AD(6,4) = 1226225.d0/4194304.d0 bnd8AD(6,5) = 1226225.d0/4194304.d0 bnd8AD(6,6) = -245245.d0/4194304.d0 bnd8AD(6,7) = 49049.d0/4194304.d0 bnd8AD(6,8) = -5005.d0/4194304.d0 bnd8AD(7,1) = 1365.d0/4194304.d0 bnd8AD(7,2) = -13377.d0/4194304.d0 bnd8AD(7,3) = 66885.d0/4194304.d0 bnd8AD(7,4) = -334425.d0/4194304.d0 bnd8AD(7,5) = -334425.d0/4194304.d0 bnd8AD(7,6) = 66885.d0/4194304.d0 bnd8AD(7,7) = -13377.d0/4194304.d0 bnd8AD(7,8) = 1365.d0/4194304.d0 bnd8AD(8,1) = -165.d0/4194304.d0 bnd8AD(8,2) = 1617.d0/4194304.d0 bnd8AD(8,3) = -8085.d0/4194304.d0 bnd8AD(8,4) = 40425.d0/4194304.d0 bnd8AD(8,5) = 40425.d0/4194304.d0 bnd8AD(8,6) = -8085.d0/4194304.d0 bnd8AD(8,7) = 1617.d0/4194304.d0 bnd8AD(8,8) = -165.d0/4194304.d0 ! bnd8BB(1,1) = 1089.d0/4194304.d0 bnd8BB(1,2) = -22869.d0/4194304.d0 bnd8BB(1,3) = -68607.d0/4194304.d0 bnd8BB(1,4) = 38115.d0/4194304.d0 bnd8BB(1,5) = -22869.d0/4194304.d0 bnd8BB(1,6) = 9801.d0/4194304.d0 bnd8BB(1,7) = -2541.d0/4194304.d0 bnd8BB(1,8) = 297.d0/4194304.d0 bnd8BB(2,1) = -22869.d0/4194304.d0 bnd8BB(2,2) = 480249.d0/4194304.d0 bnd8BB(2,3) = 1440747.d0/4194304.d0 bnd8BB(2,4) = -800415.d0/4194304.d0 bnd8BB(2,5) = 480249.d0/4194304.d0 bnd8BB(2,6) = -205821.d0/4194304.d0 bnd8BB(2,7) = 53361.d0/4194304.d0 bnd8BB(2,8) = -6237.d0/4194304.d0 bnd8BB(3,1) = -68607.d0/4194304.d0 bnd8BB(3,2) = 1440747.d0/4194304.d0 bnd8BB(3,3) = 4322241.d0/4194304.d0 bnd8BB(3,4) = -2401245.d0/4194304.d0 bnd8BB(3,5) = 1440747.d0/4194304.d0 bnd8BB(3,6) = -617463.d0/4194304.d0 bnd8BB(3,7) = 160083.d0/4194304.d0 bnd8BB(3,8) = -18711.d0/4194304.d0 bnd8BB(4,1) = 38115.d0/4194304.d0 bnd8BB(4,2) = -800415.d0/4194304.d0 bnd8BB(4,3) = -2401245.d0/4194304.d0 bnd8BB(4,4) = 1334025.d0/4194304.d0 bnd8BB(4,5) = -800415.d0/4194304.d0 bnd8BB(4,6) = 343035.d0/4194304.d0 bnd8BB(4,7) = -88935.d0/4194304.d0 bnd8BB(4,8) = 10395.d0/4194304.d0 bnd8BB(5,1) = -22869.d0/4194304.d0 bnd8BB(5,2) = 480249.d0/4194304.d0 bnd8BB(5,3) = 1440747.d0/4194304.d0 bnd8BB(5,4) = -800415.d0/4194304.d0 bnd8BB(5,5) = 480249.d0/4194304.d0 bnd8BB(5,6) = -205821.d0/4194304.d0 bnd8BB(5,7) = 53361.d0/4194304.d0 bnd8BB(5,8) = -6237.d0/4194304.d0 bnd8BB(6,1) = 9801.d0/4194304.d0 bnd8BB(6,2) = -205821.d0/4194304.d0 bnd8BB(6,3) = -617463.d0/4194304.d0 bnd8BB(6,4) = 343035.d0/4194304.d0 bnd8BB(6,5) = -205821.d0/4194304.d0 bnd8BB(6,6) = 88209.d0/4194304.d0 bnd8BB(6,7) = -22869.d0/4194304.d0 bnd8BB(6,8) = 2673.d0/4194304.d0 bnd8BB(7,1) = -2541.d0/4194304.d0 bnd8BB(7,2) = 53361.d0/4194304.d0 bnd8BB(7,3) = 160083.d0/4194304.d0 bnd8BB(7,4) = -88935.d0/4194304.d0 bnd8BB(7,5) = 53361.d0/4194304.d0 bnd8BB(7,6) = -22869.d0/4194304.d0 bnd8BB(7,7) = 5929.d0/4194304.d0 bnd8BB(7,8) = -693.d0/4194304.d0 bnd8BB(8,1) = 297.d0/4194304.d0 bnd8BB(8,2) = -6237.d0/4194304.d0 bnd8BB(8,3) = -18711.d0/4194304.d0 bnd8BB(8,4) = 10395.d0/4194304.d0 bnd8BB(8,5) = -6237.d0/4194304.d0 bnd8BB(8,6) = 2673.d0/4194304.d0 bnd8BB(8,7) = -693.d0/4194304.d0 bnd8BB(8,8) = 81.d0/4194304.d0 ! bnd8BC(1,1) = -297.d0/4194304.d0 bnd8BC(1,2) = 3465.d0/4194304.d0 bnd8BC(1,3) = -31185.d0/4194304.d0 bnd8BC(1,4) = -51975.d0/4194304.d0 bnd8BC(1,5) = 17325.d0/4194304.d0 bnd8BC(1,6) = -6237.d0/4194304.d0 bnd8BC(1,7) = 1485.d0/4194304.d0 bnd8BC(1,8) = -165.d0/4194304.d0 bnd8BC(2,1) = 6237.d0/4194304.d0 bnd8BC(2,2) = -72765.d0/4194304.d0 bnd8BC(2,3) = 654885.d0/4194304.d0 bnd8BC(2,4) = 1091475.d0/4194304.d0 bnd8BC(2,5) = -363825.d0/4194304.d0 bnd8BC(2,6) = 130977.d0/4194304.d0 bnd8BC(2,7) = -31185.d0/4194304.d0 bnd8BC(2,8) = 3465.d0/4194304.d0 bnd8BC(3,1) = 18711.d0/4194304.d0 bnd8BC(3,2) = -218295.d0/4194304.d0 bnd8BC(3,3) = 1964655.d0/4194304.d0 bnd8BC(3,4) = 3274425.d0/4194304.d0 bnd8BC(3,5) = -1091475.d0/4194304.d0 bnd8BC(3,6) = 392931.d0/4194304.d0 bnd8BC(3,7) = -93555.d0/4194304.d0 bnd8BC(3,8) = 10395.d0/4194304.d0 bnd8BC(4,1) = -10395.d0/4194304.d0 bnd8BC(4,2) = 121275.d0/4194304.d0 bnd8BC(4,3) = -1091475.d0/4194304.d0 bnd8BC(4,4) = -1819125.d0/4194304.d0 bnd8BC(4,5) = 606375.d0/4194304.d0 bnd8BC(4,6) = -218295.d0/4194304.d0 bnd8BC(4,7) = 51975.d0/4194304.d0 bnd8BC(4,8) = -5775.d0/4194304.d0 bnd8BC(5,1) = 6237.d0/4194304.d0 bnd8BC(5,2) = -72765.d0/4194304.d0 bnd8BC(5,3) = 654885.d0/4194304.d0 bnd8BC(5,4) = 1091475.d0/4194304.d0 bnd8BC(5,5) = -363825.d0/4194304.d0 bnd8BC(5,6) = 130977.d0/4194304.d0 bnd8BC(5,7) = -31185.d0/4194304.d0 bnd8BC(5,8) = 3465.d0/4194304.d0 bnd8BC(6,1) = -2673.d0/4194304.d0 bnd8BC(6,2) = 31185.d0/4194304.d0 bnd8BC(6,3) = -280665.d0/4194304.d0 bnd8BC(6,4) = -467775.d0/4194304.d0 bnd8BC(6,5) = 155925.d0/4194304.d0 bnd8BC(6,6) = -56133.d0/4194304.d0 bnd8BC(6,7) = 13365.d0/4194304.d0 bnd8BC(6,8) = -1485.d0/4194304.d0 bnd8BC(7,1) = 693.d0/4194304.d0 bnd8BC(7,2) = -8085.d0/4194304.d0 bnd8BC(7,3) = 72765.d0/4194304.d0 bnd8BC(7,4) = 121275.d0/4194304.d0 bnd8BC(7,5) = -40425.d0/4194304.d0 bnd8BC(7,6) = 14553.d0/4194304.d0 bnd8BC(7,7) = -3465.d0/4194304.d0 bnd8BC(7,8) = 385.d0/4194304.d0 bnd8BC(8,1) = -81.d0/4194304.d0 bnd8BC(8,2) = 945.d0/4194304.d0 bnd8BC(8,3) = -8505.d0/4194304.d0 bnd8BC(8,4) = -14175.d0/4194304.d0 bnd8BC(8,5) = 4725.d0/4194304.d0 bnd8BC(8,6) = -1701.d0/4194304.d0 bnd8BC(8,7) = 405.d0/4194304.d0 bnd8BC(8,8) = -45.d0/4194304.d0 ! bnd8BD(1,1) = 165.d0/4194304.d0 bnd8BD(1,2) = -1617.d0/4194304.d0 bnd8BD(1,3) = 8085.d0/4194304.d0 bnd8BD(1,4) = -40425.d0/4194304.d0 bnd8BD(1,5) = -40425.d0/4194304.d0 bnd8BD(1,6) = 8085.d0/4194304.d0 bnd8BD(1,7) = -1617.d0/4194304.d0 bnd8BD(1,8) = 165.d0/4194304.d0 bnd8BD(2,1) = -3465.d0/4194304.d0 bnd8BD(2,2) = 33957.d0/4194304.d0 bnd8BD(2,3) = -169785.d0/4194304.d0 bnd8BD(2,4) = 848925.d0/4194304.d0 bnd8BD(2,5) = 848925.d0/4194304.d0 bnd8BD(2,6) = -169785.d0/4194304.d0 bnd8BD(2,7) = 33957.d0/4194304.d0 bnd8BD(2,8) = -3465.d0/4194304.d0 bnd8BD(3,1) = -10395.d0/4194304.d0 bnd8BD(3,2) = 101871.d0/4194304.d0 bnd8BD(3,3) = -509355.d0/4194304.d0 bnd8BD(3,4) = 2546775.d0/4194304.d0 bnd8BD(3,5) = 2546775.d0/4194304.d0 bnd8BD(3,6) = -509355.d0/4194304.d0 bnd8BD(3,7) = 101871.d0/4194304.d0 bnd8BD(3,8) = -10395.d0/4194304.d0 bnd8BD(4,1) = 5775.d0/4194304.d0 bnd8BD(4,2) = -56595.d0/4194304.d0 bnd8BD(4,3) = 282975.d0/4194304.d0 bnd8BD(4,4) = -1414875.d0/4194304.d0 bnd8BD(4,5) = -1414875.d0/4194304.d0 bnd8BD(4,6) = 282975.d0/4194304.d0 bnd8BD(4,7) = -56595.d0/4194304.d0 bnd8BD(4,8) = 5775.d0/4194304.d0 bnd8BD(5,1) = -3465.d0/4194304.d0 bnd8BD(5,2) = 33957.d0/4194304.d0 bnd8BD(5,3) = -169785.d0/4194304.d0 bnd8BD(5,4) = 848925.d0/4194304.d0 bnd8BD(5,5) = 848925.d0/4194304.d0 bnd8BD(5,6) = -169785.d0/4194304.d0 bnd8BD(5,7) = 33957.d0/4194304.d0 bnd8BD(5,8) = -3465.d0/4194304.d0 bnd8BD(6,1) = 1485.d0/4194304.d0 bnd8BD(6,2) = -14553.d0/4194304.d0 bnd8BD(6,3) = 72765.d0/4194304.d0 bnd8BD(6,4) = -363825.d0/4194304.d0 bnd8BD(6,5) = -363825.d0/4194304.d0 bnd8BD(6,6) = 72765.d0/4194304.d0 bnd8BD(6,7) = -14553.d0/4194304.d0 bnd8BD(6,8) = 1485.d0/4194304.d0 bnd8BD(7,1) = -385.d0/4194304.d0 bnd8BD(7,2) = 3773.d0/4194304.d0 bnd8BD(7,3) = -18865.d0/4194304.d0 bnd8BD(7,4) = 94325.d0/4194304.d0 bnd8BD(7,5) = 94325.d0/4194304.d0 bnd8BD(7,6) = -18865.d0/4194304.d0 bnd8BD(7,7) = 3773.d0/4194304.d0 bnd8BD(7,8) = -385.d0/4194304.d0 bnd8BD(8,1) = 45.d0/4194304.d0 bnd8BD(8,2) = -441.d0/4194304.d0 bnd8BD(8,3) = 2205.d0/4194304.d0 bnd8BD(8,4) = -11025.d0/4194304.d0 bnd8BD(8,5) = -11025.d0/4194304.d0 bnd8BD(8,6) = 2205.d0/4194304.d0 bnd8BD(8,7) = -441.d0/4194304.d0 bnd8BD(8,8) = 45.d0/4194304.d0 ! bnd8CC(1,1) = 81.d0/4194304.d0 bnd8CC(1,2) = -945.d0/4194304.d0 bnd8CC(1,3) = 8505.d0/4194304.d0 bnd8CC(1,4) = 14175.d0/4194304.d0 bnd8CC(1,5) = -4725.d0/4194304.d0 bnd8CC(1,6) = 1701.d0/4194304.d0 bnd8CC(1,7) = -405.d0/4194304.d0 bnd8CC(1,8) = 45.d0/4194304.d0 bnd8CC(2,1) = -945.d0/4194304.d0 bnd8CC(2,2) = 11025.d0/4194304.d0 bnd8CC(2,3) = -99225.d0/4194304.d0 bnd8CC(2,4) = -165375.d0/4194304.d0 bnd8CC(2,5) = 55125.d0/4194304.d0 bnd8CC(2,6) = -19845.d0/4194304.d0 bnd8CC(2,7) = 4725.d0/4194304.d0 bnd8CC(2,8) = -525.d0/4194304.d0 bnd8CC(3,1) = 8505.d0/4194304.d0 bnd8CC(3,2) = -99225.d0/4194304.d0 bnd8CC(3,3) = 893025.d0/4194304.d0 bnd8CC(3,4) = 1488375.d0/4194304.d0 bnd8CC(3,5) = -496125.d0/4194304.d0 bnd8CC(3,6) = 178605.d0/4194304.d0 bnd8CC(3,7) = -42525.d0/4194304.d0 bnd8CC(3,8) = 4725.d0/4194304.d0 bnd8CC(4,1) = 14175.d0/4194304.d0 bnd8CC(4,2) = -165375.d0/4194304.d0 bnd8CC(4,3) = 1488375.d0/4194304.d0 bnd8CC(4,4) = 2480625.d0/4194304.d0 bnd8CC(4,5) = -826875.d0/4194304.d0 bnd8CC(4,6) = 297675.d0/4194304.d0 bnd8CC(4,7) = -70875.d0/4194304.d0 bnd8CC(4,8) = 7875.d0/4194304.d0 bnd8CC(5,1) = -4725.d0/4194304.d0 bnd8CC(5,2) = 55125.d0/4194304.d0 bnd8CC(5,3) = -496125.d0/4194304.d0 bnd8CC(5,4) = -826875.d0/4194304.d0 bnd8CC(5,5) = 275625.d0/4194304.d0 bnd8CC(5,6) = -99225.d0/4194304.d0 bnd8CC(5,7) = 23625.d0/4194304.d0 bnd8CC(5,8) = -2625.d0/4194304.d0 bnd8CC(6,1) = 1701.d0/4194304.d0 bnd8CC(6,2) = -19845.d0/4194304.d0 bnd8CC(6,3) = 178605.d0/4194304.d0 bnd8CC(6,4) = 297675.d0/4194304.d0 bnd8CC(6,5) = -99225.d0/4194304.d0 bnd8CC(6,6) = 35721.d0/4194304.d0 bnd8CC(6,7) = -8505.d0/4194304.d0 bnd8CC(6,8) = 945.d0/4194304.d0 bnd8CC(7,1) = -405.d0/4194304.d0 bnd8CC(7,2) = 4725.d0/4194304.d0 bnd8CC(7,3) = -42525.d0/4194304.d0 bnd8CC(7,4) = -70875.d0/4194304.d0 bnd8CC(7,5) = 23625.d0/4194304.d0 bnd8CC(7,6) = -8505.d0/4194304.d0 bnd8CC(7,7) = 2025.d0/4194304.d0 bnd8CC(7,8) = -225.d0/4194304.d0 bnd8CC(8,1) = 45.d0/4194304.d0 bnd8CC(8,2) = -525.d0/4194304.d0 bnd8CC(8,3) = 4725.d0/4194304.d0 bnd8CC(8,4) = 7875.d0/4194304.d0 bnd8CC(8,5) = -2625.d0/4194304.d0 bnd8CC(8,6) = 945.d0/4194304.d0 bnd8CC(8,7) = -225.d0/4194304.d0 bnd8CC(8,8) = 25.d0/4194304.d0 ! bnd8CD(1,1) = -45.d0/4194304.d0 bnd8CD(1,2) = 441.d0/4194304.d0 bnd8CD(1,3) = -2205.d0/4194304.d0 bnd8CD(1,4) = 11025.d0/4194304.d0 bnd8CD(1,5) = 11025.d0/4194304.d0 bnd8CD(1,6) = -2205.d0/4194304.d0 bnd8CD(1,7) = 441.d0/4194304.d0 bnd8CD(1,8) = -45.d0/4194304.d0 bnd8CD(2,1) = 525.d0/4194304.d0 bnd8CD(2,2) = -5145.d0/4194304.d0 bnd8CD(2,3) = 25725.d0/4194304.d0 bnd8CD(2,4) = -128625.d0/4194304.d0 bnd8CD(2,5) = -128625.d0/4194304.d0 bnd8CD(2,6) = 25725.d0/4194304.d0 bnd8CD(2,7) = -5145.d0/4194304.d0 bnd8CD(2,8) = 525.d0/4194304.d0 bnd8CD(3,1) = -4725.d0/4194304.d0 bnd8CD(3,2) = 46305.d0/4194304.d0 bnd8CD(3,3) = -231525.d0/4194304.d0 bnd8CD(3,4) = 1157625.d0/4194304.d0 bnd8CD(3,5) = 1157625.d0/4194304.d0 bnd8CD(3,6) = -231525.d0/4194304.d0 bnd8CD(3,7) = 46305.d0/4194304.d0 bnd8CD(3,8) = -4725.d0/4194304.d0 bnd8CD(4,1) = -7875.d0/4194304.d0 bnd8CD(4,2) = 77175.d0/4194304.d0 bnd8CD(4,3) = -385875.d0/4194304.d0 bnd8CD(4,4) = 1929375.d0/4194304.d0 bnd8CD(4,5) = 1929375.d0/4194304.d0 bnd8CD(4,6) = -385875.d0/4194304.d0 bnd8CD(4,7) = 77175.d0/4194304.d0 bnd8CD(4,8) = -7875.d0/4194304.d0 bnd8CD(5,1) = 2625.d0/4194304.d0 bnd8CD(5,2) = -25725.d0/4194304.d0 bnd8CD(5,3) = 128625.d0/4194304.d0 bnd8CD(5,4) = -643125.d0/4194304.d0 bnd8CD(5,5) = -643125.d0/4194304.d0 bnd8CD(5,6) = 128625.d0/4194304.d0 bnd8CD(5,7) = -25725.d0/4194304.d0 bnd8CD(5,8) = 2625.d0/4194304.d0 bnd8CD(6,1) = -945.d0/4194304.d0 bnd8CD(6,2) = 9261.d0/4194304.d0 bnd8CD(6,3) = -46305.d0/4194304.d0 bnd8CD(6,4) = 231525.d0/4194304.d0 bnd8CD(6,5) = 231525.d0/4194304.d0 bnd8CD(6,6) = -46305.d0/4194304.d0 bnd8CD(6,7) = 9261.d0/4194304.d0 bnd8CD(6,8) = -945.d0/4194304.d0 bnd8CD(7,1) = 225.d0/4194304.d0 bnd8CD(7,2) = -2205.d0/4194304.d0 bnd8CD(7,3) = 11025.d0/4194304.d0 bnd8CD(7,4) = -55125.d0/4194304.d0 bnd8CD(7,5) = -55125.d0/4194304.d0 bnd8CD(7,6) = 11025.d0/4194304.d0 bnd8CD(7,7) = -2205.d0/4194304.d0 bnd8CD(7,8) = 225.d0/4194304.d0 bnd8CD(8,1) = -25.d0/4194304.d0 bnd8CD(8,2) = 245.d0/4194304.d0 bnd8CD(8,3) = -1225.d0/4194304.d0 bnd8CD(8,4) = 6125.d0/4194304.d0 bnd8CD(8,5) = 6125.d0/4194304.d0 bnd8CD(8,6) = -1225.d0/4194304.d0 bnd8CD(8,7) = 245.d0/4194304.d0 bnd8CD(8,8) = -25.d0/4194304.d0 ! bnd8I(1,1) = 25.d0/4194304.d0 bnd8I(1,2) = -245.d0/4194304.d0 bnd8I(1,3) = 1225.d0/4194304.d0 bnd8I(1,4) = -6125.d0/4194304.d0 bnd8I(1,5) = -6125.d0/4194304.d0 bnd8I(1,6) = 1225.d0/4194304.d0 bnd8I(1,7) = -245.d0/4194304.d0 bnd8I(1,8) = 25.d0/4194304.d0 bnd8I(2,1) = -245.d0/4194304.d0 bnd8I(2,2) = 2401.d0/4194304.d0 bnd8I(2,3) = -12005.d0/4194304.d0 bnd8I(2,4) = 60025.d0/4194304.d0 bnd8I(2,5) = 60025.d0/4194304.d0 bnd8I(2,6) = -12005.d0/4194304.d0 bnd8I(2,7) = 2401.d0/4194304.d0 bnd8I(2,8) = -245.d0/4194304.d0 bnd8I(3,1) = 1225.d0/4194304.d0 bnd8I(3,2) = -12005.d0/4194304.d0 bnd8I(3,3) = 60025.d0/4194304.d0 bnd8I(3,4) = -300125.d0/4194304.d0 bnd8I(3,5) = -300125.d0/4194304.d0 bnd8I(3,6) = 60025.d0/4194304.d0 bnd8I(3,7) = -12005.d0/4194304.d0 bnd8I(3,8) = 1225.d0/4194304.d0 bnd8I(4,1) = -6125.d0/4194304.d0 bnd8I(4,2) = 60025.d0/4194304.d0 bnd8I(4,3) = -300125.d0/4194304.d0 bnd8I(4,4) = 1500625.d0/4194304.d0 bnd8I(4,5) = 1500625.d0/4194304.d0 bnd8I(4,6) = -300125.d0/4194304.d0 bnd8I(4,7) = 60025.d0/4194304.d0 bnd8I(4,8) = -6125.d0/4194304.d0 bnd8I(5,1) = -6125.d0/4194304.d0 bnd8I(5,2) = 60025.d0/4194304.d0 bnd8I(5,3) = -300125.d0/4194304.d0 bnd8I(5,4) = 1500625.d0/4194304.d0 bnd8I(5,5) = 1500625.d0/4194304.d0 bnd8I(5,6) = -300125.d0/4194304.d0 bnd8I(5,7) = 60025.d0/4194304.d0 bnd8I(5,8) = -6125.d0/4194304.d0 bnd8I(6,1) = 1225.d0/4194304.d0 bnd8I(6,2) = -12005.d0/4194304.d0 bnd8I(6,3) = 60025.d0/4194304.d0 bnd8I(6,4) = -300125.d0/4194304.d0 bnd8I(6,5) = -300125.d0/4194304.d0 bnd8I(6,6) = 60025.d0/4194304.d0 bnd8I(6,7) = -12005.d0/4194304.d0 bnd8I(6,8) = 1225.d0/4194304.d0 bnd8I(7,1) = -245.d0/4194304.d0 bnd8I(7,2) = 2401.d0/4194304.d0 bnd8I(7,3) = -12005.d0/4194304.d0 bnd8I(7,4) = 60025.d0/4194304.d0 bnd8I(7,5) = 60025.d0/4194304.d0 bnd8I(7,6) = -12005.d0/4194304.d0 bnd8I(7,7) = 2401.d0/4194304.d0 bnd8I(7,8) = -245.d0/4194304.d0 bnd8I(8,1) = 25.d0/4194304.d0 bnd8I(8,2) = -245.d0/4194304.d0 bnd8I(8,3) = 1225.d0/4194304.d0 bnd8I(8,4) = -6125.d0/4194304.d0 bnd8I(8,5) = -6125.d0/4194304.d0 bnd8I(8,6) = 1225.d0/4194304.d0 bnd8I(8,7) = -245.d0/4194304.d0 bnd8I(8,8) = 25.d0/4194304.d0 ! bnd10I(1,1) = 1225.d0/4294967296.d0 bnd10I(1,2) = -14175.d0/4294967296.d0 bnd10I(1,3) = 79380.d0/4294967296.d0 bnd10I(1,4) = -308700.d0/4294967296.d0 bnd10I(1,5) = 1389150.d0/4294967296.d0 bnd10I(1,6) = 1389150.d0/4294967296.d0 bnd10I(1,7) = -308700.d0/4294967296.d0 bnd10I(1,8) = 79380.d0/4294967296.d0 bnd10I(1,9) = -14175.d0/4294967296.d0 bnd10I(1,10) = 1225.d0/4294967296.d0 bnd10I(2,1) = -14175.d0/4294967296.d0 bnd10I(2,2) = 164025.d0/4294967296.d0 bnd10I(2,3) = -918540.d0/4294967296.d0 bnd10I(2,4) = 3572100.d0/4294967296.d0 bnd10I(2,5) = -16074450.d0/4294967296.d0 bnd10I(2,6) = -16074450.d0/4294967296.d0 bnd10I(2,7) = 3572100.d0/4294967296.d0 bnd10I(2,8) = -918540.d0/4294967296.d0 bnd10I(2,9) = 164025.d0/4294967296.d0 bnd10I(2,10) = -14175.d0/4294967296.d0 bnd10I(3,1) = 79380.d0/4294967296.d0 bnd10I(3,2) = -918540.d0/4294967296.d0 bnd10I(3,3) = 5143824.d0/4294967296.d0 bnd10I(3,4) = -20003760.d0/4294967296.d0 bnd10I(3,5) = 90016920.d0/4294967296.d0 bnd10I(3,6) = 90016920.d0/4294967296.d0 bnd10I(3,7) = -20003760.d0/4294967296.d0 bnd10I(3,8) = 5143824.d0/4294967296.d0 bnd10I(3,9) = -918540.d0/4294967296.d0 bnd10I(3,10) = 79380.d0/4294967296.d0 bnd10I(4,1) = -308700.d0/4294967296.d0 bnd10I(4,2) = 3572100.d0/4294967296.d0 bnd10I(4,3) = -20003760.d0/4294967296.d0 bnd10I(4,4) = 77792400.d0/4294967296.d0 bnd10I(4,5) = -350065800.d0/4294967296.d0 bnd10I(4,6) = -350065800.d0/4294967296.d0 bnd10I(4,7) = 77792400.d0/4294967296.d0 bnd10I(4,8) = -20003760.d0/4294967296.d0 bnd10I(4,9) = 3572100.d0/4294967296.d0 bnd10I(4,10) = -308700.d0/4294967296.d0 bnd10I(5,1) = 1389150.d0/4294967296.d0 bnd10I(5,2) = -16074450.d0/4294967296.d0 bnd10I(5,3) = 90016920.d0/4294967296.d0 bnd10I(5,4) = -350065800.d0/4294967296.d0 bnd10I(5,5) = 1575296100.d0/4294967296.d0 bnd10I(5,6) = 1575296100.d0/4294967296.d0 bnd10I(5,7) = -350065800.d0/4294967296.d0 bnd10I(5,8) = 90016920.d0/4294967296.d0 bnd10I(5,9) = -16074450.d0/4294967296.d0 bnd10I(5,10) = 1389150.d0/4294967296.d0 bnd10I(6,1) = 1389150.d0/4294967296.d0 bnd10I(6,2) = -16074450.d0/4294967296.d0 bnd10I(6,3) = 90016920.d0/4294967296.d0 bnd10I(6,4) = -350065800.d0/4294967296.d0 bnd10I(6,5) = 1575296100.d0/4294967296.d0 bnd10I(6,6) = 1575296100.d0/4294967296.d0 bnd10I(6,7) = -350065800.d0/4294967296.d0 bnd10I(6,8) = 90016920.d0/4294967296.d0 bnd10I(6,9) = -16074450.d0/4294967296.d0 bnd10I(6,10) = 1389150.d0/4294967296.d0 bnd10I(7,1) = -308700.d0/4294967296.d0 bnd10I(7,2) = 3572100.d0/4294967296.d0 bnd10I(7,3) = -20003760.d0/4294967296.d0 bnd10I(7,4) = 77792400.d0/4294967296.d0 bnd10I(7,5) = -350065800.d0/4294967296.d0 bnd10I(7,6) = -350065800.d0/4294967296.d0 bnd10I(7,7) = 77792400.d0/4294967296.d0 bnd10I(7,8) = -20003760.d0/4294967296.d0 bnd10I(7,9) = 3572100.d0/4294967296.d0 bnd10I(7,10) = -308700.d0/4294967296.d0 bnd10I(8,1) = 79380.d0/4294967296.d0 bnd10I(8,2) = -918540.d0/4294967296.d0 bnd10I(8,3) = 5143824.d0/4294967296.d0 bnd10I(8,4) = -20003760.d0/4294967296.d0 bnd10I(8,5) = 90016920.d0/4294967296.d0 bnd10I(8,6) = 90016920.d0/4294967296.d0 bnd10I(8,7) = -20003760.d0/4294967296.d0 bnd10I(8,8) = 5143824.d0/4294967296.d0 bnd10I(8,9) = -918540.d0/4294967296.d0 bnd10I(8,10) = 79380.d0/4294967296.d0 bnd10I(9,1) = -14175.d0/4294967296.d0 bnd10I(9,2) = 164025.d0/4294967296.d0 bnd10I(9,3) = -918540.d0/4294967296.d0 bnd10I(9,4) = 3572100.d0/4294967296.d0 bnd10I(9,5) = -16074450.d0/4294967296.d0 bnd10I(9,6) = -16074450.d0/4294967296.d0 bnd10I(9,7) = 3572100.d0/4294967296.d0 bnd10I(9,8) = -918540.d0/4294967296.d0 bnd10I(9,9) = 164025.d0/4294967296.d0 bnd10I(9,10) = -14175.d0/4294967296.d0 bnd10I(10,1) = 1225.d0/4294967296.d0 bnd10I(10,2) = -14175.d0/4294967296.d0 bnd10I(10,3) = 79380.d0/4294967296.d0 bnd10I(10,4) = -308700.d0/4294967296.d0 bnd10I(10,5) = 1389150.d0/4294967296.d0 bnd10I(10,6) = 1389150.d0/4294967296.d0 bnd10I(10,7) = -308700.d0/4294967296.d0 bnd10I(10,8) = 79380.d0/4294967296.d0 bnd10I(10,9) = -14175.d0/4294967296.d0 bnd10I(10,10) = 1225.d0/4294967296.d0 ! return end !====================================================================