TABLE OF CONTENTS
ABINIT/contstr25 [ Functions ]
NAME
contstr25
FUNCTION
Carries out specialized metric tensor operations needed for contraction of the 2nd strain derivative of the l=0,1,2,3 nonlocal Kleinman-Bylander pseudopotential operation. Derivatives are wrt a pair of cartesian strain components. Full advantage is taken of the full permutational symmetry of these tensors.
COPYRIGHT
Copyright (C) 1998-2022 ABINIT group (DRH) This file is distributed under the terms of the GNU General Public License, see ~abinit/COPYING or http://www.gnu.org/copyleft/gpl.txt . For the initials of contributors, see ~abinit/doc/developers/contributors.txt.
INPUTS
istr1=1,...6 specifies cartesian strain component 11,22,33,32,31,21 istr2=seccond strain component rank=angular momentum gm(3,3)=metric tensor (array is symmetric but stored as 3x3) gprimd(3,3)=reciprocal space dimensional primitive translations aa(2,*)=unique elements of complex right-hand tensor bb(2,*)=unique elements of complex left-hand tensor
OUTPUT
e2nl=contraction for nonlocal 2nd-order strain derivative energy
NOTES
All tensors are stored in a compressed storage mode defined below; input and output conform to this scheme. When tensor elements occur repeatedly due to symmetry, the WEIGHT IS INCLUDED in the output tensor element to simplify later contractions with other tensors of the same rank and form, i.e. the next contraction is then simply a dot product over the unique elements.
SOURCE
42 #if defined HAVE_CONFIG_H 43 #include "config.h" 44 #endif 45 46 #include "abi_common.h" 47 48 49 module m_contstr25 50 contains 51 !!** 52 53 54 55 subroutine contstr25(istr1,istr2,rank,gm,gprimd,e2nl,aa,bb) 56 57 use defs_basis 58 use m_abicore 59 implicit none 60 61 !Arguments ------------------------------------ 62 !scalars 63 integer,intent(in) :: istr1,istr2,rank 64 real(dp),intent(out) :: e2nl 65 !arrays 66 real(dp),intent(in) :: aa(2,((rank+3)*(rank+4))/2),bb(2,((rank+3)*(rank+4))/2) 67 real(dp),intent(in) :: gm(3,3),gprimd(3,3) 68 69 !Local variables------------------------------- 70 !scalars 71 integer,parameter :: mrank=3 72 integer :: ii,jj,ka,kb,kd,kg 73 !arrays 74 integer,save :: idx(12)=(/1,1,2,2,3,3,3,2,3,1,2,1/) 75 real(dp) :: d2gm(3,3),dgm01(3,3),dgm10(3,3),tmp(2) 76 real(dp),allocatable :: cm(:,:) 77 78 ! ************************************************************************* 79 ABI_MALLOC(cm,(((mrank+3)*(mrank+4))/2,((mrank+3)*(mrank+4))/2)) 80 81 ka=idx(2*istr1-1);kb=idx(2*istr1);kg=idx(2*istr2-1);kd=idx(2*istr2) 82 83 do ii = 1,3 84 dgm01(:,ii)=-(gprimd(ka,:)*gprimd(kb,ii)+gprimd(kb,:)*gprimd(ka,ii)) 85 dgm10(:,ii)=-(gprimd(kg,:)*gprimd(kd,ii)+gprimd(kd,:)*gprimd(kg,ii)) 86 end do 87 88 d2gm(:,:)=0.d0 89 do ii = 1,3 90 if(ka==kg) d2gm(:,ii)=d2gm(:,ii)& 91 & +gprimd(kb,:)*gprimd(kd,ii)+gprimd(kd,:)*gprimd(kb,ii) 92 if(ka==kd) d2gm(:,ii)=d2gm(:,ii)& 93 & +gprimd(kb,:)*gprimd(kg,ii)+gprimd(kg,:)*gprimd(kb,ii) 94 if(kb==kg) d2gm(:,ii)=d2gm(:,ii)& 95 & +gprimd(ka,:)*gprimd(kd,ii)+gprimd(kd,:)*gprimd(ka,ii) 96 if(kb==kd) d2gm(:,ii)=d2gm(:,ii)& 97 & +gprimd(ka,:)*gprimd(kg,ii)+gprimd(kg,:)*gprimd(ka,ii) 98 end do 99 d2gm(:,:)=0.5d0*d2gm(:,:) 100 101 ! 102 !The code below was written by a Mathematica program and formatted by 103 !a combination of editing scripts. It is not intended to be read 104 !by human beings, and certainly not to be modified by one. Conceivably 105 !it could be shortened somewhat by identifying common subexpressions. 106 ! 107 108 if(rank==0)then 109 cm(1,1)=2*dgm01(1,1)*dgm10(1,1) 110 cm(2,1)=dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2) 111 cm(3,1)=dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3) 112 cm(4,1)=2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)) 113 cm(5,1)=2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3)) 114 cm(6,1)=2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2)) 115 cm(1,2)=dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2) 116 cm(2,2)=2*dgm01(2,2)*dgm10(2,2) 117 cm(3,2)=dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3) 118 cm(4,2)=2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)) 119 cm(5,2)=2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)) 120 cm(6,2)=2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 121 cm(1,3)=dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3) 122 cm(2,3)=dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3) 123 cm(3,3)=2*dgm01(3,3)*dgm10(3,3) 124 cm(4,3)=2*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 125 cm(5,3)=2*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 126 cm(6,3)=2*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 127 cm(1,4)=2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)) 128 cm(2,4)=2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)) 129 cm(3,4)=2*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 130 cm(4,4)=8*dgm01(2,3)*dgm10(2,3) 131 cm(5,4)=4*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3)) 132 cm(6,4)=4*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3)) 133 cm(1,5)=2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3)) 134 cm(2,5)=2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)) 135 cm(3,5)=2*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 136 cm(4,5)=4*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3)) 137 cm(5,5)=8*dgm01(1,3)*dgm10(1,3) 138 cm(6,5)=4*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3)) 139 cm(1,6)=2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2)) 140 cm(2,6)=2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 141 cm(3,6)=2*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 142 cm(4,6)=4*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3)) 143 cm(5,6)=4*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3)) 144 cm(6,6)=8*dgm01(1,2)*dgm10(1,2) 145 elseif(rank==1)then 146 cm(1,1)=2*gm(1,1)*dgm01(1,1)*dgm10(1,1) 147 cm(2,1)=2*gm(1,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 148 & +gm(1,1)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)) 149 cm(3,1)=2*gm(1,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 150 & +gm(1,1)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 151 cm(4,1)=2*(gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 152 & +gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+gm(1,1)& 153 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))) 154 cm(5,1)=2*(gm(1,3)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(dgm01(1,3)& 155 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))) 156 cm(6,1)=2*(gm(1,2)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(dgm01(1,2)& 157 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))) 158 cm(7,1)=gm(1,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)) 159 cm(8,1)=2*gm(1,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 160 & +gm(1,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 161 cm(9,1)=gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 162 & +2*gm(1,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)) 163 cm(10,1)=gm(1,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 164 cm(1,2)=2*gm(1,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 165 & +gm(1,1)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)) 166 cm(2,2)=2*(4*gm(2,2)*dgm01(1,2)*dgm10(1,2)+gm(1,1)*dgm01(2,2)& 167 & *dgm10(2,2)+2*gm(1,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 168 cm(3,2)=2*gm(1,3)*dgm01(2,2)*dgm10(1,3)+4*gm(2,3)*(dgm01(1,3)& 169 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+2*gm(1,3)*dgm01(1,3)*dgm10(2,2)& 170 & +gm(1,1)*dgm01(3,3)*dgm10(2,2)+gm(1,1)*dgm01(2,2)*dgm10(3,3)& 171 & +2*gm(1,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 172 cm(4,2)=2*(4*gm(2,3)*dgm01(1,2)*dgm10(1,2)+gm(1,3)*dgm01(2,2)& 173 & *dgm10(1,2)+2*gm(1,2)*dgm01(2,3)*dgm10(1,2)+gm(1,2)*dgm01(2,2)& 174 & *dgm10(1,3)+2*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 175 & +gm(1,3)*dgm01(1,2)*dgm10(2,2)+gm(1,2)*dgm01(1,3)*dgm10(2,2)& 176 & +gm(1,1)*dgm01(2,3)*dgm10(2,2)+2*gm(1,2)*dgm01(1,2)*dgm10(2,3)& 177 & +gm(1,1)*dgm01(2,2)*dgm10(2,3)) 178 cm(5,2)=2*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 179 & +gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+2*(2*gm(1,2)& 180 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+gm(1,1)*(dgm01(2,2)& 181 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))) 182 cm(6,2)=2*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 183 & +gm(1,2)*(dgm01(2,2)*dgm10(1,1)+8*dgm01(1,2)*dgm10(1,2)+dgm01(1,1)& 184 & *dgm10(2,2))+2*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 185 cm(7,2)=2*(gm(1,2)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(dgm01(2,2)& 186 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 187 cm(8,2)=2*gm(1,3)*dgm01(2,3)*dgm10(2,2)+gm(1,2)*dgm01(3,3)*dgm10(2,2)& 188 & +2*gm(1,3)*dgm01(2,2)*dgm10(2,3)+4*gm(2,3)*(dgm01(2,3)*dgm10(1,2)& 189 & +dgm01(1,2)*dgm10(2,3))+gm(1,2)*dgm01(2,2)*dgm10(3,3)+2*gm(2,2)& 190 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 191 cm(9,2)=2*(gm(1,3)*dgm01(2,2)*dgm10(2,2)+gm(1,2)*dgm01(2,3)& 192 & *dgm10(2,2)+gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 193 & +gm(1,2)*dgm01(2,2)*dgm10(2,3)+2*gm(2,2)*(dgm01(2,3)*dgm10(1,2)& 194 & +dgm01(1,2)*dgm10(2,3))) 195 cm(10,2)=2*gm(2,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 196 & +gm(1,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 197 cm(1,3)=2*gm(1,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 198 & +gm(1,1)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 199 cm(2,3)=2*gm(1,3)*dgm01(2,2)*dgm10(1,3)+4*gm(2,3)*(dgm01(1,3)& 200 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+2*gm(1,3)*dgm01(1,3)*dgm10(2,2)& 201 & +gm(1,1)*dgm01(3,3)*dgm10(2,2)+gm(1,1)*dgm01(2,2)*dgm10(3,3)& 202 & +2*gm(1,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 203 cm(3,3)=2*(4*gm(3,3)*dgm01(1,3)*dgm10(1,3)+gm(1,1)*dgm01(3,3)& 204 & *dgm10(3,3)+2*gm(1,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 205 cm(4,3)=2*(4*gm(2,3)*dgm01(1,3)*dgm10(1,3)+gm(1,2)*dgm01(3,3)& 206 & *dgm10(1,3)+2*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 207 & +gm(1,1)*dgm01(3,3)*dgm10(2,3)+gm(1,2)*dgm01(1,3)*dgm10(3,3)& 208 & +gm(1,1)*dgm01(2,3)*dgm10(3,3)+gm(1,3)*(dgm01(3,3)*dgm10(1,2)& 209 & +2*dgm01(2,3)*dgm10(1,3)+2*dgm01(1,3)*dgm10(2,3)+dgm01(1,2)*dgm10(3,3))) 210 cm(5,3)=2*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 211 & +gm(1,3)*(dgm01(3,3)*dgm10(1,1)+8*dgm01(1,3)*dgm10(1,3)+dgm01(1,1)& 212 & *dgm10(3,3))+2*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 213 cm(6,3)=2*gm(2,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 214 & +gm(1,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 215 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+gm(1,1)*(dgm01(3,3)& 216 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))) 217 cm(7,3)=2*gm(2,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 218 & +gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 219 cm(8,3)=2*(gm(1,3)*dgm01(3,3)*dgm10(2,3)+2*gm(3,3)*(dgm01(2,3)& 220 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,3)*dgm01(2,3)*dgm10(3,3)& 221 & +gm(1,2)*dgm01(3,3)*dgm10(3,3)+gm(2,3)*(dgm01(3,3)*dgm10(1,3)& 222 & +dgm01(1,3)*dgm10(3,3))) 223 cm(9,3)=gm(1,3)*dgm01(3,3)*dgm10(2,2)+2*gm(3,3)*(dgm01(2,2)& 224 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+2*gm(1,2)*dgm01(3,3)*dgm10(2,3)& 225 & +4*gm(2,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,3)& 226 & *dgm01(2,2)*dgm10(3,3)+2*gm(1,2)*dgm01(2,3)*dgm10(3,3) 227 cm(10,3)=2*(gm(1,3)*dgm01(3,3)*dgm10(3,3)+gm(3,3)*(dgm01(3,3)& 228 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 229 cm(1,4)=2*(gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 230 & +gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+gm(1,1)& 231 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))) 232 cm(2,4)=2*(4*gm(2,3)*dgm01(1,2)*dgm10(1,2)+gm(1,3)*dgm01(2,2)& 233 & *dgm10(1,2)+2*gm(1,2)*dgm01(2,3)*dgm10(1,2)+gm(1,2)*dgm01(2,2)& 234 & *dgm10(1,3)+2*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 235 & +gm(1,3)*dgm01(1,2)*dgm10(2,2)+gm(1,2)*dgm01(1,3)*dgm10(2,2)& 236 & +gm(1,1)*dgm01(2,3)*dgm10(2,2)+2*gm(1,2)*dgm01(1,2)*dgm10(2,3)& 237 & +gm(1,1)*dgm01(2,2)*dgm10(2,3)) 238 cm(3,4)=2*(4*gm(2,3)*dgm01(1,3)*dgm10(1,3)+gm(1,2)*dgm01(3,3)& 239 & *dgm10(1,3)+2*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 240 & +gm(1,1)*dgm01(3,3)*dgm10(2,3)+gm(1,2)*dgm01(1,3)*dgm10(3,3)& 241 & +gm(1,1)*dgm01(2,3)*dgm10(3,3)+gm(1,3)*(dgm01(3,3)*dgm10(1,2)& 242 & +2*dgm01(2,3)*dgm10(1,3)+2*dgm01(1,3)*dgm10(2,3)+dgm01(1,2)*dgm10(3,3))) 243 cm(4,4)=8*(gm(3,3)*dgm01(1,2)*dgm10(1,2)+gm(1,3)*dgm01(2,3)& 244 & *dgm10(1,2)+gm(2,2)*dgm01(1,3)*dgm10(1,3)+gm(1,2)*dgm01(2,3)& 245 & *dgm10(1,3)+gm(2,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 246 & +gm(1,3)*dgm01(1,2)*dgm10(2,3)+gm(1,2)*dgm01(1,3)*dgm10(2,3)& 247 & +gm(1,1)*dgm01(2,3)*dgm10(2,3)) 248 cm(5,4)=2*(gm(1,3)*dgm01(2,3)*dgm10(1,1)+2*gm(1,3)*dgm01(1,3)& 249 & *dgm10(1,2)+gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 250 & +2*gm(1,3)*dgm01(1,2)*dgm10(1,3)+4*gm(1,2)*dgm01(1,3)*dgm10(1,3)& 251 & +2*gm(1,1)*dgm01(2,3)*dgm10(1,3)+gm(2,3)*(dgm01(1,3)*dgm10(1,1)& 252 & +dgm01(1,1)*dgm10(1,3))+gm(1,3)*dgm01(1,1)*dgm10(2,3)+2*gm(1,1)& 253 & *dgm01(1,3)*dgm10(2,3)) 254 cm(6,4)=2*(gm(1,2)*dgm01(2,3)*dgm10(1,1)+4*gm(1,3)*dgm01(1,2)& 255 & *dgm10(1,2)+2*gm(1,2)*dgm01(1,3)*dgm10(1,2)+2*gm(1,1)*dgm01(2,3)& 256 & *dgm10(1,2)+gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 257 & +2*gm(1,2)*dgm01(1,2)*dgm10(1,3)+gm(2,2)*(dgm01(1,3)*dgm10(1,1)& 258 & +dgm01(1,1)*dgm10(1,3))+gm(1,2)*dgm01(1,1)*dgm10(2,3)+2*gm(1,1)& 259 & *dgm01(1,2)*dgm10(2,3)) 260 cm(7,4)=2*(gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 261 & +gm(2,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+gm(1,2)& 262 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 263 cm(8,4)=2*(gm(2,2)*dgm01(3,3)*dgm10(1,3)+4*gm(1,3)*dgm01(2,3)& 264 & *dgm10(2,3)+gm(1,2)*dgm01(3,3)*dgm10(2,3)+2*gm(3,3)*(dgm01(2,3)& 265 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+gm(2,2)*dgm01(1,3)*dgm10(3,3)& 266 & +gm(1,2)*dgm01(2,3)*dgm10(3,3)+gm(2,3)*(dgm01(3,3)*dgm10(1,2)& 267 & +2*dgm01(2,3)*dgm10(1,3)+2*dgm01(1,3)*dgm10(2,3)+dgm01(1,2)*dgm10(3,3))) 268 cm(9,4)=2*(2*gm(2,2)*dgm01(2,3)*dgm10(1,3)+gm(1,3)*dgm01(2,3)& 269 & *dgm10(2,2)+gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 270 & +2*gm(2,2)*dgm01(1,3)*dgm10(2,3)+gm(1,3)*dgm01(2,2)*dgm10(2,3)& 271 & +4*gm(1,2)*dgm01(2,3)*dgm10(2,3)+gm(2,3)*(2*dgm01(2,3)*dgm10(1,2)& 272 & +dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)+2*dgm01(1,2)*dgm10(2,3))) 273 cm(10,4)=2*(gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 274 & +gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+gm(1,3)& 275 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 276 cm(1,5)=2*(gm(1,3)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(dgm01(1,3)& 277 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))) 278 cm(2,5)=2*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 279 & +gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+2*(2*gm(1,2)& 280 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+gm(1,1)*(dgm01(2,2)& 281 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))) 282 cm(3,5)=2*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 283 & +gm(1,3)*(dgm01(3,3)*dgm10(1,1)+8*dgm01(1,3)*dgm10(1,3)+dgm01(1,1)& 284 & *dgm10(3,3))+2*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 285 cm(4,5)=2*(gm(1,3)*dgm01(2,3)*dgm10(1,1)+2*gm(1,3)*dgm01(1,3)& 286 & *dgm10(1,2)+gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 287 & +2*gm(1,3)*dgm01(1,2)*dgm10(1,3)+4*gm(1,2)*dgm01(1,3)*dgm10(1,3)& 288 & +2*gm(1,1)*dgm01(2,3)*dgm10(1,3)+gm(2,3)*(dgm01(1,3)*dgm10(1,1)& 289 & +dgm01(1,1)*dgm10(1,3))+gm(1,3)*dgm01(1,1)*dgm10(2,3)+2*gm(1,1)& 290 & *dgm01(1,3)*dgm10(2,3)) 291 cm(5,5)=2*(gm(3,3)*dgm01(1,1)*dgm10(1,1)+2*(2*gm(1,1)*dgm01(1,3)& 292 & *dgm10(1,3)+gm(1,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3)))) 293 cm(6,5)=2*(gm(2,3)*dgm01(1,1)*dgm10(1,1)+gm(1,2)*dgm01(1,3)& 294 & *dgm10(1,1)+2*gm(1,1)*dgm01(1,3)*dgm10(1,2)+gm(1,3)*(dgm01(1,2)& 295 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+gm(1,2)*dgm01(1,1)*dgm10(1,3)& 296 & +2*gm(1,1)*dgm01(1,2)*dgm10(1,3)) 297 cm(7,5)=gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 298 & +2*gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)) 299 cm(8,5)=2*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 300 & +gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 301 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,2)*(dgm01(3,3)& 302 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 303 cm(9,5)=gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 304 & +2*(gm(1,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+gm(2,3)& 305 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+2*gm(1,2)*(dgm01(2,3)& 306 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))) 307 cm(10,5)=gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 308 & +2*gm(1,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 309 cm(1,6)=2*(gm(1,2)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(dgm01(1,2)& 310 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))) 311 cm(2,6)=2*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 312 & +gm(1,2)*(dgm01(2,2)*dgm10(1,1)+8*dgm01(1,2)*dgm10(1,2)+dgm01(1,1)& 313 & *dgm10(2,2))+2*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 314 cm(3,6)=2*gm(2,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 315 & +gm(1,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 316 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+gm(1,1)*(dgm01(3,3)& 317 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))) 318 cm(4,6)=2*(gm(1,2)*dgm01(2,3)*dgm10(1,1)+4*gm(1,3)*dgm01(1,2)& 319 & *dgm10(1,2)+2*gm(1,2)*dgm01(1,3)*dgm10(1,2)+2*gm(1,1)*dgm01(2,3)& 320 & *dgm10(1,2)+gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 321 & +2*gm(1,2)*dgm01(1,2)*dgm10(1,3)+gm(2,2)*(dgm01(1,3)*dgm10(1,1)& 322 & +dgm01(1,1)*dgm10(1,3))+gm(1,2)*dgm01(1,1)*dgm10(2,3)+2*gm(1,1)& 323 & *dgm01(1,2)*dgm10(2,3)) 324 cm(5,6)=2*(gm(2,3)*dgm01(1,1)*dgm10(1,1)+gm(1,2)*dgm01(1,3)& 325 & *dgm10(1,1)+2*gm(1,1)*dgm01(1,3)*dgm10(1,2)+gm(1,3)*(dgm01(1,2)& 326 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+gm(1,2)*dgm01(1,1)*dgm10(1,3)& 327 & +2*gm(1,1)*dgm01(1,2)*dgm10(1,3)) 328 cm(6,6)=2*(gm(2,2)*dgm01(1,1)*dgm10(1,1)+2*(2*gm(1,1)*dgm01(1,2)& 329 & *dgm10(1,2)+gm(1,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2)))) 330 cm(7,6)=gm(2,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 331 & +2*gm(1,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 332 cm(8,6)=2*gm(2,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 333 & +gm(2,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 334 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+gm(1,2)*(dgm01(3,3)& 335 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))) 336 cm(9,6)=gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 337 & +2*(gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+gm(2,2)& 338 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+2*gm(1,2)*(dgm01(2,3)& 339 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))) 340 cm(10,6)=gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 341 & +2*gm(1,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 342 cm(1,7)=gm(1,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)) 343 cm(2,7)=2*(gm(1,2)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(dgm01(2,2)& 344 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 345 cm(3,7)=2*gm(2,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 346 & +gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 347 cm(4,7)=2*(gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 348 & +gm(2,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+gm(1,2)& 349 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 350 cm(5,7)=gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 351 & +2*gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)) 352 cm(6,7)=gm(2,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 353 & +2*gm(1,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)) 354 cm(7,7)=2*gm(2,2)*dgm01(2,2)*dgm10(2,2) 355 cm(8,7)=2*gm(2,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 356 & +gm(2,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 357 cm(9,7)=2*(gm(2,3)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(dgm01(2,3)& 358 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 359 cm(10,7)=gm(2,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 360 cm(1,8)=2*gm(1,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 361 & +gm(1,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 362 cm(2,8)=2*gm(1,3)*dgm01(2,3)*dgm10(2,2)+gm(1,2)*dgm01(3,3)*dgm10(2,2)& 363 & +2*gm(1,3)*dgm01(2,2)*dgm10(2,3)+4*gm(2,3)*(dgm01(2,3)*dgm10(1,2)& 364 & +dgm01(1,2)*dgm10(2,3))+gm(1,2)*dgm01(2,2)*dgm10(3,3)+2*gm(2,2)& 365 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 366 cm(3,8)=2*(gm(1,3)*dgm01(3,3)*dgm10(2,3)+2*gm(3,3)*(dgm01(2,3)& 367 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,3)*dgm01(2,3)*dgm10(3,3)& 368 & +gm(1,2)*dgm01(3,3)*dgm10(3,3)+gm(2,3)*(dgm01(3,3)*dgm10(1,3)& 369 & +dgm01(1,3)*dgm10(3,3))) 370 cm(4,8)=2*(gm(2,2)*dgm01(3,3)*dgm10(1,3)+4*gm(1,3)*dgm01(2,3)& 371 & *dgm10(2,3)+gm(1,2)*dgm01(3,3)*dgm10(2,3)+2*gm(3,3)*(dgm01(2,3)& 372 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+gm(2,2)*dgm01(1,3)*dgm10(3,3)& 373 & +gm(1,2)*dgm01(2,3)*dgm10(3,3)+gm(2,3)*(dgm01(3,3)*dgm10(1,2)& 374 & +2*dgm01(2,3)*dgm10(1,3)+2*dgm01(1,3)*dgm10(2,3)+dgm01(1,2)*dgm10(3,3))) 375 cm(5,8)=2*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 376 & +gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 377 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,2)*(dgm01(3,3)& 378 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 379 cm(6,8)=2*gm(2,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 380 & +gm(2,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+2*(2*gm(1,3)& 381 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+gm(1,2)*(dgm01(3,3)& 382 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))) 383 cm(7,8)=2*gm(2,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 384 & +gm(2,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 385 cm(8,8)=2*(4*gm(3,3)*dgm01(2,3)*dgm10(2,3)+gm(2,2)*dgm01(3,3)& 386 & *dgm10(3,3)+2*gm(2,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 387 cm(9,8)=2*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 388 & +gm(2,3)*(dgm01(3,3)*dgm10(2,2)+8*dgm01(2,3)*dgm10(2,3)+dgm01(2,2)& 389 & *dgm10(3,3))+2*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 390 cm(10,8)=2*(gm(2,3)*dgm01(3,3)*dgm10(3,3)+gm(3,3)*(dgm01(3,3)& 391 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 392 cm(1,9)=gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 393 & +2*gm(1,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)) 394 cm(2,9)=2*(gm(1,3)*dgm01(2,2)*dgm10(2,2)+gm(1,2)*dgm01(2,3)& 395 & *dgm10(2,2)+gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 396 & +gm(1,2)*dgm01(2,2)*dgm10(2,3)+2*gm(2,2)*(dgm01(2,3)*dgm10(1,2)& 397 & +dgm01(1,2)*dgm10(2,3))) 398 cm(3,9)=gm(1,3)*dgm01(3,3)*dgm10(2,2)+2*gm(3,3)*(dgm01(2,2)& 399 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+2*gm(1,2)*dgm01(3,3)*dgm10(2,3)& 400 & +4*gm(2,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(1,3)& 401 & *dgm01(2,2)*dgm10(3,3)+2*gm(1,2)*dgm01(2,3)*dgm10(3,3) 402 cm(4,9)=2*(2*gm(2,2)*dgm01(2,3)*dgm10(1,3)+gm(1,3)*dgm01(2,3)& 403 & *dgm10(2,2)+gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 404 & +2*gm(2,2)*dgm01(1,3)*dgm10(2,3)+gm(1,3)*dgm01(2,2)*dgm10(2,3)& 405 & +4*gm(1,2)*dgm01(2,3)*dgm10(2,3)+gm(2,3)*(2*dgm01(2,3)*dgm10(1,2)& 406 & +dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)+2*dgm01(1,2)*dgm10(2,3))) 407 cm(5,9)=gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 408 & +2*(gm(1,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+gm(2,3)& 409 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+2*gm(1,2)*(dgm01(2,3)& 410 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))) 411 cm(6,9)=gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 412 & +2*(gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+gm(2,2)& 413 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+2*gm(1,2)*(dgm01(2,3)& 414 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))) 415 cm(7,9)=2*(gm(2,3)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(dgm01(2,3)& 416 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 417 cm(8,9)=2*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 418 & +gm(2,3)*(dgm01(3,3)*dgm10(2,2)+8*dgm01(2,3)*dgm10(2,3)+dgm01(2,2)& 419 & *dgm10(3,3))+2*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 420 cm(9,9)=2*(gm(3,3)*dgm01(2,2)*dgm10(2,2)+2*(2*gm(2,2)*dgm01(2,3)& 421 & *dgm10(2,3)+gm(2,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))) 422 cm(10,9)=gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))& 423 & +2*gm(2,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 424 cm(1,10)=gm(1,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)) 425 cm(2,10)=2*gm(2,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 426 & +gm(1,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 427 cm(3,10)=2*(gm(1,3)*dgm01(3,3)*dgm10(3,3)+gm(3,3)*(dgm01(3,3)& 428 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 429 cm(4,10)=2*(gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 430 & +gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+gm(1,3)& 431 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 432 cm(5,10)=gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 433 & +2*gm(1,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)) 434 cm(6,10)=gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 435 & +2*gm(1,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)) 436 cm(7,10)=gm(2,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)) 437 cm(8,10)=2*(gm(2,3)*dgm01(3,3)*dgm10(3,3)+gm(3,3)*(dgm01(3,3)& 438 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 439 cm(9,10)=gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))& 440 & +2*gm(2,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)) 441 cm(10,10)=2*gm(3,3)*dgm01(3,3)*dgm10(3,3) 442 elseif(rank==2)then 443 cm(1,1)=2*gm(1,1)**2*dgm01(1,1)*dgm10(1,1) 444 cm(2,1)=(48*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(1,1)*dgm10(1,1)& 445 & +384*gm(1,1)*gm(1,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 446 & +96*gm(1,1)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)))& 447 & /96.d0 448 cm(3,1)=(48*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 449 & +384*gm(1,1)*gm(1,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 450 & +96*gm(1,1)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))& 451 & /96.d0 452 cm(4,1)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,1)& 453 & *dgm10(1,1)+192*gm(1,1)*gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 454 & *dgm10(1,2))+192*gm(1,1)*gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 455 & *dgm10(1,3))+96*gm(1,1)**2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)& 456 & *dgm10(2,3)))/48.d0 457 cm(5,1)=gm(1,1)*(4*gm(1,3)*dgm01(1,1)*dgm10(1,1)+2*gm(1,1)*(dgm01(1,3)& 458 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))) 459 cm(6,1)=gm(1,1)*(4*gm(1,2)*dgm01(1,1)*dgm10(1,1)+2*gm(1,1)*(dgm01(1,2)& 460 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))) 461 cm(7,1)=(72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(1,2)*dgm10(1,1)& 462 & +dgm01(1,1)*dgm10(1,2))+288*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,1)& 463 & +dgm01(1,1)*dgm10(2,2)))/144.d0 464 cm(8,1)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 465 & +dgm01(1,1)*dgm10(1,2))+48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 466 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)*gm(1,3)& 467 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,2)& 468 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/48.d0 469 cm(9,1)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,2)& 470 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 471 & *gm(2,2))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)& 472 & *gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+192*gm(1,1)& 473 & *gm(1,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)))/48.d0 474 cm(10,1)=(72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 475 & +dgm01(1,1)*dgm10(1,3))+288*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,1)& 476 & +dgm01(1,1)*dgm10(3,3)))/144.d0 477 cm(11,1)=((6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,1)& 478 & +dgm01(1,1)*dgm10(2,2)))/4.d0 479 cm(12,1)=((6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 480 & +dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 481 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+(6*gm(1,2)**2-2*gm(1,1)& 482 & *gm(2,2))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/4.d0 483 cm(13,1)=(4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)& 484 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 485 & *gm(2,2))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)))/8.d0 486 cm(14,1)=(4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 487 & +dgm01(1,1)*dgm10(2,3))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 488 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/8.d0 489 cm(15,1)=((6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 490 & +dgm01(1,1)*dgm10(3,3)))/4.d0 491 cm(1,2)=(48*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(1,1)*dgm10(1,1)& 492 & +384*gm(1,1)*gm(1,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 493 & +96*gm(1,1)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2)))& 494 & /96.d0 495 cm(2,2)=(32*gm(2,2)**2*dgm01(1,1)*dgm10(1,1)+128*(1*gm(1,2)& 496 & **2+3*gm(1,1)*gm(2,2))*dgm01(1,2)*dgm10(1,2)+128*gm(1,2)*gm(2,2)& 497 & *(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+32*gm(1,1)**2*dgm01(2,2)& 498 & *dgm10(2,2)+8*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,1)& 499 & +dgm01(1,1)*dgm10(2,2))+128*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,2)& 500 & +dgm01(1,2)*dgm10(2,2)))/16.d0 501 cm(3,2)=(8*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 502 & +16*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 503 & +dgm01(1,1)*dgm10(1,2))+16*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 504 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+64*(1*gm(1,2)& 505 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 506 & *dgm10(1,3))+4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 507 & +dgm01(1,1)*dgm10(2,2))+64*gm(1,1)*gm(1,3)*(dgm01(2,2)*dgm10(1,3)& 508 & +dgm01(1,3)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)& 509 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+64*gm(1,1)*gm(1,2)*(dgm01(3,3)& 510 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+16*gm(1,1)**2*(dgm01(3,3)& 511 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/16.d0 512 cm(4,2)=(32*gm(2,2)*gm(2,3)*dgm01(1,1)*dgm10(1,1)+64*(1*gm(1,2)& 513 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,2)*dgm10(1,2)+16*(3*gm(1,3)& 514 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 515 & +8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)& 516 & +dgm01(1,1)*dgm10(1,2))+32*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,1)& 517 & +dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(1,3)& 518 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 519 & *gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+32*gm(1,1)& 520 & *gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+32*gm(1,1)& 521 & *gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+4*(6*gm(1,2)& 522 & **2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 523 & +64*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 524 & +16*gm(1,1)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))& 525 & /8.d0 526 cm(5,2)=(24*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,1)& 527 & *dgm10(1,1)+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,2)& 528 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+12*(6*gm(1,2)**2-2*gm(1,1)& 529 & *gm(2,2))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)& 530 & *gm(1,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)& 531 & *gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+48*gm(1,1)& 532 & **2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)))/24.d0 533 cm(6,2)=5*gm(1,2)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 534 & +gm(1,2)*(4*gm(2,2)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(2*dgm01(2,2)& 535 & *dgm10(1,1)+16*dgm01(1,2)*dgm10(1,2)+2*dgm01(1,1)*dgm10(2,2)))& 536 & +gm(1,1)*(5*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 537 & +2*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 538 cm(7,2)=2*gm(2,2)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 539 & +gm(2,2)*(gm(1,2)*(2*dgm01(2,2)*dgm10(1,1)+16*dgm01(1,2)*dgm10(1,2)& 540 & +2*dgm01(1,1)*dgm10(2,2))+5*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 541 & *dgm10(2,2)))+gm(1,2)*(4*gm(1,1)*dgm01(2,2)*dgm10(2,2)+5*gm(1,2)& 542 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 543 cm(8,2)=(32*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,2)& 544 & *dgm10(1,2)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 545 & +dgm01(1,1)*dgm10(1,2))+32*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,1)& 546 & +dgm01(1,1)*dgm10(1,3))+32*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 547 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 548 & *gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+8*(6*gm(1,2)& 549 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 550 & *dgm10(2,2))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)& 551 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 552 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+32*gm(1,1)& 553 & *gm(1,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*gm(1,2)& 554 & *gm(2,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+16*(1*gm(1,2)& 555 & **2+3*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 556 & +16*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 557 & /8.d0 558 cm(9,2)=(64*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,2)*dgm10(1,2)& 559 & +32*gm(2,2)*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 560 & +16*gm(2,2)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 561 & +64*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 562 & +32*gm(1,1)*gm(1,3)*dgm01(2,2)*dgm10(2,2)+4*(-2*gm(1,3)*gm(2,2)& 563 & +6*gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 564 & +8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)& 565 & +dgm01(1,2)*dgm10(2,2))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 566 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 567 & *gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+32*gm(1,2)& 568 & *gm(2,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)& 569 & **2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 570 & +32*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))& 571 & /8.d0 572 cm(10,2)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 573 & +dgm01(1,1)*dgm10(1,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 574 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(6*gm(1,3)& 575 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 576 & +12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 577 & +dgm01(1,1)*dgm10(3,3))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 578 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,1)*gm(1,3)& 579 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 580 cm(11,2)=(48*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(2,2)*dgm10(2,2)& 581 & +96*gm(2,2)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 582 & +384*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)))& 583 & /96.d0 584 cm(12,2)=(8*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(2,2)*dgm10(2,2)& 585 & +4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 586 & *dgm10(2,2))+16*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 587 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+64*gm(2,2)*gm(2,3)*(dgm01(2,3)& 588 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+64*(3*gm(1,3)*gm(2,2)+gm(1,2)& 589 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(6*gm(1,2)& 590 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 591 & *dgm10(2,3))+16*gm(2,2)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 592 & *dgm10(3,3))+64*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 593 & *dgm10(3,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(2,2)& 594 & +dgm01(2,2)*dgm10(3,3)))/16.d0 595 cm(13,2)=(24*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(2,2)& 596 & *dgm10(2,2)+48*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 597 & *dgm10(2,2))+48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 598 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*gm(2,2)**2*(dgm01(2,3)& 599 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+192*gm(1,2)*gm(2,2)*(dgm01(2,3)& 600 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 601 & *gm(2,2))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/24.d0 602 cm(14,2)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 603 & +dgm01(1,1)*dgm10(2,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 604 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(6*gm(1,3)& 605 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 606 & +48*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 607 & +48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 608 & +dgm01(1,2)*dgm10(3,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 609 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 610 cm(15,2)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 611 & +dgm01(1,1)*dgm10(3,3))+4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 612 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+(6*gm(1,3)**2-2*gm(1,1)& 613 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/4.d0 614 cm(1,3)=(48*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 615 & +384*gm(1,1)*gm(1,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 616 & +96*gm(1,1)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))& 617 & /96.d0 618 cm(2,3)=(8*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 619 & +16*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 620 & +dgm01(1,1)*dgm10(1,2))+16*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 621 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+64*(1*gm(1,2)& 622 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 623 & *dgm10(1,3))+4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 624 & +dgm01(1,1)*dgm10(2,2))+64*gm(1,1)*gm(1,3)*(dgm01(2,2)*dgm10(1,3)& 625 & +dgm01(1,3)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)& 626 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+64*gm(1,1)*gm(1,2)*(dgm01(3,3)& 627 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+16*gm(1,1)**2*(dgm01(3,3)& 628 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/16.d0 629 cm(3,3)=(32*gm(3,3)**2*dgm01(1,1)*dgm10(1,1)+128*(1*gm(1,3)& 630 & **2+3*gm(1,1)*gm(3,3))*dgm01(1,3)*dgm10(1,3)+128*gm(1,3)*gm(3,3)& 631 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+32*gm(1,1)**2*dgm01(3,3)& 632 & *dgm10(3,3)+8*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 633 & +dgm01(1,1)*dgm10(3,3))+128*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,3)& 634 & +dgm01(1,3)*dgm10(3,3)))/16.d0 635 cm(4,3)=(32*gm(2,3)*gm(3,3)*dgm01(1,1)*dgm10(1,1)+32*gm(1,3)& 636 & *gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+64*(1*gm(1,2)& 637 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,3)*dgm10(1,3)+8*(6*gm(1,3)& 638 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 639 & *dgm10(1,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,3)& 640 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,3)**2+3*gm(1,1)& 641 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,3)& 642 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 643 & +64*gm(1,1)*gm(1,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 644 & +4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 645 & +dgm01(1,1)*dgm10(3,3))+32*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,2)& 646 & +dgm01(1,2)*dgm10(3,3))+32*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(1,3)& 647 & +dgm01(1,3)*dgm10(3,3))+16*gm(1,1)**2*(dgm01(3,3)*dgm10(2,3)& 648 & +dgm01(2,3)*dgm10(3,3)))/8.d0 649 cm(5,3)=5*gm(1,3)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 650 & +gm(1,3)*(4*gm(3,3)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(2*dgm01(3,3)& 651 & *dgm10(1,1)+16*dgm01(1,3)*dgm10(1,3)+2*dgm01(1,1)*dgm10(3,3)))& 652 & +gm(1,1)*(5*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 653 & +2*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 654 cm(6,3)=(24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,1)& 655 & *dgm10(1,1)+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 656 & +dgm01(1,1)*dgm10(1,2))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 657 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)*gm(1,3)& 658 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,2)& 659 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+48*gm(1,1)**2*(dgm01(3,3)& 660 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/24.d0 661 cm(7,3)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 662 & +dgm01(1,1)*dgm10(1,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 663 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(6*gm(1,3)& 664 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 665 & *dgm10(2,2))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)& 666 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+12*(6*gm(1,2)**2-2*gm(1,1)& 667 & *gm(2,2))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,1)& 668 & *gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 669 cm(8,3)=(16*gm(3,3)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 670 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 671 & +32*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 672 & +64*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 673 & +32*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 674 & +32*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 675 & *dgm10(2,3))+32*gm(1,1)*gm(1,2)*dgm01(3,3)*dgm10(3,3)+4*(6*gm(1,3)& 676 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 677 & *dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)& 678 & +dgm01(1,2)*dgm10(3,3))+8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 679 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+16*(1*gm(1,2)& 680 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 681 & *dgm10(3,3))+32*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 682 & *dgm10(3,3)))/8.d0 683 cm(9,3)=(32*gm(2,3)*gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 684 & *dgm10(1,2))+32*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,3)& 685 & *dgm10(1,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 686 & +dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 687 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+16*gm(1,3)*gm(3,3)& 688 & *(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+16*(1*gm(1,3)& 689 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 690 & +8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 691 & +dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 692 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(6*gm(1,2)*gm(1,3)& 693 & -2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 694 & +4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 695 & *dgm10(3,3))+16*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 696 & *dgm10(3,3))+32*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 697 & *dgm10(3,3)))/8.d0 698 cm(10,3)=2*gm(3,3)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 699 & +gm(3,3)*(gm(1,3)*(2*dgm01(3,3)*dgm10(1,1)+16*dgm01(1,3)*dgm10(1,3)& 700 & +2*dgm01(1,1)*dgm10(3,3))+5*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 701 & *dgm10(3,3)))+gm(1,3)*(4*gm(1,1)*dgm01(3,3)*dgm10(3,3)+5*gm(1,3)& 702 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 703 cm(11,3)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 704 & +dgm01(1,1)*dgm10(2,2))+4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 705 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(6*gm(1,2)**2-2*gm(1,1)& 706 & *gm(2,2))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/4.d0 707 cm(12,3)=(16*gm(3,3)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 708 & +64*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 709 & +64*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 710 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 711 & +dgm01(1,3)*dgm10(2,3))+8*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(3,3)& 712 & *dgm10(3,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 713 & +dgm01(1,1)*dgm10(3,3))+16*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 714 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 715 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+16*(6*gm(1,2)& 716 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 717 & *dgm10(3,3)))/16.d0 718 cm(13,3)=(48*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 719 & *dgm10(2,2))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)& 720 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+12*(6*gm(2,3)**2-2*gm(2,2)& 721 & *gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*(-2*gm(1,3)& 722 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 723 & *dgm10(2,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)& 724 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 725 & *gm(2,2))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 726 cm(14,3)=(48*gm(3,3)**2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 727 & +192*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 728 & +24*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(3,3)*dgm10(3,3)& 729 & +48*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 730 & +48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,3)& 731 & +dgm01(1,3)*dgm10(3,3))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)& 732 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 733 cm(15,3)=(48*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(3,3)*dgm10(3,3)& 734 & +96*gm(3,3)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 735 & +384*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))& 736 & /96.d0 737 cm(1,4)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,1)& 738 & *dgm10(1,1)+192*gm(1,1)*gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 739 & *dgm10(1,2))+192*gm(1,1)*gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 740 & *dgm10(1,3))+96*gm(1,1)**2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)& 741 & *dgm10(2,3)))/48.d0 742 cm(2,4)=(32*gm(2,2)*gm(2,3)*dgm01(1,1)*dgm10(1,1)+64*(1*gm(1,2)& 743 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,2)*dgm10(1,2)+16*(3*gm(1,3)& 744 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 745 & +8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)& 746 & +dgm01(1,1)*dgm10(1,2))+32*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,1)& 747 & +dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(1,3)& 748 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 749 & *gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+32*gm(1,1)& 750 & *gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+32*gm(1,1)& 751 & *gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+4*(6*gm(1,2)& 752 & **2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 753 & +64*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 754 & +16*gm(1,1)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))& 755 & /8.d0 756 cm(3,4)=(32*gm(2,3)*gm(3,3)*dgm01(1,1)*dgm10(1,1)+32*gm(1,3)& 757 & *gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+64*(1*gm(1,2)& 758 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,3)*dgm10(1,3)+8*(6*gm(1,3)& 759 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 760 & *dgm10(1,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,3)& 761 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,3)**2+3*gm(1,1)& 762 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,3)& 763 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 764 & +64*gm(1,1)*gm(1,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 765 & +4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 766 & +dgm01(1,1)*dgm10(3,3))+32*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,2)& 767 & +dgm01(1,2)*dgm10(3,3))+32*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(1,3)& 768 & +dgm01(1,3)*dgm10(3,3))+16*gm(1,1)**2*(dgm01(3,3)*dgm10(2,3)& 769 & +dgm01(2,3)*dgm10(3,3)))/8.d0 770 cm(4,4)=(8*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 771 & +32*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*dgm01(1,2)*dgm10(1,2)+16*(1*gm(1,3)& 772 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 773 & *dgm10(1,2))+32*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*dgm01(1,3)*dgm10(1,3)& 774 & +16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,3)*dgm10(1,1)& 775 & +dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 776 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+32*gm(1,1)**2*dgm01(2,3)& 777 & *dgm10(2,3)+8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)& 778 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+64*gm(1,1)*gm(1,3)*(dgm01(2,3)& 779 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+64*gm(1,1)*gm(1,2)*(dgm01(2,3)& 780 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3)))/4.d0 781 cm(5,4)=(24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,1)& 782 & *dgm10(1,1)+24*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 783 & +dgm01(1,1)*dgm10(1,2))+192*gm(1,1)*gm(1,2)*dgm01(1,3)*dgm10(1,3)& 784 & +12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,1)& 785 & +dgm01(1,1)*dgm10(1,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 786 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)*gm(1,3)& 787 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,3)& 788 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*gm(1,1)**2*(dgm01(2,3)& 789 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3)))/12.d0 790 cm(6,4)=(24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,1)*dgm10(1,1)& 791 & +192*gm(1,1)*gm(1,3)*dgm01(1,2)*dgm10(1,2)+12*(6*gm(1,2)*gm(1,3)& 792 & -2*gm(1,1)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 793 & +24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)& 794 & +dgm01(1,1)*dgm10(1,2))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(1,3)& 795 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)*gm(1,2)*(dgm01(1,3)& 796 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,2)*(dgm01(2,3)& 797 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*gm(1,1)**2*(dgm01(2,3)& 798 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))/12.d0 799 cm(7,4)=(48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,2)& 800 & *dgm10(1,2)+48*gm(2,2)*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 801 & *dgm10(1,2))+96*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 802 & *dgm10(1,3))+12*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 803 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 804 & *gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+24*(1*gm(1,2)& 805 & **2+3*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 806 & +12*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 807 & *dgm10(2,3))+48*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 808 & *dgm10(2,3)))/12.d0 809 cm(8,4)=(64*gm(1,3)*gm(3,3)*dgm01(1,2)*dgm10(1,2)+16*gm(2,3)& 810 & *gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+32*(3*gm(1,3)& 811 & *gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,3)*dgm10(1,3)+8*(1*gm(2,3)& 812 & **2+3*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 813 & +8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)& 814 & +dgm01(1,2)*dgm10(1,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 815 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+64*gm(1,1)*gm(1,3)& 816 & *dgm01(2,3)*dgm10(2,3)+8*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 817 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 818 & *gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,3)& 819 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 820 & +8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)& 821 & +dgm01(1,3)*dgm10(2,3))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 822 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+4*(3*gm(1,3)*gm(2,2)& 823 & +gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 824 & +8*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 825 & +dgm01(1,2)*dgm10(3,3))+8*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(3,3)& 826 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+16*gm(1,1)*gm(1,2)*(dgm01(3,3)& 827 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/4.d0 828 cm(9,4)=(32*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,2)& 829 & *dgm10(1,2)+8*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 830 & +dgm01(1,1)*dgm10(1,2))+64*gm(1,2)*gm(2,2)*dgm01(1,3)*dgm10(1,3)& 831 & +16*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 832 & +16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)& 833 & +dgm01(1,2)*dgm10(1,3))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 834 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(1*gm(1,3)*gm(2,3)& 835 & +3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 836 & +8*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 837 & *dgm10(2,2))+8*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)& 838 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+64*gm(1,1)*gm(1,2)*dgm01(2,3)& 839 & *dgm10(2,3)+8*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 840 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 841 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,2)& 842 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 843 & *dgm10(2,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,3)& 844 & +dgm01(1,3)*dgm10(2,3))+16*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 845 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+16*gm(1,1)*gm(1,3)*(dgm01(2,3)& 846 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/4.d0 847 cm(10,4)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,3)& 848 & *dgm10(1,3)+48*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 849 & *dgm10(1,3))+96*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 850 & *dgm10(1,3))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 851 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+12*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 852 & *gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+24*(1*gm(1,3)& 853 & **2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 854 & +24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)& 855 & +dgm01(1,3)*dgm10(3,3))+48*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,3)& 856 & +dgm01(2,3)*dgm10(3,3)))/12.d0 857 cm(11,4)=(96*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 858 & *dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)& 859 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+192*gm(1,2)*gm(2,2)*(dgm01(2,2)& 860 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 861 & *gm(2,2))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/48.d0 862 cm(12,4)=(16*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 863 & *dgm10(2,2))+32*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 864 & *dgm10(2,2))+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 865 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+32*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 866 & *gm(2,3))*dgm01(2,3)*dgm10(2,3)+16*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))& 867 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,3)& 868 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 869 & *dgm10(2,3))+32*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 870 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 871 & *gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*gm(2,2)& 872 & *gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+8*(-2*gm(1,3)& 873 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 874 & *dgm10(3,3))+32*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 875 & *dgm10(3,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(2,3)& 876 & +dgm01(2,3)*dgm10(3,3)))/8.d0 877 cm(13,4)=(12*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 878 & +dgm01(1,1)*dgm10(2,2))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 879 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+24*(3*gm(1,3)& 880 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 881 & +24*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(2,3)*dgm10(2,3)+48*gm(2,2)& 882 & *gm(2,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+24*(-2*gm(1,3)& 883 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 884 & *dgm10(2,3))+96*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 885 & *dgm10(2,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)& 886 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/12.d0 887 cm(14,4)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 888 & +48*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 889 & +96*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 890 & +24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 891 & +dgm01(1,3)*dgm10(2,3))+12*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(3,3)& 892 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 893 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+24*(3*gm(1,3)& 894 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 895 & +12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)& 896 & +dgm01(2,3)*dgm10(3,3)))/12.d0 897 cm(15,4)=(96*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 898 & *dgm10(3,3))+192*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 899 & *dgm10(3,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)& 900 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+24*(6*gm(1,3)**2-2*gm(1,1)& 901 & *gm(3,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/48.d0 902 cm(1,5)=gm(1,1)*(4*gm(1,3)*dgm01(1,1)*dgm10(1,1)+2*gm(1,1)*(dgm01(1,3)& 903 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))) 904 cm(2,5)=(24*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,1)& 905 & *dgm10(1,1)+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,2)& 906 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+12*(6*gm(1,2)**2-2*gm(1,1)& 907 & *gm(2,2))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)& 908 & *gm(1,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)& 909 & *gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+48*gm(1,1)& 910 & **2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)))/24.d0 911 cm(3,5)=5*gm(1,3)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 912 & +gm(1,3)*(4*gm(3,3)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(2*dgm01(3,3)& 913 & *dgm10(1,1)+16*dgm01(1,3)*dgm10(1,3)+2*dgm01(1,1)*dgm10(3,3)))& 914 & +gm(1,1)*(5*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 915 & +2*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 916 cm(4,5)=(24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,1)& 917 & *dgm10(1,1)+24*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 918 & +dgm01(1,1)*dgm10(1,2))+192*gm(1,1)*gm(1,2)*dgm01(1,3)*dgm10(1,3)& 919 & +12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,1)& 920 & +dgm01(1,1)*dgm10(1,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 921 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)*gm(1,3)& 922 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,3)& 923 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*gm(1,1)**2*(dgm01(2,3)& 924 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3)))/12.d0 925 cm(5,5)=(72*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*dgm01(1,1)*dgm10(1,1)& 926 & +288*gm(1,1)**2*dgm01(1,3)*dgm10(1,3)+288*gm(1,1)*gm(1,3)*(dgm01(1,3)& 927 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3)))/36.d0 928 cm(6,5)=(72*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,1)& 929 & *dgm10(1,1)+144*gm(1,1)*gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 930 & *dgm10(1,2))+144*gm(1,1)*gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 931 & *dgm10(1,3))+144*gm(1,1)**2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 932 & *dgm10(1,3)))/36.d0 933 cm(7,5)=(36*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(1,2)& 934 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+36*(6*gm(1,2)**2-2*gm(1,1)& 935 & *gm(2,2))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+36*(1*gm(1,2)& 936 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 937 & *dgm10(2,2))+144*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 938 & *dgm10(2,2)))/36.d0 939 cm(8,5)=(48*gm(1,3)*gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 940 & *dgm10(1,2))+48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,3)& 941 & *dgm10(1,3)+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,3)& 942 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+12*(6*gm(1,3)**2-2*gm(1,1)& 943 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+24*(1*gm(1,3)& 944 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 945 & +96*gm(1,1)*gm(1,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 946 & +12*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 947 & +dgm01(1,1)*dgm10(3,3))+48*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(1,3)& 948 & +dgm01(1,3)*dgm10(3,3)))/12.d0 949 cm(9,5)=(24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,2)& 950 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 951 & *gm(2,2))*dgm01(1,3)*dgm10(1,3)+12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 952 & *gm(2,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+24*(6*gm(1,2)& 953 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 954 & *dgm10(1,3))+12*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)& 955 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+48*gm(1,1)*gm(1,3)*(dgm01(2,2)& 956 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 957 & *gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)& 958 & *gm(1,2)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3)))/12.d0 959 cm(10,5)=(72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 960 & +144*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 961 & +36*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 962 & *dgm10(3,3))+144*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 963 & *dgm10(3,3)))/36.d0 964 cm(11,5)=(4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)& 965 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 966 & *gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)))/8.d0 967 cm(12,5)=(48*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 968 & *dgm10(2,2))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)& 969 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 970 & *gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*(6*gm(1,2)& 971 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 972 & *dgm10(2,3))+12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)& 973 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 974 & *gm(2,2))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))/24.d0 975 cm(13,5)=(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 976 & +dgm01(1,1)*dgm10(2,2))+(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 977 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(-2*gm(1,3)*gm(2,2)& 978 & +6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 979 & +(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 980 & *dgm10(2,3)) 981 cm(14,5)=(144*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)& 982 & *dgm10(2,3))+36*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 983 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+36*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 984 & *gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+36*(6*gm(1,2)& 985 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 986 & *dgm10(3,3)))/36.d0 987 cm(15,5)=(288*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 988 & *dgm10(3,3))+72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)& 989 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))/144.d0 990 cm(1,6)=gm(1,1)*(4*gm(1,2)*dgm01(1,1)*dgm10(1,1)+2*gm(1,1)*(dgm01(1,2)& 991 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))) 992 cm(2,6)=5*gm(1,2)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 993 & +gm(1,2)*(4*gm(2,2)*dgm01(1,1)*dgm10(1,1)+gm(1,1)*(2*dgm01(2,2)& 994 & *dgm10(1,1)+16*dgm01(1,2)*dgm10(1,2)+2*dgm01(1,1)*dgm10(2,2)))& 995 & +gm(1,1)*(5*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 996 & +2*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 997 cm(3,6)=(24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,1)& 998 & *dgm10(1,1)+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 999 & +dgm01(1,1)*dgm10(1,2))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1000 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)*gm(1,3)& 1001 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,2)& 1002 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+48*gm(1,1)**2*(dgm01(3,3)& 1003 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/24.d0 1004 cm(4,6)=(24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,1)*dgm10(1,1)& 1005 & +192*gm(1,1)*gm(1,3)*dgm01(1,2)*dgm10(1,2)+12*(6*gm(1,2)*gm(1,3)& 1006 & -2*gm(1,1)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1007 & +24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(1,2)*dgm10(1,1)& 1008 & +dgm01(1,1)*dgm10(1,2))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(1,3)& 1009 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)*gm(1,2)*(dgm01(1,3)& 1010 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+48*gm(1,1)*gm(1,2)*(dgm01(2,3)& 1011 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*gm(1,1)**2*(dgm01(2,3)& 1012 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))/12.d0 1013 cm(5,6)=(72*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(1,1)& 1014 & *dgm10(1,1)+144*gm(1,1)*gm(1,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 1015 & *dgm10(1,2))+144*gm(1,1)*gm(1,2)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 1016 & *dgm10(1,3))+144*gm(1,1)**2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1017 & *dgm10(1,3)))/36.d0 1018 cm(6,6)=(72*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*dgm01(1,1)*dgm10(1,1)& 1019 & +288*gm(1,1)**2*dgm01(1,2)*dgm10(1,2)+288*gm(1,1)*gm(1,2)*(dgm01(1,2)& 1020 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2)))/36.d0 1021 cm(7,6)=(72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(1,2)*dgm10(1,2)& 1022 & +144*gm(1,2)*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1023 & +36*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1024 & *dgm10(2,2))+144*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1025 & *dgm10(2,2)))/36.d0 1026 cm(8,6)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,2)*dgm10(1,2)& 1027 & +12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1028 & +dgm01(1,1)*dgm10(1,2))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1029 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+24*(6*gm(1,2)& 1030 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1031 & *dgm10(1,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)& 1032 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,3)*(dgm01(2,3)& 1033 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(1*gm(1,2)**2+3*gm(1,1)& 1034 & *gm(2,2))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+48*gm(1,1)& 1035 & *gm(1,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/12.d0 1036 cm(9,6)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,2)& 1037 & *dgm10(1,2)+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,2)& 1038 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+48*gm(1,2)*gm(2,2)*(dgm01(1,3)& 1039 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 1040 & *gm(2,2))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(1*gm(1,2)& 1041 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1042 & *dgm10(2,2))+48*gm(1,1)*gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1043 & *dgm10(2,2))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 1044 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,2)*(dgm01(2,3)& 1045 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))/12.d0 1046 cm(10,6)=(36*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)& 1047 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+36*(6*gm(1,3)**2-2*gm(1,1)& 1048 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+36*(1*gm(1,2)& 1049 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1050 & *dgm10(3,3))+144*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1051 & *dgm10(3,3)))/36.d0 1052 cm(11,6)=(288*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1053 & *dgm10(2,2))+72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,2)& 1054 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2)))/144.d0 1055 cm(12,6)=(12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 1056 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+12*(6*gm(1,3)**2-2*gm(1,1)& 1057 & *gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*(3*gm(1,3)& 1058 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1059 & +48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)& 1060 & +dgm01(1,2)*dgm10(2,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,1)& 1061 & +dgm01(1,1)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)& 1062 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/24.d0 1063 cm(13,6)=(36*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 1064 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+36*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 1065 & *gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+144*gm(1,2)& 1066 & *gm(2,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+36*(6*gm(1,2)& 1067 & **2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))& 1068 & /36.d0 1069 cm(14,6)=(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 1070 & +dgm01(1,1)*dgm10(2,3))+(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 1071 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+(3*gm(1,3)*gm(2,2)+gm(1,2)& 1072 & *gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+(6*gm(1,2)& 1073 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1074 & *dgm10(3,3)) 1075 cm(15,6)=(4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)& 1076 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1077 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/8.d0 1078 cm(1,7)=(72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(1,2)*dgm10(1,1)& 1079 & +dgm01(1,1)*dgm10(1,2))+288*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,1)& 1080 & +dgm01(1,1)*dgm10(2,2)))/144.d0 1081 cm(2,7)=2*gm(2,2)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1082 & +gm(2,2)*(gm(1,2)*(2*dgm01(2,2)*dgm10(1,1)+16*dgm01(1,2)*dgm10(1,2)& 1083 & +2*dgm01(1,1)*dgm10(2,2))+5*gm(1,1)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1084 & *dgm10(2,2)))+gm(1,2)*(4*gm(1,1)*dgm01(2,2)*dgm10(2,2)+5*gm(1,2)& 1085 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 1086 cm(3,7)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1087 & +dgm01(1,1)*dgm10(1,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1088 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(6*gm(1,3)& 1089 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1090 & *dgm10(2,2))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)& 1091 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+12*(6*gm(1,2)**2-2*gm(1,1)& 1092 & *gm(2,2))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,1)& 1093 & *gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 1094 cm(4,7)=(48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,2)& 1095 & *dgm10(1,2)+48*gm(2,2)*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 1096 & *dgm10(1,2))+96*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1097 & *dgm10(1,3))+12*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 1098 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 1099 & *gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+24*(1*gm(1,2)& 1100 & **2+3*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1101 & +12*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1102 & *dgm10(2,3))+48*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1103 & *dgm10(2,3)))/12.d0 1104 cm(5,7)=(36*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(1,2)& 1105 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+36*(6*gm(1,2)**2-2*gm(1,1)& 1106 & *gm(2,2))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+36*(1*gm(1,2)& 1107 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1108 & *dgm10(2,2))+144*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1109 & *dgm10(2,2)))/36.d0 1110 cm(6,7)=(72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(1,2)*dgm10(1,2)& 1111 & +144*gm(1,2)*gm(2,2)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1112 & +36*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1113 & *dgm10(2,2))+144*gm(1,1)*gm(1,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1114 & *dgm10(2,2)))/36.d0 1115 cm(7,7)=(288*gm(2,2)**2*dgm01(1,2)*dgm10(1,2)+72*(1*gm(1,2)& 1116 & **2+3*gm(1,1)*gm(2,2))*dgm01(2,2)*dgm10(2,2)+288*gm(1,2)*gm(2,2)& 1117 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)))/36.d0 1118 cm(8,7)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,2)*dgm10(1,2)& 1119 & +96*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1120 & +12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1121 & +dgm01(1,2)*dgm10(2,2))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1122 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(-2*gm(1,3)& 1123 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1124 & *dgm10(2,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)& 1125 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1126 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+12*(1*gm(1,2)**2+3*gm(1,1)& 1127 & *gm(2,2))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/12.d0 1128 cm(9,7)=(192*gm(2,2)*gm(2,3)*dgm01(1,2)*dgm10(1,2)+48*gm(2,2)& 1129 & **2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+24*(1*gm(1,2)& 1130 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(2,2)*dgm10(2,2)+24*(3*gm(1,3)& 1131 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1132 & +12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)& 1133 & +dgm01(1,2)*dgm10(2,2))+48*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,3)& 1134 & +dgm01(1,3)*dgm10(2,2))+96*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(1,2)& 1135 & +dgm01(1,2)*dgm10(2,3))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 1136 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/12.d0 1137 cm(10,7)=(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)& 1138 & +dgm01(1,2)*dgm10(1,3))+(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1139 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(-2*gm(1,3)*gm(2,2)& 1140 & +6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1141 & +(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)& 1142 & +dgm01(2,2)*dgm10(3,3)) 1143 cm(11,7)=gm(2,2)*(4*gm(1,2)*dgm01(2,2)*dgm10(2,2)+2*gm(2,2)& 1144 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 1145 cm(12,7)=(24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(2,2)& 1146 & *dgm10(2,2)+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1147 & +dgm01(1,2)*dgm10(2,2))+192*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(1,2)& 1148 & +dgm01(1,2)*dgm10(2,3))+48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1149 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+48*gm(2,2)**2*(dgm01(3,3)& 1150 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1151 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 1152 cm(13,7)=(72*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(2,2)& 1153 & *dgm10(2,2)+144*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1154 & *dgm10(2,2))+144*gm(2,2)**2*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1155 & *dgm10(2,3))+144*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1156 & *dgm10(2,3)))/36.d0 1157 cm(14,7)=(36*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)& 1158 & +dgm01(1,2)*dgm10(2,3))+36*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1159 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+144*gm(2,2)*gm(2,3)& 1160 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+36*(3*gm(1,3)& 1161 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 1162 & /36.d0 1163 cm(15,7)=(4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)& 1164 & +dgm01(1,2)*dgm10(3,3))+4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1165 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/8.d0 1166 cm(1,8)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1167 & +dgm01(1,1)*dgm10(1,2))+48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1168 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+192*gm(1,1)*gm(1,3)& 1169 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,2)& 1170 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/48.d0 1171 cm(2,8)=(32*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,2)& 1172 & *dgm10(1,2)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1173 & +dgm01(1,1)*dgm10(1,2))+32*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,1)& 1174 & +dgm01(1,1)*dgm10(1,3))+32*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1175 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1176 & *gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+8*(6*gm(1,2)& 1177 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1178 & *dgm10(2,2))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)& 1179 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 1180 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+32*gm(1,1)& 1181 & *gm(1,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*gm(1,2)& 1182 & *gm(2,2)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+16*(1*gm(1,2)& 1183 & **2+3*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1184 & +16*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 1185 & /8.d0 1186 cm(3,8)=(16*gm(3,3)**2*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1187 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 1188 & +32*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1189 & +64*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1190 & +32*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1191 & +32*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1192 & *dgm10(2,3))+32*gm(1,1)*gm(1,2)*dgm01(3,3)*dgm10(3,3)+4*(6*gm(1,3)& 1193 & *gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1194 & *dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)& 1195 & +dgm01(1,2)*dgm10(3,3))+8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1196 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+16*(1*gm(1,2)& 1197 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1198 & *dgm10(3,3))+32*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1199 & *dgm10(3,3)))/8.d0 1200 cm(4,8)=(64*gm(1,3)*gm(3,3)*dgm01(1,2)*dgm10(1,2)+16*gm(2,3)& 1201 & *gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+32*(3*gm(1,3)& 1202 & *gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,3)*dgm10(1,3)+8*(1*gm(2,3)& 1203 & **2+3*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1204 & +8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)& 1205 & +dgm01(1,2)*dgm10(1,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1206 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+64*gm(1,1)*gm(1,3)& 1207 & *dgm01(2,3)*dgm10(2,3)+8*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1208 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1209 & *gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,3)& 1210 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1211 & +8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)& 1212 & +dgm01(1,3)*dgm10(2,3))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1213 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+4*(3*gm(1,3)*gm(2,2)& 1214 & +gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 1215 & +8*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 1216 & +dgm01(1,2)*dgm10(3,3))+8*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(3,3)& 1217 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+16*gm(1,1)*gm(1,2)*(dgm01(3,3)& 1218 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/4.d0 1219 cm(5,8)=(48*gm(1,3)*gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 1220 & *dgm10(1,2))+48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,3)& 1221 & *dgm10(1,3)+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,3)& 1222 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+12*(6*gm(1,3)**2-2*gm(1,1)& 1223 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+24*(1*gm(1,3)& 1224 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1225 & +96*gm(1,1)*gm(1,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 1226 & +12*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 1227 & +dgm01(1,1)*dgm10(3,3))+48*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(1,3)& 1228 & +dgm01(1,3)*dgm10(3,3)))/12.d0 1229 cm(6,8)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,2)*dgm10(1,2)& 1230 & +12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1231 & +dgm01(1,1)*dgm10(1,2))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1232 & *(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+24*(6*gm(1,2)& 1233 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1234 & *dgm10(1,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)& 1235 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,3)*(dgm01(2,3)& 1236 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(1*gm(1,2)**2+3*gm(1,1)& 1237 & *gm(2,2))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+48*gm(1,1)& 1238 & *gm(1,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/12.d0 1239 cm(7,8)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,2)*dgm10(1,2)& 1240 & +96*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1241 & +12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1242 & +dgm01(1,2)*dgm10(2,2))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1243 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(-2*gm(1,3)& 1244 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1245 & *dgm10(2,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)& 1246 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1247 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+12*(1*gm(1,2)**2+3*gm(1,1)& 1248 & *gm(2,2))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/12.d0 1249 cm(8,8)=(32*gm(3,3)**2*dgm01(1,2)*dgm10(1,2)+32*(1*gm(2,3)**2+3*gm(2,2)& 1250 & *gm(3,3))*dgm01(1,3)*dgm10(1,3)+64*gm(2,3)*gm(3,3)*(dgm01(1,3)& 1251 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+32*(1*gm(1,3)**2+3*gm(1,1)& 1252 & *gm(3,3))*dgm01(2,3)*dgm10(2,3)+64*gm(1,3)*gm(3,3)*(dgm01(2,3)& 1253 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+32*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1254 & *gm(3,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(1*gm(1,2)& 1255 & **2+3*gm(1,1)*gm(2,2))*dgm01(3,3)*dgm10(3,3)+8*(6*gm(1,3)*gm(2,3)& 1256 & -2*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1257 & +16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)& 1258 & +dgm01(1,3)*dgm10(3,3))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1259 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/4.d0 1260 cm(9,8)=(64*gm(2,3)*gm(3,3)*dgm01(1,2)*dgm10(1,2)+64*gm(2,2)& 1261 & *gm(2,3)*dgm01(1,3)*dgm10(1,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))& 1262 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+16*(1*gm(2,3)& 1263 & **2+3*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1264 & +16*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1265 & +8*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)& 1266 & +dgm01(1,3)*dgm10(2,2))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1267 & *dgm01(2,3)*dgm10(2,3)+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1268 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,3)& 1269 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1270 & *dgm10(2,3))+16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 1271 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 1272 & *gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(1*gm(1,3)& 1273 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1274 & +8*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 1275 & +dgm01(1,2)*dgm10(3,3))+16*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,3)& 1276 & +dgm01(1,3)*dgm10(3,3))+4*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1277 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+8*(1*gm(1,2)**2+3*gm(1,1)& 1278 & *gm(2,2))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/4.d0 1279 cm(10,8)=(48*gm(3,3)**2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1280 & +72*gm(1,1)*gm(2,3)*dgm01(3,3)*dgm10(3,3)+24*gm(1,3)**2*(dgm01(3,3)& 1281 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3))+gm(3,3)*(192*gm(2,3)*dgm01(1,3)& 1282 & *dgm10(1,3)+48*gm(1,2)*dgm01(3,3)*dgm10(1,3)+72*gm(1,1)*dgm01(3,3)& 1283 & *dgm10(2,3)+48*gm(1,2)*dgm01(1,3)*dgm10(3,3)+72*gm(1,1)*dgm01(2,3)& 1284 & *dgm10(3,3)+gm(1,3)*(48*dgm01(3,3)*dgm10(1,2)+96*dgm01(2,3)*dgm10(1,3)& 1285 & +96*dgm01(1,3)*dgm10(2,3)+48*dgm01(1,2)*dgm10(3,3)))+gm(1,3)& 1286 & *(24*gm(1,2)*dgm01(3,3)*dgm10(3,3)+96*gm(2,3)*(dgm01(3,3)*dgm10(1,3)& 1287 & +dgm01(1,3)*dgm10(3,3))))/12.d0 1288 cm(11,8)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1289 & +dgm01(1,2)*dgm10(2,2))+192*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,3)& 1290 & +dgm01(1,3)*dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1291 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+96*gm(1,2)*gm(2,2)& 1292 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/48.d0 1293 cm(12,8)=(16*gm(3,3)**2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1294 & +32*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1295 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 1296 & +64*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1297 & +32*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1298 & *dgm10(2,3))+32*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1299 & *dgm10(2,3))+32*gm(1,2)*gm(2,2)*dgm01(3,3)*dgm10(3,3)+4*(6*gm(2,3)& 1300 & **2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1301 & +32*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 1302 & +4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)& 1303 & +dgm01(2,2)*dgm10(3,3))+16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1304 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))+8*(-2*gm(1,3)& 1305 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1306 & *dgm10(3,3)))/8.d0 1307 cm(13,8)=(48*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1308 & *dgm10(2,2))+24*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,2)& 1309 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 1310 & *gm(2,3))*dgm01(2,3)*dgm10(2,3)+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))& 1311 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+96*gm(2,2)*gm(2,3)& 1312 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+24*(1*gm(1,3)& 1313 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1314 & *dgm10(2,3))+12*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)& 1315 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1316 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1317 cm(14,8)=(48*gm(3,3)**2*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1318 & +72*gm(1,3)*gm(2,2)*dgm01(3,3)*dgm10(3,3)+24*gm(2,3)**2*(dgm01(3,3)& 1319 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+gm(3,3)*(192*gm(1,3)*dgm01(2,3)& 1320 & *dgm10(2,3)+48*gm(1,2)*dgm01(3,3)*dgm10(2,3)+48*gm(1,2)*dgm01(2,3)& 1321 & *dgm10(3,3)+gm(2,3)*(48*dgm01(3,3)*dgm10(1,2)+96*dgm01(2,3)*dgm10(1,3)& 1322 & +96*dgm01(1,3)*dgm10(2,3)+48*dgm01(1,2)*dgm10(3,3))+72*gm(2,2)& 1323 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))+gm(2,3)*(24*gm(1,2)& 1324 & *dgm01(3,3)*dgm10(3,3)+96*gm(1,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1325 & *dgm10(3,3))))/12.d0 1326 cm(15,8)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(3,3)& 1327 & *dgm10(3,3)+96*gm(3,3)**2*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1328 & +192*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 1329 & +192*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))& 1330 & /48.d0 1331 cm(1,9)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,2)& 1332 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 1333 & *gm(2,2))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+96*gm(1,1)& 1334 & *gm(1,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+192*gm(1,1)& 1335 & *gm(1,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)))/48.d0 1336 cm(2,9)=(64*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(1,2)*dgm10(1,2)& 1337 & +32*gm(2,2)*gm(2,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)*dgm10(1,2))& 1338 & +16*gm(2,2)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1339 & +64*gm(1,2)*gm(2,2)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1340 & +32*gm(1,1)*gm(1,3)*dgm01(2,2)*dgm10(2,2)+4*(-2*gm(1,3)*gm(2,2)& 1341 & +6*gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 1342 & +8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)& 1343 & +dgm01(1,2)*dgm10(2,2))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1344 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 1345 & *gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+32*gm(1,2)& 1346 & *gm(2,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)& 1347 & **2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1348 & +32*gm(1,1)*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))& 1349 & /8.d0 1350 cm(3,9)=(32*gm(2,3)*gm(3,3)*(dgm01(1,2)*dgm10(1,1)+dgm01(1,1)& 1351 & *dgm10(1,2))+32*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(1,3)& 1352 & *dgm10(1,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 1353 & +dgm01(1,1)*dgm10(1,3))+32*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1354 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+16*gm(1,3)*gm(3,3)& 1355 & *(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+16*(1*gm(1,3)& 1356 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1357 & +8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 1358 & +dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1359 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(6*gm(1,2)*gm(1,3)& 1360 & -2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1361 & +4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1362 & *dgm10(3,3))+16*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1363 & *dgm10(3,3))+32*gm(1,1)*gm(1,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1364 & *dgm10(3,3)))/8.d0 1365 cm(4,9)=(32*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(1,2)& 1366 & *dgm10(1,2)+8*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(1,2)*dgm10(1,1)& 1367 & +dgm01(1,1)*dgm10(1,2))+64*gm(1,2)*gm(2,2)*dgm01(1,3)*dgm10(1,3)& 1368 & +16*gm(2,2)*gm(2,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1369 & +16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)& 1370 & +dgm01(1,2)*dgm10(1,3))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1371 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+4*(1*gm(1,3)*gm(2,3)& 1372 & +3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 1373 & +8*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1374 & *dgm10(2,2))+8*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)& 1375 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+64*gm(1,1)*gm(1,2)*dgm01(2,3)& 1376 & *dgm10(2,3)+8*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 1377 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+8*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 1378 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,2)& 1379 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1380 & *dgm10(2,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,3)& 1381 & +dgm01(1,3)*dgm10(2,3))+16*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 1382 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+16*gm(1,1)*gm(1,3)*(dgm01(2,3)& 1383 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/4.d0 1384 cm(5,9)=(24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(1,2)& 1385 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 1386 & *gm(2,2))*dgm01(1,3)*dgm10(1,3)+12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 1387 & *gm(2,3))*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+24*(6*gm(1,2)& 1388 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1389 & *dgm10(1,3))+12*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(2,2)& 1390 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+48*gm(1,1)*gm(1,3)*(dgm01(2,2)& 1391 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 1392 & *gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)& 1393 & *gm(1,2)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3)))/12.d0 1394 cm(6,9)=(48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(1,2)& 1395 & *dgm10(1,2)+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(1,2)& 1396 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,2))+48*gm(1,2)*gm(2,2)*(dgm01(1,3)& 1397 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 1398 & *gm(2,2))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(1*gm(1,2)& 1399 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1400 & *dgm10(2,2))+48*gm(1,1)*gm(1,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1401 & *dgm10(2,2))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 1402 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+96*gm(1,1)*gm(1,2)*(dgm01(2,3)& 1403 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))/12.d0 1404 cm(7,9)=(192*gm(2,2)*gm(2,3)*dgm01(1,2)*dgm10(1,2)+48*gm(2,2)& 1405 & **2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+24*(1*gm(1,2)& 1406 & *gm(1,3)+3*gm(1,1)*gm(2,3))*dgm01(2,2)*dgm10(2,2)+24*(3*gm(1,3)& 1407 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1408 & +12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,2)& 1409 & +dgm01(1,2)*dgm10(2,2))+48*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,3)& 1410 & +dgm01(1,3)*dgm10(2,2))+96*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(1,2)& 1411 & +dgm01(1,2)*dgm10(2,3))+24*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*(dgm01(2,3)& 1412 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/12.d0 1413 cm(8,9)=(64*gm(2,3)*gm(3,3)*dgm01(1,2)*dgm10(1,2)+64*gm(2,2)& 1414 & *gm(2,3)*dgm01(1,3)*dgm10(1,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))& 1415 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+16*(1*gm(2,3)& 1416 & **2+3*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1417 & +16*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1418 & +8*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)& 1419 & +dgm01(1,3)*dgm10(2,2))+32*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1420 & *dgm01(2,3)*dgm10(2,3)+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1421 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(1*gm(1,3)& 1422 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1423 & *dgm10(2,3))+16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 1424 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 1425 & *gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+8*(1*gm(1,3)& 1426 & **2+3*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1427 & +8*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 1428 & +dgm01(1,2)*dgm10(3,3))+16*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,3)& 1429 & +dgm01(1,3)*dgm10(3,3))+4*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1430 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+8*(1*gm(1,2)**2+3*gm(1,1)& 1431 & *gm(2,2))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/4.d0 1432 cm(9,9)=(32*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*dgm01(1,2)*dgm10(1,2)& 1433 & +32*gm(2,2)**2*dgm01(1,3)*dgm10(1,3)+64*gm(2,2)*gm(2,3)*(dgm01(1,3)& 1434 & *dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+8*(1*gm(1,3)**2+3*gm(1,1)& 1435 & *gm(3,3))*dgm01(2,2)*dgm10(2,2)+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1436 & *gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+8*(-2*gm(1,3)& 1437 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1438 & *dgm10(2,2))+32*(1*gm(1,2)**2+3*gm(1,1)*gm(2,2))*dgm01(2,3)*dgm10(2,3)& 1439 & +32*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)& 1440 & +dgm01(1,2)*dgm10(2,3))+64*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(1,3)& 1441 & +dgm01(1,3)*dgm10(2,3))+16*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1442 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/4.d0 1443 cm(10,9)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 1444 & +96*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1445 & +48*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1446 & +24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 1447 & +dgm01(1,3)*dgm10(2,3))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1448 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+12*(-2*gm(1,3)& 1449 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1450 & *dgm10(3,3))+12*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)& 1451 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 1452 & *gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1453 cm(11,9)=(48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(2,2)& 1454 & *dgm10(2,2)+192*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1455 & *dgm10(2,2))+96*gm(2,2)**2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1456 & *dgm10(2,2))+192*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1457 & *dgm10(2,3)))/48.d0 1458 cm(12,9)=(32*gm(1,3)*gm(3,3)*dgm01(2,2)*dgm10(2,2)+32*gm(2,3)& 1459 & *gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+4*(6*gm(2,3)& 1460 & **2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1461 & +64*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(2,3)*dgm10(2,3)& 1462 & +32*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1463 & *dgm10(2,3))+64*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1464 & *dgm10(2,3))+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)& 1465 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1466 & *gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+32*gm(2,2)& 1467 & *gm(2,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+16*gm(2,2)& 1468 & **2*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+4*(-2*gm(1,3)& 1469 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1470 & *dgm10(3,3))+32*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1471 & *dgm10(3,3)))/8.d0 1472 cm(13,9)=(72*gm(1,2)*gm(3,3)*dgm01(2,2)*dgm10(2,2)+24*gm(2,3)& 1473 & **2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*gm(2,2)& 1474 & **2*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(2,3)*(24*gm(1,3)& 1475 & *dgm01(2,2)*dgm10(2,2)+gm(2,2)*(96*dgm01(2,3)*dgm10(1,2)+48*dgm01(2,2)& 1476 & *dgm10(1,3)+48*dgm01(1,3)*dgm10(2,2)+96*dgm01(1,2)*dgm10(2,3))& 1477 & +96*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))+gm(2,2)& 1478 & *(72*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+192*gm(1,2)& 1479 & *dgm01(2,3)*dgm10(2,3)+48*gm(1,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1480 & *dgm10(2,3))))/12.d0 1481 cm(14,9)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(2,3)& 1482 & *dgm10(2,3)+96*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1483 & *dgm10(2,3))+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)& 1484 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+48*gm(1,3)*gm(3,3)*(dgm01(2,3)& 1485 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+24*(1*gm(2,3)**2+3*gm(2,2)& 1486 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(2,2)& 1487 & *gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+12*(1*gm(1,3)& 1488 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1489 & *dgm10(3,3))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)& 1490 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1491 cm(15,9)=(192*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1492 & *dgm10(3,3))+24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)& 1493 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+96*gm(1,3)*gm(3,3)*(dgm01(3,3)& 1494 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)& 1495 & *gm(3,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/48.d0 1496 cm(1,10)=(72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 1497 & +dgm01(1,1)*dgm10(1,3))+288*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,1)& 1498 & +dgm01(1,1)*dgm10(3,3)))/144.d0 1499 cm(2,10)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,1)& 1500 & +dgm01(1,1)*dgm10(1,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1501 & *(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+12*(6*gm(1,3)& 1502 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1503 & +12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)& 1504 & +dgm01(1,1)*dgm10(3,3))+48*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))& 1505 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,1)*gm(1,3)& 1506 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 1507 cm(3,10)=2*gm(3,3)**2*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1508 & +gm(3,3)*(gm(1,3)*(2*dgm01(3,3)*dgm10(1,1)+16*dgm01(1,3)*dgm10(1,3)& 1509 & +2*dgm01(1,1)*dgm10(3,3))+5*gm(1,1)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1510 & *dgm10(3,3)))+gm(1,3)*(4*gm(1,1)*dgm01(3,3)*dgm10(3,3)+5*gm(1,3)& 1511 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 1512 cm(4,10)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(1,3)& 1513 & *dgm10(1,3)+48*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)& 1514 & *dgm10(1,3))+96*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)& 1515 & *dgm10(1,3))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 1516 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+12*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1517 & *gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+24*(1*gm(1,3)& 1518 & **2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1519 & +24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)& 1520 & +dgm01(1,3)*dgm10(3,3))+48*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(2,3)& 1521 & +dgm01(2,3)*dgm10(3,3)))/12.d0 1522 cm(5,10)=(72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 1523 & +144*gm(1,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,1)+dgm01(1,1)*dgm10(1,3))& 1524 & +36*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1525 & *dgm10(3,3))+144*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1526 & *dgm10(3,3)))/36.d0 1527 cm(6,10)=(36*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(1,3)& 1528 & *dgm10(1,1)+dgm01(1,1)*dgm10(1,3))+36*(6*gm(1,3)**2-2*gm(1,1)& 1529 & *gm(3,3))*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))+36*(1*gm(1,2)& 1530 & *gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1531 & *dgm10(3,3))+144*gm(1,1)*gm(1,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1532 & *dgm10(3,3)))/36.d0 1533 cm(7,10)=(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(1,3)*dgm10(1,2)& 1534 & +dgm01(1,2)*dgm10(1,3))+(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1535 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(-2*gm(1,3)*gm(2,2)& 1536 & +6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1537 & +(1*gm(1,2)*gm(1,3)+3*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)& 1538 & +dgm01(2,2)*dgm10(3,3)) 1539 cm(8,10)=(48*gm(3,3)**2*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1540 & +72*gm(1,1)*gm(2,3)*dgm01(3,3)*dgm10(3,3)+24*gm(1,3)**2*(dgm01(3,3)& 1541 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3))+gm(3,3)*(192*gm(2,3)*dgm01(1,3)& 1542 & *dgm10(1,3)+48*gm(1,2)*dgm01(3,3)*dgm10(1,3)+72*gm(1,1)*dgm01(3,3)& 1543 & *dgm10(2,3)+48*gm(1,2)*dgm01(1,3)*dgm10(3,3)+72*gm(1,1)*dgm01(2,3)& 1544 & *dgm10(3,3)+gm(1,3)*(48*dgm01(3,3)*dgm10(1,2)+96*dgm01(2,3)*dgm10(1,3)& 1545 & +96*dgm01(1,3)*dgm10(2,3)+48*dgm01(1,2)*dgm10(3,3)))+gm(1,3)& 1546 & *(24*gm(1,2)*dgm01(3,3)*dgm10(3,3)+96*gm(2,3)*(dgm01(3,3)*dgm10(1,3)& 1547 & +dgm01(1,3)*dgm10(3,3))))/12.d0 1548 cm(9,10)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(1,3)*dgm10(1,3)& 1549 & +96*gm(2,3)*gm(3,3)*(dgm01(1,3)*dgm10(1,2)+dgm01(1,2)*dgm10(1,3))& 1550 & +48*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1551 & +24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 1552 & +dgm01(1,3)*dgm10(2,3))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1553 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+12*(-2*gm(1,3)& 1554 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1555 & *dgm10(3,3))+12*(1*gm(1,3)**2+3*gm(1,1)*gm(3,3))*(dgm01(3,3)& 1556 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+24*(1*gm(1,2)*gm(1,3)+3*gm(1,1)& 1557 & *gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1558 cm(10,10)=(288*gm(3,3)**2*dgm01(1,3)*dgm10(1,3)+72*(1*gm(1,3)& 1559 & **2+3*gm(1,1)*gm(3,3))*dgm01(3,3)*dgm10(3,3)+288*gm(1,3)*gm(3,3)& 1560 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))/36.d0 1561 cm(11,10)=(4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)& 1562 & +dgm01(1,3)*dgm10(2,2))+4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1563 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/8.d0 1564 cm(12,10)=(48*gm(3,3)**2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1565 & +192*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 1566 & +24*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(3,3)*dgm10(3,3)& 1567 & +12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1568 & *dgm10(3,3))+48*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1569 & *dgm10(3,3))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)& 1570 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 1571 cm(13,10)=(144*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1572 & *dgm10(2,2))+36*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)& 1573 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+36*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1574 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+36*(-2*gm(1,3)& 1575 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1576 & *dgm10(3,3)))/36.d0 1577 cm(14,10)=(144*gm(3,3)**2*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1578 & *dgm10(2,3))+72*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(3,3)& 1579 & *dgm10(3,3)+144*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1580 & *dgm10(3,3))+144*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1581 & *dgm10(3,3)))/36.d0 1582 cm(15,10)=gm(3,3)*(4*gm(1,3)*dgm01(3,3)*dgm10(3,3)+2*gm(3,3)& 1583 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 1584 cm(1,11)=((6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,2)*dgm10(1,1)& 1585 & +dgm01(1,1)*dgm10(2,2)))/4.d0 1586 cm(2,11)=(48*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(2,2)*dgm10(2,2)& 1587 & +96*gm(2,2)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 1588 & +384*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2)))& 1589 & /96.d0 1590 cm(3,11)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 1591 & +dgm01(1,1)*dgm10(2,2))+4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1592 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(6*gm(1,2)**2-2*gm(1,1)& 1593 & *gm(2,2))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/4.d0 1594 cm(4,11)=(96*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1595 & *dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)& 1596 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+192*gm(1,2)*gm(2,2)*(dgm01(2,2)& 1597 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+24*(6*gm(1,2)**2-2*gm(1,1)& 1598 & *gm(2,2))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/48.d0 1599 cm(5,11)=(4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,2)& 1600 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 1601 & *gm(2,2))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2)))/8.d0 1602 cm(6,11)=(288*gm(1,2)*gm(2,2)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1603 & *dgm10(2,2))+72*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,2)& 1604 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2)))/144.d0 1605 cm(7,11)=gm(2,2)*(4*gm(1,2)*dgm01(2,2)*dgm10(2,2)+2*gm(2,2)& 1606 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))) 1607 cm(8,11)=(24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1608 & +dgm01(1,2)*dgm10(2,2))+192*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,3)& 1609 & +dgm01(1,3)*dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1610 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+96*gm(1,2)*gm(2,2)& 1611 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/48.d0 1612 cm(9,11)=(48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(2,2)& 1613 & *dgm10(2,2)+192*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1614 & *dgm10(2,2))+96*gm(2,2)**2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1615 & *dgm10(2,2))+192*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1616 & *dgm10(2,3)))/48.d0 1617 cm(10,11)=(4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)& 1618 & +dgm01(1,3)*dgm10(2,2))+4*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1619 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/8.d0 1620 cm(11,11)=2*gm(2,2)**2*dgm01(2,2)*dgm10(2,2) 1621 cm(12,11)=(48*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(2,2)*dgm10(2,2)& 1622 & +384*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1623 & +96*gm(2,2)**2*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 1624 & /96.d0 1625 cm(13,11)=gm(2,2)*(4*gm(2,3)*dgm01(2,2)*dgm10(2,2)+2*gm(2,2)& 1626 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 1627 cm(14,11)=(72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)& 1628 & +dgm01(2,2)*dgm10(2,3))+288*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(2,2)& 1629 & +dgm01(2,2)*dgm10(3,3)))/144.d0 1630 cm(15,11)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)& 1631 & +dgm01(2,2)*dgm10(3,3)))/4.d0 1632 cm(1,12)=((6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 1633 & +dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1634 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+(6*gm(1,2)**2-2*gm(1,1)& 1635 & *gm(2,2))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/4.d0 1636 cm(2,12)=(8*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(2,2)*dgm10(2,2)& 1637 & +4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1638 & *dgm10(2,2))+16*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 1639 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+64*gm(2,2)*gm(2,3)*(dgm01(2,3)& 1640 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+64*(3*gm(1,3)*gm(2,2)+gm(1,2)& 1641 & *gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+16*(6*gm(1,2)& 1642 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1643 & *dgm10(2,3))+16*gm(2,2)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1644 & *dgm10(3,3))+64*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1645 & *dgm10(3,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(2,2)& 1646 & +dgm01(2,2)*dgm10(3,3)))/16.d0 1647 cm(3,12)=(16*gm(3,3)**2*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)*dgm10(2,2))& 1648 & +64*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1649 & +64*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1650 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 1651 & +dgm01(1,3)*dgm10(2,3))+8*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(3,3)& 1652 & *dgm10(3,3)+4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 1653 & +dgm01(1,1)*dgm10(3,3))+16*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))& 1654 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1655 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+16*(6*gm(1,2)& 1656 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1657 & *dgm10(3,3)))/16.d0 1658 cm(4,12)=(16*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1659 & *dgm10(2,2))+32*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1660 & *dgm10(2,2))+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 1661 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+32*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 1662 & *gm(2,3))*dgm01(2,3)*dgm10(2,3)+16*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))& 1663 & *(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+32*(1*gm(1,3)& 1664 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1665 & *dgm10(2,3))+32*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)& 1666 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1667 & *gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*gm(2,2)& 1668 & *gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+8*(-2*gm(1,3)& 1669 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1670 & *dgm10(3,3))+32*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1671 & *dgm10(3,3))+4*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)*dgm10(2,3)& 1672 & +dgm01(2,3)*dgm10(3,3)))/8.d0 1673 cm(5,12)=(48*gm(1,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1674 & *dgm10(2,2))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,2)& 1675 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1676 & *gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*(6*gm(1,2)& 1677 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1678 & *dgm10(2,3))+12*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)& 1679 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 1680 & *gm(2,2))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))/24.d0 1681 cm(6,12)=(12*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,2)& 1682 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+12*(6*gm(1,3)**2-2*gm(1,1)& 1683 & *gm(3,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*(3*gm(1,3)& 1684 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1685 & +48*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)& 1686 & +dgm01(1,2)*dgm10(2,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(1,1)& 1687 & +dgm01(1,1)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(3,3)& 1688 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/24.d0 1689 cm(7,12)=(24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(2,2)& 1690 & *dgm10(2,2)+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,2)& 1691 & +dgm01(1,2)*dgm10(2,2))+192*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(1,2)& 1692 & +dgm01(1,2)*dgm10(2,3))+48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1693 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+48*gm(2,2)**2*(dgm01(3,3)& 1694 & *dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1695 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 1696 cm(8,12)=(16*gm(3,3)**2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))& 1697 & +32*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1698 & +64*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 1699 & +64*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1700 & +32*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1701 & *dgm10(2,3))+32*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1702 & *dgm10(2,3))+32*gm(1,2)*gm(2,2)*dgm01(3,3)*dgm10(3,3)+4*(6*gm(2,3)& 1703 & **2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1704 & +32*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 1705 & +4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)& 1706 & +dgm01(2,2)*dgm10(3,3))+16*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))& 1707 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))+8*(-2*gm(1,3)& 1708 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1709 & *dgm10(3,3)))/8.d0 1710 cm(9,12)=(32*gm(1,3)*gm(3,3)*dgm01(2,2)*dgm10(2,2)+32*gm(2,3)& 1711 & *gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+4*(6*gm(2,3)& 1712 & **2-2*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1713 & +64*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(2,3)*dgm10(2,3)& 1714 & +32*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1715 & *dgm10(2,3))+64*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1716 & *dgm10(2,3))+8*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)& 1717 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+16*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1718 & *gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+32*gm(2,2)& 1719 & *gm(2,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+16*gm(2,2)& 1720 & **2*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+4*(-2*gm(1,3)& 1721 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1722 & *dgm10(3,3))+32*gm(1,2)*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1723 & *dgm10(3,3)))/8.d0 1724 cm(10,12)=(48*gm(3,3)**2*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1725 & +192*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 1726 & +24*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)*gm(2,3))*dgm01(3,3)*dgm10(3,3)& 1727 & +12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1728 & *dgm10(3,3))+48*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1729 & *dgm10(3,3))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)& 1730 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 1731 cm(11,12)=(48*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(2,2)*dgm10(2,2)& 1732 & +384*gm(2,2)*gm(2,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1733 & +96*gm(2,2)**2*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 1734 & /96.d0 1735 cm(12,12)=(32*gm(3,3)**2*dgm01(2,2)*dgm10(2,2)+128*(1*gm(2,3)& 1736 & **2+3*gm(2,2)*gm(3,3))*dgm01(2,3)*dgm10(2,3)+128*gm(2,3)*gm(3,3)& 1737 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+32*gm(2,2)**2*dgm01(3,3)& 1738 & *dgm10(3,3)+8*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)& 1739 & +dgm01(2,2)*dgm10(3,3))+128*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(2,3)& 1740 & +dgm01(2,3)*dgm10(3,3)))/16.d0 1741 cm(13,12)=5*gm(2,3)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1742 & +gm(2,3)*(4*gm(3,3)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(2*dgm01(3,3)& 1743 & *dgm10(2,2)+16*dgm01(2,3)*dgm10(2,3)+2*dgm01(2,2)*dgm10(3,3)))& 1744 & +gm(2,2)*(5*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1745 & +2*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1746 cm(14,12)=2*gm(3,3)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1747 & +gm(3,3)*(gm(2,3)*(2*dgm01(3,3)*dgm10(2,2)+16*dgm01(2,3)*dgm10(2,3)& 1748 & +2*dgm01(2,2)*dgm10(3,3))+5*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1749 & *dgm10(3,3)))+gm(2,3)*(4*gm(2,2)*dgm01(3,3)*dgm10(3,3)+5*gm(2,3)& 1750 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1751 cm(15,12)=(48*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(3,3)*dgm10(3,3)& 1752 & +96*gm(3,3)**2*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))& 1753 & +384*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))& 1754 & /96.d0 1755 cm(1,13)=(4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,2)& 1756 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+4*(6*gm(1,2)**2-2*gm(1,1)& 1757 & *gm(2,2))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3)))/8.d0 1758 cm(2,13)=(24*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(2,2)& 1759 & *dgm10(2,2)+48*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1760 & *dgm10(2,2))+48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 1761 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*gm(2,2)**2*(dgm01(2,3)& 1762 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+192*gm(1,2)*gm(2,2)*(dgm01(2,3)& 1763 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 1764 & *gm(2,2))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/24.d0 1765 cm(3,13)=(48*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,1)+dgm01(1,1)& 1766 & *dgm10(2,2))+48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)& 1767 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+12*(6*gm(2,3)**2-2*gm(2,2)& 1768 & *gm(3,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+48*(-2*gm(1,3)& 1769 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1770 & *dgm10(2,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)& 1771 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+12*(6*gm(1,2)**2-2*gm(1,1)& 1772 & *gm(2,2))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 1773 cm(4,13)=(12*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 1774 & +dgm01(1,1)*dgm10(2,2))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))& 1775 & *(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+24*(3*gm(1,3)& 1776 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))& 1777 & +24*(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*dgm01(2,3)*dgm10(2,3)+48*gm(2,2)& 1778 & *gm(2,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+24*(-2*gm(1,3)& 1779 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1780 & *dgm10(2,3))+96*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1781 & *dgm10(2,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(2,3)& 1782 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/12.d0 1783 cm(5,13)=(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,2)*dgm10(1,1)& 1784 & +dgm01(1,1)*dgm10(2,2))+(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1785 & *(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+(-2*gm(1,3)*gm(2,2)& 1786 & +6*gm(1,2)*gm(2,3))*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1787 & +(6*gm(1,2)**2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1788 & *dgm10(2,3)) 1789 cm(6,13)=(36*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(2,2)& 1790 & *dgm10(1,1)+dgm01(1,1)*dgm10(2,2))+36*(6*gm(1,2)*gm(1,3)-2*gm(1,1)& 1791 & *gm(2,3))*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+144*gm(1,2)& 1792 & *gm(2,2)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))+36*(6*gm(1,2)& 1793 & **2-2*gm(1,1)*gm(2,2))*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3)))& 1794 & /36.d0 1795 cm(7,13)=(72*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*dgm01(2,2)& 1796 & *dgm10(2,2)+144*gm(2,2)*gm(2,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1797 & *dgm10(2,2))+144*gm(2,2)**2*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1798 & *dgm10(2,3))+144*gm(1,2)*gm(2,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1799 & *dgm10(2,3)))/36.d0 1800 cm(8,13)=(48*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)& 1801 & *dgm10(2,2))+24*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(2,2)& 1802 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,2))+48*(-2*gm(1,3)*gm(2,2)+6*gm(1,2)& 1803 & *gm(2,3))*dgm01(2,3)*dgm10(2,3)+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))& 1804 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+96*gm(2,2)*gm(2,3)& 1805 & *(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+24*(1*gm(1,3)& 1806 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1807 & *dgm10(2,3))+12*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)& 1808 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+48*gm(1,2)*gm(2,2)*(dgm01(3,3)& 1809 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1810 cm(9,13)=(72*gm(1,2)*gm(3,3)*dgm01(2,2)*dgm10(2,2)+24*gm(2,3)& 1811 & **2*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+48*gm(2,2)& 1812 & **2*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+gm(2,3)*(24*gm(1,3)& 1813 & *dgm01(2,2)*dgm10(2,2)+gm(2,2)*(96*dgm01(2,3)*dgm10(1,2)+48*dgm01(2,2)& 1814 & *dgm10(1,3)+48*dgm01(1,3)*dgm10(2,2)+96*dgm01(1,2)*dgm10(2,3))& 1815 & +96*gm(1,2)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))+gm(2,2)& 1816 & *(72*gm(3,3)*(dgm01(2,2)*dgm10(1,2)+dgm01(1,2)*dgm10(2,2))+192*gm(1,2)& 1817 & *dgm01(2,3)*dgm10(2,3)+48*gm(1,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)& 1818 & *dgm10(2,3))))/12.d0 1819 cm(10,13)=(144*gm(2,3)*gm(3,3)*(dgm01(2,2)*dgm10(1,3)+dgm01(1,3)& 1820 & *dgm10(2,2))+36*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)& 1821 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+36*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1822 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+36*(-2*gm(1,3)& 1823 & *gm(2,2)+6*gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1824 & *dgm10(3,3)))/36.d0 1825 cm(11,13)=gm(2,2)*(4*gm(2,3)*dgm01(2,2)*dgm10(2,2)+2*gm(2,2)& 1826 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))) 1827 cm(12,13)=5*gm(2,3)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1828 & +gm(2,3)*(4*gm(3,3)*dgm01(2,2)*dgm10(2,2)+gm(2,2)*(2*dgm01(3,3)& 1829 & *dgm10(2,2)+16*dgm01(2,3)*dgm10(2,3)+2*dgm01(2,2)*dgm10(3,3)))& 1830 & +gm(2,2)*(5*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1831 & +2*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1832 cm(13,13)=(72*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*dgm01(2,2)*dgm10(2,2)& 1833 & +288*gm(2,2)**2*dgm01(2,3)*dgm10(2,3)+288*gm(2,2)*gm(2,3)*(dgm01(2,3)& 1834 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3)))/36.d0 1835 cm(14,13)=(72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 1836 & +144*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1837 & +36*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1838 & *dgm10(3,3))+144*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1839 & *dgm10(3,3)))/36.d0 1840 cm(15,13)=(288*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1841 & *dgm10(3,3))+72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)& 1842 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/144.d0 1843 cm(1,14)=(4*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 1844 & +dgm01(1,1)*dgm10(2,3))+4*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1845 & *(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3)))/8.d0 1846 cm(2,14)=(12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 1847 & +dgm01(1,1)*dgm10(2,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1848 & *(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+12*(6*gm(1,3)& 1849 & **2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1850 & +48*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 1851 & +48*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)& 1852 & +dgm01(1,2)*dgm10(3,3))+12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))& 1853 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/24.d0 1854 cm(3,14)=(48*gm(3,3)**2*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1855 & +192*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)*dgm10(2,3))& 1856 & +24*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*dgm01(3,3)*dgm10(3,3)& 1857 & +48*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 1858 & +48*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,3)& 1859 & +dgm01(1,3)*dgm10(3,3))+12*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)& 1860 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/24.d0 1861 cm(4,14)=(24*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 1862 & +48*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)*dgm10(2,3))& 1863 & +96*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1864 & +24*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,3)& 1865 & +dgm01(1,3)*dgm10(2,3))+12*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(3,3)& 1866 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+24*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1867 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+24*(3*gm(1,3)& 1868 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 1869 & +12*(6*gm(1,2)*gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(2,3)& 1870 & +dgm01(2,3)*dgm10(3,3)))/12.d0 1871 cm(5,14)=(144*gm(1,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,1)+dgm01(1,1)& 1872 & *dgm10(2,3))+36*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 1873 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+36*(1*gm(1,3)*gm(2,3)+3*gm(1,2)& 1874 & *gm(3,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+36*(6*gm(1,2)& 1875 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1876 & *dgm10(3,3)))/36.d0 1877 cm(6,14)=(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,1)& 1878 & +dgm01(1,1)*dgm10(2,3))+(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(2,3)& 1879 & *dgm10(1,2)+dgm01(1,2)*dgm10(2,3))+(3*gm(1,3)*gm(2,2)+gm(1,2)& 1880 & *gm(2,3))*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+(6*gm(1,2)& 1881 & *gm(1,3)-2*gm(1,1)*gm(2,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1882 & *dgm10(3,3)) 1883 cm(7,14)=(36*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(1,2)& 1884 & +dgm01(1,2)*dgm10(2,3))+36*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1885 & *(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+144*gm(2,2)*gm(2,3)& 1886 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+36*(3*gm(1,3)& 1887 & *gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))& 1888 & /36.d0 1889 cm(8,14)=(48*gm(3,3)**2*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)*dgm10(2,3))& 1890 & +72*gm(1,3)*gm(2,2)*dgm01(3,3)*dgm10(3,3)+24*gm(2,3)**2*(dgm01(3,3)& 1891 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+gm(3,3)*(192*gm(1,3)*dgm01(2,3)& 1892 & *dgm10(2,3)+48*gm(1,2)*dgm01(3,3)*dgm10(2,3)+48*gm(1,2)*dgm01(2,3)& 1893 & *dgm10(3,3)+gm(2,3)*(48*dgm01(3,3)*dgm10(1,2)+96*dgm01(2,3)*dgm10(1,3)& 1894 & +96*dgm01(1,3)*dgm10(2,3)+48*dgm01(1,2)*dgm10(3,3))+72*gm(2,2)& 1895 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))+gm(2,3)*(24*gm(1,2)& 1896 & *dgm01(3,3)*dgm10(3,3)+96*gm(1,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1897 & *dgm10(3,3))))/12.d0 1898 cm(9,14)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(2,3)& 1899 & *dgm10(2,3)+96*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(1,2)+dgm01(1,2)& 1900 & *dgm10(2,3))+12*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)& 1901 & *dgm10(1,3)+dgm01(1,3)*dgm10(2,3))+48*gm(1,3)*gm(3,3)*(dgm01(2,3)& 1902 & *dgm10(2,2)+dgm01(2,2)*dgm10(2,3))+24*(1*gm(2,3)**2+3*gm(2,2)& 1903 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+48*gm(2,2)& 1904 & *gm(2,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+12*(1*gm(1,3)& 1905 & *gm(2,3)+3*gm(1,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1906 & *dgm10(3,3))+24*(3*gm(1,3)*gm(2,2)+gm(1,2)*gm(2,3))*(dgm01(3,3)& 1907 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/12.d0 1908 cm(10,14)=(144*gm(3,3)**2*(dgm01(2,3)*dgm10(1,3)+dgm01(1,3)& 1909 & *dgm10(2,3))+72*(1*gm(1,3)*gm(2,3)+3*gm(1,2)*gm(3,3))*dgm01(3,3)& 1910 & *dgm10(3,3)+144*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)& 1911 & *dgm10(3,3))+144*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1912 & *dgm10(3,3)))/36.d0 1913 cm(11,14)=(72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(2,3)*dgm10(2,2)& 1914 & +dgm01(2,2)*dgm10(2,3))+288*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(2,2)& 1915 & +dgm01(2,2)*dgm10(3,3)))/144.d0 1916 cm(12,14)=2*gm(3,3)**2*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1917 & +gm(3,3)*(gm(2,3)*(2*dgm01(3,3)*dgm10(2,2)+16*dgm01(2,3)*dgm10(2,3)& 1918 & +2*dgm01(2,2)*dgm10(3,3))+5*gm(2,2)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1919 & *dgm10(3,3)))+gm(2,3)*(4*gm(2,2)*dgm01(3,3)*dgm10(3,3)+5*gm(2,3)& 1920 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1921 cm(13,14)=(72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(2,3)*dgm10(2,3)& 1922 & +144*gm(2,3)*gm(3,3)*(dgm01(2,3)*dgm10(2,2)+dgm01(2,2)*dgm10(2,3))& 1923 & +36*(1*gm(2,3)**2+3*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1924 & *dgm10(3,3))+144*gm(2,2)*gm(2,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)& 1925 & *dgm10(3,3)))/36.d0 1926 cm(14,14)=(288*gm(3,3)**2*dgm01(2,3)*dgm10(2,3)+72*(1*gm(2,3)& 1927 & **2+3*gm(2,2)*gm(3,3))*dgm01(3,3)*dgm10(3,3)+288*gm(2,3)*gm(3,3)& 1928 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/36.d0 1929 cm(15,14)=gm(3,3)*(4*gm(2,3)*dgm01(3,3)*dgm10(3,3)+2*gm(3,3)& 1930 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1931 cm(1,15)=((6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 1932 & +dgm01(1,1)*dgm10(3,3)))/4.d0 1933 cm(2,15)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,1)& 1934 & +dgm01(1,1)*dgm10(3,3))+4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1935 & *(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))+(6*gm(1,3)**2-2*gm(1,1)& 1936 & *gm(3,3))*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/4.d0 1937 cm(3,15)=(48*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*dgm01(3,3)*dgm10(3,3)& 1938 & +96*gm(3,3)**2*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)*dgm10(3,3))& 1939 & +384*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))& 1940 & /96.d0 1941 cm(4,15)=(96*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1942 & *dgm10(3,3))+192*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1943 & *dgm10(3,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)& 1944 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+24*(6*gm(1,3)**2-2*gm(1,1)& 1945 & *gm(3,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/48.d0 1946 cm(5,15)=(288*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,1)+dgm01(1,1)& 1947 & *dgm10(3,3))+72*(6*gm(1,3)**2-2*gm(1,1)*gm(3,3))*(dgm01(3,3)& 1948 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3)))/144.d0 1949 cm(6,15)=(4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*(dgm01(3,3)& 1950 & *dgm10(1,1)+dgm01(1,1)*dgm10(3,3))+4*(6*gm(1,3)**2-2*gm(1,1)& 1951 & *gm(3,3))*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3)))/8.d0 1952 cm(7,15)=(4*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(1,2)& 1953 & +dgm01(1,2)*dgm10(3,3))+4*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))& 1954 & *(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3)))/8.d0 1955 cm(8,15)=(48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)*gm(3,3))*dgm01(3,3)& 1956 & *dgm10(3,3)+96*gm(3,3)**2*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)*dgm10(3,3))& 1957 & +192*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))& 1958 & +192*gm(1,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))& 1959 & /48.d0 1960 cm(9,15)=(192*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(1,2)+dgm01(1,2)& 1961 & *dgm10(3,3))+24*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)& 1962 & *dgm10(1,3)+dgm01(1,3)*dgm10(3,3))+96*gm(1,3)*gm(3,3)*(dgm01(3,3)& 1963 & *dgm10(2,2)+dgm01(2,2)*dgm10(3,3))+48*(6*gm(1,3)*gm(2,3)-2*gm(1,2)& 1964 & *gm(3,3))*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/48.d0 1965 cm(10,15)=gm(3,3)*(4*gm(1,3)*dgm01(3,3)*dgm10(3,3)+2*gm(3,3)& 1966 & *(dgm01(3,3)*dgm10(1,3)+dgm01(1,3)*dgm10(3,3))) 1967 cm(11,15)=((6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)*dgm10(2,2)& 1968 & +dgm01(2,2)*dgm10(3,3)))/4.d0 1969 cm(12,15)=(48*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*dgm01(3,3)*dgm10(3,3)& 1970 & +96*gm(3,3)**2*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)*dgm10(3,3))& 1971 & +384*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))& 1972 & /96.d0 1973 cm(13,15)=(288*gm(2,3)*gm(3,3)*(dgm01(3,3)*dgm10(2,2)+dgm01(2,2)& 1974 & *dgm10(3,3))+72*(6*gm(2,3)**2-2*gm(2,2)*gm(3,3))*(dgm01(3,3)& 1975 & *dgm10(2,3)+dgm01(2,3)*dgm10(3,3)))/144.d0 1976 cm(14,15)=gm(3,3)*(4*gm(2,3)*dgm01(3,3)*dgm10(3,3)+2*gm(3,3)& 1977 & *(dgm01(3,3)*dgm10(2,3)+dgm01(2,3)*dgm10(3,3))) 1978 cm(15,15)=2*gm(3,3)**2*dgm01(3,3)*dgm10(3,3) 1979 end if 1980 ! 1981 !contraction to output scalar 1982 ! 1983 e2nl=0.d0 1984 do jj=1,((rank+3)*(rank+4))/2 1985 tmp(:)=0.d0 1986 do ii=1,((rank+3)*(rank+4))/2 1987 tmp(:)=tmp(:)+aa(:,ii)*cm(ii,jj) 1988 end do 1989 e2nl=e2nl+tmp(1)*bb(1,jj)+tmp(2)*bb(2,jj) 1990 end do 1991 e2nl=0.25d0*e2nl 1992 1993 ABI_FREE(cm) 1994 1995 end subroutine contstr25