Rosetta 3.5
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
NatbiasStrandPairPotential.cc
Go to the documentation of this file.
1 // -*- mode:c++;tab-width:2;indent-tabs-mode:t;show-trailing-whitespace:t;rm-trailing-spaces:t -*-
2 // vi: set ts=2 noet:
3 //
4 // (c) Copyright Rosetta Commons Member Institutions.
5 // (c) This file is part of the Rosetta software suite and is made available under license.
6 // (c) The Rosetta software is developed by the contributing members of the Rosetta Commons.
7 // (c) For more information, see http://www.rosettacommons.org. Questions about this can be
8 // (c) addressed to University of Washington UW TechTransfer, email: license@u.washington.edu.
9 
10 /// @file
11 /// @brief
12 /// @author Nobuyasu Koga ( nobuyasu@uw.edu )
13 
14 // Unit header
15 #include <protocols/fldsgn/potentials/sspot/NatbiasStrandPairPotential.hh>
16 
17 // Project headers
18 #include <core/pose/Pose.hh>
19 // AUTO-REMOVED #include <core/conformation/Conformation.hh>
20 #include <core/scoring/Energies.hh>
21 #include <core/scoring/EnergyGraph.hh>
22 #include <protocols/fldsgn/topology/SS_Info2.hh>
23 #include <protocols/fldsgn/topology/BB_Pos.hh>
24 // AUTO-REMOVED #include <protocols/fldsgn/topology/DimerPairing.hh>
25 #include <protocols/fldsgn/potentials/sspot/util.hh>
26 #include <protocols/fldsgn/topology/StrandPairing.hh>
27 
28 #include <core/pose/symmetry/util.hh>
29 // AUTO-REMOVED #include <core/scoring/symmetry/SymmetricEnergies.hh>
30 
31 #include <basic/Tracer.hh>
32 #include <basic/database/open.hh>
33 
34 // Numeric
35 #include <numeric/conversions.hh>
36 
37 // utility
38 #include <utility/io/izstream.hh>
39 
40 // C++ headers
41 #include <cmath>
42 #include <iostream>
43 
44 #include <ObjexxFCL/format.hh>
45 
46 #include <core/scoring/symmetry/SymmetricEnergies.fwd.hh>
47 #include <utility/vector1.hh>
48 #include <ObjexxFCL/FArray3D.hh>
49 
50 //Auto Headers
51 #include <core/conformation/Residue.hh>
52 #include <core/kinematics/Jump.hh>
53 
54 
55 
56 
57 static basic::Tracer TR("protocols.fldsgn.potentials.sspot.NatbiasStrandPairPotential",basic::t_info);
58 
59 namespace protocols {
60 namespace fldsgn {
61 namespace potentials {
62 namespace sspot {
63 
64 /// @brief default constructor
65 NatbiasStrandPairPotential::NatbiasStrandPairPotential():
66  strand_dist_cutoff_( 6.5 ),
67  phithetascore_( 2, 3, 36, 36 ),
68  dotscore_( 6 ),
69  rsigma_dot_( 12, 18, 2, 2, rsigma_dot_initializer ),
70  native_spairset_( NULL )
71 {
72  load_phi_theta_bins();
73  load_dotscore_bins();
74 }
75 
76 /// @brief value constructor
77 NatbiasStrandPairPotential::NatbiasStrandPairPotential( StrandPairingSetOP const spairset ):
78  strand_dist_cutoff_( 6.5 ),
79  phithetascore_( 2, 3, 36, 36 ),
80  dotscore_( 6 ),
81  rsigma_dot_( 12, 18, 2, 2, rsigma_dot_initializer ),
82  native_spairset_( spairset )
83 {
84  load_phi_theta_bins();
85  load_dotscore_bins();
86 }
87 
88 /// @brief destructor
89 NatbiasStrandPairPotential::~NatbiasStrandPairPotential()
90 {}
91 
92 /// @brief
93 void
94 NatbiasStrandPairPotential::set_native_spairset( StrandPairingSetOP const spairset ){
95  native_spairset_ = spairset;
96 }
97 
98 /// @brief return score of phitheta
99 core::Real
100 NatbiasStrandPairPotential::calc_phithetascore( Real const phi, Real const theta ) const
101 {
102  Size iphi = static_cast< Size >( 1 + ( phi + 180.0 )/10 );
103  if ( iphi > 36 ) {
104  iphi = 36;
105  } else if ( iphi < 1 ) {
106  iphi = 1;
107  }
108  Size itheta = static_cast< Size >( 1 + ( theta/5 ) );
109  if ( itheta > 36 ) {
110  itheta = 36;
111  } else if ( itheta < 1 ) {
112  itheta = 1;
113  }
114  return phithetascore_( 2, 3, int( iphi ), int( itheta ) );
115 }
116 
117 /// @brief return dot score
118 core::Real
119 NatbiasStrandPairPotential::calc_dotscore( Real const dpall ) const
120 {
121  int idot;
122  if ( dpall > 0.33 && dpall <= 0.66 ) {
123  idot = 2;
124  } else if ( dpall > 0.66 && dpall <= 1.00 ) {
125  idot = 3;
126  } else if ( dpall > 1.00 && dpall <= 1.33 ) {
127  idot = 4;
128  } else if ( dpall > 1.33 && dpall <= 1.60 ) {
129  idot = 5;
130  } else if ( dpall > 1.60 && dpall <= 1.80 ) {
131  idot = 6;
132  } else if ( dpall > 1.80 && dpall <= 2.00 ) {
133  idot = 5;
134  } else {
135  idot = 1;
136  }
137  return dotscore_( idot );
138 }
139 
140 /// @brief return rsigma score
141 core::Real
142 NatbiasStrandPairPotential::calc_rsigmascore( Real sig, Real dist, Size const sign1, Size const sign2 ) const
143 {
144  if ( sig > 179.0 ) {
145  sig = 179.0;
146  } else if ( sig < 0.0 ) {
147  sig = 0.0;
148  }
149  Size isig = static_cast< Size >( sig / 10 ) + 1;
150  if ( dist > 6.4 ) {
151  dist = 6.4;
152  } else if ( dist < 3.5 ) {
153  dist = 3.5;
154  }
155  Size idist = static_cast< Size >( ( dist - 3.5 ) / 0.25 ) + 1;
156 
157  // FIX THIS !!!!!!!!!!!!!!
158  // The definition of dimer signs (sign1,sign2) appears inverted (1 should be 2, vice versa).
159  Real tempscore_rsigma = rsigma_dot_( int( idist ), int( isig ), 3 - sign1, 3 - sign2 );
160 
161  // Modify sigma potential to no longer give an rsigma bonus to strands in wrong register.
162  if ( sign1 == 1 & sign2 == 1 & sig < 110. & sig > 70. ) tempscore_rsigma = 0.0;
163  if ( sign1 == 1 & sign2 == 2 & (sig < 75. | sig > 95.) ) tempscore_rsigma = 0.0;
164  if ( sign1 == 2 & sign2 == 1 & (sig < 90. | sig > 110.)) tempscore_rsigma = 0.0;
165  if ( sign1 == 2 & sign2 == 2 & sig < 120. & sig > 80. ) tempscore_rsigma = 0.0;
166 
167  return tempscore_rsigma;
168 }
169 
170 /// @brief calculate sum of dot product of the co vectors of strand dimers ss1 and ss2
171 /// @brief with the vector connecting the midpoints of the dimer vectors (vdist)
172 /// @brief also determine return the sign of the dot products for each dimer
173 /// @brief to determine which direction the CO groups point
174 void
175 NatbiasStrandPairPotential::pair_dp(
176  Size const & ss1,
177  Size const & ss2,
178  BB_Pos const & bb_pos,
179  Real & dp,
180  Vector const & mid_vector,
181  Size & sign1,
182  Size & sign2
183 ) const
184 {
185  // length of C=O bond
186  static Real const dist_co_inv = { 1.0f / 1.231015f };
187 
188  // unit vector of mid_vector
189  Vector const u_midvec( mid_vector.normalized_or_zero() );
190 
191  //
192  Real dp1( 0.0 );
193  Real sdp1( 0.0 );
194  for ( Size i=ss1; i<=ss1+1; ++i ) {
195  Vector temp;
196  if ( i == ss1+1 ) {
197  temp = dist_co_inv * ( bb_pos.C(i) - bb_pos.O(i) );
198  } else {
199  temp = dist_co_inv * ( bb_pos.O(i) - bb_pos.C(i) );
200  }
201  Real const tempdot = temp.dot( u_midvec );
202  dp1 += std::abs(tempdot);
203  sdp1 += tempdot;
204  }
205  dp1 *= 0.5;
206 
207  //
208  Real dp2( 0.0 );
209  Real sdp2( 0.0 );
210  for ( Size i=ss2; i<=ss2+1; ++i ) {
211  Vector temp;
212  if ( i == ss2+1 ) {
213  temp = dist_co_inv * ( bb_pos.C(i) - bb_pos.O(i) );
214  } else {
215  temp = dist_co_inv * ( bb_pos.O(i) - bb_pos.C(i) );
216  }
217  Real const tempdot = temp.dot( u_midvec );
218  dp2 += std::abs(tempdot);
219  sdp2 += tempdot;
220  }
221  dp2 *= 0.5;
222 
223  dp = dp1 + dp2;
224 
225  //js These signs tell whether the first c=o bond vector of a dimer points
226  //js at the other dimer. sign1 = 1 means that the first c=o bond of dimer1
227  //js points at dimer2. sign2 = 1 means that dimer2 points at dimer1. When
228  //js sign1 or sign2 equals 2, that dimer points away from the other dimer.
229  sign1 = ( sdp1 > 0.0 ? 2 : 1 );
230  sign2 = ( sdp2 < 0.0 ? 2 : 1 );
231 }
232 
233 
234 
235 /// @brief
236 void
237 NatbiasStrandPairPotential::score( Pose const & pose, SS_Info2 const & ss_info, Real & ss_score ) const
238 {
239  using protocols::fldsgn::potentials::sspot::spherical;
240  using protocols::fldsgn::potentials::sspot::get_foldtree_seqsep;
241  using core::pose::symmetry::is_symmetric;
242  using core::scoring::Energies;
243  using core::scoring::EnergyGraph;
244  using core::scoring::symmetry::SymmetricEnergies;
245  using protocols::fldsgn::topology::StrandPairingOP; // defined in StrandPairing.hh
246  using protocols::fldsgn::topology::Strands;
247  using protocols::fldsgn::topology::Strand;
248  using protocols::fldsgn::topology::StrandOP;
249  typedef utility::vector1< Size > VecSize;
250 
251  ss_score = 0.0;
252  Real rsigma_score = 0.0;
253 
254  Strands const & strands( ss_info.strands() );
255  BB_Pos const & bb_pos( ss_info.bb_pos() );
256 
257  // EnergyGraph eg;
258  // if( is_symmetric( pose ) ) {
259  // Energies ene = pose.energies();
260  // eg = ene.energy_graph();
261  // } else {
262  // SymmetricEnergies ene = dynamic_cast< SymmetricEnergies const & > ( pose.energies() );
263  // eg = ene.energy_graph();
264  // }
265 
266  EnergyGraph const & eg( pose.energies().energy_graph() );
267 
268  for( Size istrand=1; istrand<=strands.size(); istrand++ ){
269 
270  StrandOP const strand( strands[ istrand ] );
271  for( Size ss1=strand->begin(); ss1<strand->end(); ss1++ ){
272 
273  for( core::graph::Graph::EdgeListConstIter iru =eg.get_node( ss1 )->const_upper_edge_list_begin(),
274  irue=eg.get_node( ss1 )->const_upper_edge_list_end(); iru != irue; ++iru ) {
275 
276  Size const ss2( (*iru)->get_second_node_ind() );
277  Size jstrand = ss_info.strand_id( ss2 );
278 
279  if( pose.residue_type( ss2 ).is_upper_terminus() ) continue;
280  if( ss_info.strand_id( ss2+1 ) == 0 || istrand == jstrand || jstrand == 0 ) continue;
281 
282  StrandPairingOP const & spair = native_spairset_->strand_pairing( istrand, jstrand );
283  if( spair->s1() == 0 || spair->s2() == 0 ) continue;
284 
285  ss_score += -0.2;
286 
287  Vector const & pt1( bb_pos.N( ss1 ) );
288  Vector const & pt2( bb_pos.C( ss1+1 ) );
289  Vector const & pt3( bb_pos.N( ss2 ) );
290  Vector const & pt4( bb_pos.C( ss2+1 ) );
291 
292  // midpoint coordinates of two dimers
293  Vector cen1 = Real( 0.5 )*( pt1 + pt2 );
294  Vector cen2 = Real( 0.5 )*( pt3 + pt4 );
295 
296  // vector between midpoints
297  Vector mid_vector = cen2 -cen1;
298  Real dist_dimers = mid_vector.length();
299 
300  if( dist_dimers > 8.0 ) continue;
301 
302  // calc phi and theta score between dimers
303  Real phi, theta;
304  spherical( pt2, pt4, phi, theta, cen1, cen2, mid_vector );
305 
306  if( spair->orient() == 'P' && theta>90 ) continue;
307  if( spair->orient() == 'A' && theta<90 ) continue;
308 
309  ss_score += -0.8;
310 
311  if( dist_dimers < strand_dist_cutoff_ ) {
312 
313  if( ! spair->is_member( ss1 ) || ! spair->is_member( ss2 ) ) continue;
314  Size pairres = spair->residue_pair( ss1 );
315 
316  if( spair->rgstr_shift() != 99 ) {
317  if ( spair->orient() == 'P' ) {
318  if( ss2 != pairres && ss2 != pairres+1 && ss2 != pairres-1 ) continue;
319  } else if ( spair->orient() == 'A' ) {
320  if( ss2 != pairres && ss2 != pairres-2 && ss2 != pairres-1 ) continue;
321  } else {
322  TR << "invalid strand pair orientaion " << spair->orient() << std::endl;
323  }
324 
325  }
326 
327  Real dpall;
328  Size sign1, sign2;
329  pair_dp( ss1, ss2, bb_pos, dpall, mid_vector, sign1, sign2 );
330  if( spair->rgstr_shift() != 99 ) {
331  // pairing wo/ register shift
332  if( (sign1 == 1 && sign2 == 2) || (sign1 == 2 && sign2 == 1) ) {
333  if( (spair->orient() == 'P' && ss2 != pairres) || (spair->orient() == 'A' && ss2 != pairres-1) ) continue;
334  }
335 
336  // pairing w/ register shift +/-1
337  if( (sign1 == 1 && sign2 == 1) || (sign1 == 2 && sign2 == 2) ) {
338  if( spair->orient() == 'P' && ss2 != pairres+1 &&
339  spair->orient() == 'P' && ss2 != pairres-1 &&
340  spair->orient() == 'A' && ss2 != pairres &&
341  spair->orient() == 'A' && ss2 != pairres-2 ) continue;
342  }
343  }
344 
345  Real phithetascore = calc_phithetascore( phi, theta );
346 
347  if ( phithetascore < 0.0 ) {
348 
349  ss_score += phithetascore;
350 
351  // calc dot product and dist score
352  Real dotscore( 0.0 );
353  dotscore = calc_dotscore( dpall );
354 
355  ss_score += dotscore;
356 
357  // rsigma score
358  Real const u21z = ( pt2 - cen1 ).normalized_or_zero().dot( mid_vector.normalized_or_zero() );
359  Real sig = numeric::conversions::degrees( numeric::arccos( u21z ) ); //std::acos( sin_cos_range( u21z ) );
360  Real tempscore_rsigma = calc_rsigmascore( sig, dist_dimers, sign1, sign2 );
361  if ( tempscore_rsigma<0.0 && dotscore<0.0 ) rsigma_score += tempscore_rsigma;
362 
363  TR.Debug << istrand << " " << jstrand << " " << ss1 << " " << ss2 << " " << sign1 << " " << sign2 << " "
364  << phithetascore << " " << dotscore << " " << tempscore_rsigma << std::endl;
365 
366  } // if phithetascore < 0.0
367  } // if ( dist_dimers < strand_dist_cutoff_ )
368  } //
369  } // for ( ss1 )
370  } // for( istrand )
371 
372  ss_score *= 0.498 * 0.75;
373  ss_score += rsigma_score * 0.1;
374 
375 } // score
376 
377 
378 /// @brief load phi/theta bins for use in secondary structure scoring
379 void
380 NatbiasStrandPairPotential::load_phi_theta_bins( String const & ss_filename )
381 {
382  using ObjexxFCL::fmt::skip;
383  typedef ObjexxFCL::FArray3D< Real > FArray3D_real;
384  FArray3D_real pts_SS( 36, 36, 3 );
385 
386  FArray1D_int iptsn( 36 );
387  for ( int itheta = 1; itheta <= 36; ++itheta ) {
388  iptsn(itheta) = 100;
389  }
390 
391  // FIXME: need equivalent to open_data_file() function here
392  utility::io::izstream SS_stream;
393  basic::database::open( SS_stream, ss_filename );
394  for ( int isep = 1; isep <= 3; ++isep ) {
395  for ( int itheta = 1; itheta <= 36; ++itheta ) {
396  for ( int iph = 1; iph <= 36; ++iph ) {
397  SS_stream >> pts_SS(itheta,iph,isep) >> skip;
398  }
399  }
400  if ( isep == 1 ) SS_stream.seek_beg();
401  }
402  SS_stream.close();
403  SS_stream.clear();
404 
405  for ( int isep = 1; isep <= 3; ++isep ) {
406  Real tot = 0.0;
407  Real totn = 0.0;
408  for ( int iph = 1; iph <= 36; ++iph ) {
409  for ( int itheta = 1; itheta <= 36; ++itheta ) {
410  pts_SS( itheta, iph, isep ) += iptsn( itheta )*0.000001f; // SMALL COUNTS CORRECTION
411  tot += pts_SS( itheta, iph, isep );
412  totn += iptsn( itheta );
413  }
414  }
415  for ( int iph = 1; iph <= 36; ++iph ) {
416  for ( int itheta = 1; itheta <= 36; ++itheta ) {
417  phithetascore_( 2, isep, iph, itheta ) = -std::log( pts_SS( itheta, iph, isep )/tot) + std::log(iptsn( itheta )/totn );
418  }
419  }
420  }
421 } // load_phi_theta_bins
422 
423 /// @brief
424 void
425 NatbiasStrandPairPotential::load_dotscore_bins()
426 {
427  // triangle-2 random numbers
428  // data idsn/56,167,278,278,167,56/
429  // sort of triangle-4 random numbers
430  FArray1D_int idsn( 6 );
431  idsn( 1 ) = 5596;
432  idsn( 2 ) = 16581;
433  idsn( 3 ) = 27823;
434  idsn( 4 ) = 27823;
435  idsn( 5 ) = 16581;
436  idsn( 6 ) = 5596;
437 
438  FArray1D_int ids( 6 );
439  ids( 1 ) = 1;
440  ids( 2 ) = 48;
441  ids( 3 ) = 368;
442  ids( 4 ) = 2378;
443  ids( 5 ) = 7141;
444  ids( 6 ) = 8904;
445 
446  Real tot = 0.0;
447  Real totn = 0.0;
448  for ( int idot = 1; idot <= 6; ++idot ) {
449  tot += ids(idot);
450  totn += idsn(idot);
451  }
452  for ( int idot = 1; idot <= 6; ++idot ) {
453  if ( ids(idot) != 0 ) {
454  dotscore_(idot) = -std::log(ids(idot)/tot) + std::log(idsn(idot)/totn);
455  } else {
456  dotscore_(idot) = 0.0;
457  }
458  }
459 }
460 
461 
462 /// @brief
463 void
464 NatbiasStrandPairPotential::rsigma_dot_initializer(
465  FArray4D_real & rsigma_dot
466 )
467 {
468  // section 12
469  //js --------------------------
470  //js new rsigma stats that take into account whether the first
471  //js c=o bond vector points towards away from the other dimer
472  rsigma_dot( 1, 1, 1, 2 ) = -1.038100;
473  rsigma_dot( 1, 2, 1, 2 ) = -1.038100;
474  rsigma_dot( 1, 3, 1, 2 ) = -1.038100;
475  rsigma_dot( 1, 4, 1, 2 ) = -1.038100;
476  rsigma_dot( 1, 5, 1, 2 ) = -1.038100;
477  rsigma_dot( 1, 6, 1, 2 ) = -1.038100;
478  rsigma_dot( 1, 7, 1, 2 ) = -1.038100;
479  rsigma_dot( 1, 8, 1, 2 ) = -2.984000;
480  rsigma_dot( 1, 9, 1, 2 ) = -3.746110;
481  rsigma_dot( 1, 10, 1, 2 ) = -3.746110;
482  rsigma_dot( 1, 11, 1, 2 ) = -2.647500;
483  rsigma_dot( 1, 12, 1, 2 ) = -1.038100;
484  rsigma_dot( 1, 13, 1, 2 ) = -1.038100;
485  rsigma_dot( 1, 14, 1, 2 ) = -1.038100;
486  rsigma_dot( 1, 15, 1, 2 ) = -1.038100;
487  rsigma_dot( 1, 16, 1, 2 ) = -1.038100;
488  rsigma_dot( 1, 17, 1, 2 ) = -1.038100;
489  rsigma_dot( 1, 18, 1, 2 ) = -1.038100;
490  rsigma_dot( 2, 1, 1, 2 ) = -0.973500;
491  rsigma_dot( 2, 2, 1, 2 ) = -0.973500;
492  rsigma_dot( 2, 3, 1, 2 ) = -0.973500;
493  rsigma_dot( 2, 4, 1, 2 ) = -0.973500;
494  rsigma_dot( 2, 5, 1, 2 ) = -0.973500;
495  rsigma_dot( 2, 6, 1, 2 ) = -0.973500;
496  rsigma_dot( 2, 7, 1, 2 ) = -2.072100;
497  rsigma_dot( 2, 8, 1, 2 ) = -0.973500;
498  rsigma_dot( 2, 9, 1, 2 ) = -4.943820;
499  rsigma_dot( 2, 10, 1, 2 ) = -5.178220;
500  rsigma_dot( 2, 11, 1, 2 ) = -4.823670;
501  rsigma_dot( 2, 12, 1, 2 ) = -2.919400;
502  rsigma_dot( 2, 13, 1, 2 ) = -2.583000;
503  rsigma_dot( 2, 14, 1, 2 ) = -0.973500;
504  rsigma_dot( 2, 15, 1, 2 ) = -0.973500;
505  rsigma_dot( 2, 16, 1, 2 ) = -0.973500;
506  rsigma_dot( 2, 17, 1, 2 ) = -0.973500;
507  rsigma_dot( 2, 18, 1, 2 ) = -0.973500;
508  rsigma_dot( 3, 1, 1, 2 ) = -0.912900;
509  rsigma_dot( 3, 2, 1, 2 ) = -0.912900;
510  rsigma_dot( 3, 3, 1, 2 ) = -0.912900;
511  rsigma_dot( 3, 4, 1, 2 ) = -0.912900;
512  rsigma_dot( 3, 5, 1, 2 ) = -0.912900;
513  rsigma_dot( 3, 6, 1, 2 ) = -0.912900;
514  rsigma_dot( 3, 7, 1, 2 ) = -2.011500;
515  rsigma_dot( 3, 8, 1, 2 ) = -4.048400;
516  rsigma_dot( 3, 9, 1, 2 ) = -5.981810;
517  rsigma_dot( 3, 10, 1, 2 ) = -7.006470;
518  rsigma_dot( 3, 11, 1, 2 ) = -6.696730;
519  rsigma_dot( 3, 12, 1, 2 ) = -4.131780;
520  rsigma_dot( 3, 13, 1, 2 ) = -0.912900;
521  rsigma_dot( 3, 14, 1, 2 ) = -0.912900;
522  rsigma_dot( 3, 15, 1, 2 ) = -0.912900;
523  rsigma_dot( 3, 16, 1, 2 ) = -0.912900;
524  rsigma_dot( 3, 17, 1, 2 ) = -0.912900;
525  rsigma_dot( 3, 18, 1, 2 ) = -0.912900;
526  rsigma_dot( 4, 1, 1, 2 ) = -0.855700;
527  rsigma_dot( 4, 2, 1, 2 ) = -0.855700;
528  rsigma_dot( 4, 3, 1, 2 ) = -0.855700;
529  rsigma_dot( 4, 4, 1, 2 ) = -0.855700;
530  rsigma_dot( 4, 5, 1, 2 ) = -0.855700;
531  rsigma_dot( 4, 6, 1, 2 ) = -2.465200;
532  rsigma_dot( 4, 7, 1, 2 ) = -3.420700;
533  rsigma_dot( 4, 8, 1, 2 ) = -3.900270;
534  rsigma_dot( 4, 9, 1, 2 ) = -6.315330;
535  rsigma_dot( 4, 10, 1, 2 ) = -8.116270;
536  rsigma_dot( 4, 11, 1, 2 ) = -8.329380;
537  rsigma_dot( 4, 12, 1, 2 ) = -5.298390;
538  rsigma_dot( 4, 13, 1, 2 ) = -1.954400;
539  rsigma_dot( 4, 14, 1, 2 ) = -0.855700;
540  rsigma_dot( 4, 15, 1, 2 ) = -0.855700;
541  rsigma_dot( 4, 16, 1, 2 ) = -0.855700;
542  rsigma_dot( 4, 17, 1, 2 ) = -0.855700;
543  rsigma_dot( 4, 18, 1, 2 ) = -0.855700;
544  rsigma_dot( 5, 1, 1, 2 ) = -0.801700;
545  rsigma_dot( 5, 2, 1, 2 ) = -0.801700;
546  rsigma_dot( 5, 3, 1, 2 ) = -0.801700;
547  rsigma_dot( 5, 4, 1, 2 ) = -0.801700;
548  rsigma_dot( 5, 5, 1, 2 ) = -0.801700;
549  rsigma_dot( 5, 6, 1, 2 ) = -1.900300;
550  rsigma_dot( 5, 7, 1, 2 ) = -3.634900;
551  rsigma_dot( 5, 8, 1, 2 ) = -5.119160;
552  rsigma_dot( 5, 9, 1, 2 ) = -6.381410;
553  rsigma_dot( 5, 10, 1, 2 ) = -9.110860;
554  rsigma_dot( 5, 11, 1, 2 ) = -9.615860;
555  rsigma_dot( 5, 12, 1, 2 ) = -6.566870;
556  rsigma_dot( 5, 13, 1, 2 ) = -2.411100;
557  rsigma_dot( 5, 14, 1, 2 ) = -1.900300;
558  rsigma_dot( 5, 15, 1, 2 ) = -0.801700;
559  rsigma_dot( 5, 16, 1, 2 ) = -0.801700;
560  rsigma_dot( 5, 17, 1, 2 ) = -0.801700;
561  rsigma_dot( 5, 18, 1, 2 ) = -0.801700;
562  rsigma_dot( 6, 1, 1, 2 ) = -0.750400;
563  rsigma_dot( 6, 2, 1, 2 ) = -0.750400;
564  rsigma_dot( 6, 3, 1, 2 ) = -0.750400;
565  rsigma_dot( 6, 4, 1, 2 ) = -0.750400;
566  rsigma_dot( 6, 5, 1, 2 ) = -0.750400;
567  rsigma_dot( 6, 6, 1, 2 ) = -2.359800;
568  rsigma_dot( 6, 7, 1, 2 ) = -4.511580;
569  rsigma_dot( 6, 8, 1, 2 ) = -5.325090;
570  rsigma_dot( 6, 9, 1, 2 ) = -6.768980;
571  rsigma_dot( 6, 10, 1, 2 ) = -8.613260;
572  rsigma_dot( 6, 11, 1, 2 ) = -9.405770;
573  rsigma_dot( 6, 12, 1, 2 ) = -7.032650;
574  rsigma_dot( 6, 13, 1, 2 ) = -3.885880;
575  rsigma_dot( 6, 14, 1, 2 ) = -3.315300;
576  rsigma_dot( 6, 15, 1, 2 ) = -0.750400;
577  rsigma_dot( 6, 16, 1, 2 ) = -0.750400;
578  rsigma_dot( 6, 17, 1, 2 ) = -0.750400;
579  rsigma_dot( 6, 18, 1, 2 ) = -0.750400;
580  rsigma_dot( 7, 1, 1, 2 ) = -0.701600;
581  rsigma_dot( 7, 2, 1, 2 ) = -0.701600;
582  rsigma_dot( 7, 3, 1, 2 ) = -0.701600;
583  rsigma_dot( 7, 4, 1, 2 ) = -0.701600;
584  rsigma_dot( 7, 5, 1, 2 ) = -0.701600;
585  rsigma_dot( 7, 6, 1, 2 ) = -2.311000;
586  rsigma_dot( 7, 7, 1, 2 ) = -4.365150;
587  rsigma_dot( 7, 8, 1, 2 ) = -5.446520;
588  rsigma_dot( 7, 9, 1, 2 ) = -6.795160;
589  rsigma_dot( 7, 10, 1, 2 ) = -7.145720;
590  rsigma_dot( 7, 11, 1, 2 ) = -7.527050;
591  rsigma_dot( 7, 12, 1, 2 ) = -6.885740;
592  rsigma_dot( 7, 13, 1, 2 ) = -4.508250;
593  rsigma_dot( 7, 14, 1, 2 ) = -1.800200;
594  rsigma_dot( 7, 15, 1, 2 ) = -2.647500;
595  rsigma_dot( 7, 16, 1, 2 ) = -0.701600;
596  rsigma_dot( 7, 17, 1, 2 ) = -0.701600;
597  rsigma_dot( 7, 18, 1, 2 ) = -0.701600;
598  rsigma_dot( 8, 1, 1, 2 ) = -0.655100;
599  rsigma_dot( 8, 2, 1, 2 ) = -0.655100;
600  rsigma_dot( 8, 3, 1, 2 ) = -0.655100;
601  rsigma_dot( 8, 4, 1, 2 ) = -0.655100;
602  rsigma_dot( 8, 5, 1, 2 ) = -0.655100;
603  rsigma_dot( 8, 6, 1, 2 ) = -3.699600;
604  rsigma_dot( 8, 7, 1, 2 ) = -4.089060;
605  rsigma_dot( 8, 8, 1, 2 ) = -5.417250;
606  rsigma_dot( 8, 9, 1, 2 ) = -6.293430;
607  rsigma_dot( 8, 10, 1, 2 ) = -6.663890;
608  rsigma_dot( 8, 11, 1, 2 ) = -6.788470;
609  rsigma_dot( 8, 12, 1, 2 ) = -6.368800;
610  rsigma_dot( 8, 13, 1, 2 ) = -4.732610;
611  rsigma_dot( 8, 14, 1, 2 ) = -1.753700;
612  rsigma_dot( 8, 15, 1, 2 ) = -0.655100;
613  rsigma_dot( 8, 16, 1, 2 ) = -0.655100;
614  rsigma_dot( 8, 17, 1, 2 ) = -0.655100;
615  rsigma_dot( 8, 18, 1, 2 ) = -0.655100;
616  rsigma_dot( 9, 1, 1, 2 ) = -0.610600;
617  rsigma_dot( 9, 2, 1, 2 ) = -0.610600;
618  rsigma_dot( 9, 3, 1, 2 ) = -0.610600;
619  rsigma_dot( 9, 4, 1, 2 ) = -0.610600;
620  rsigma_dot( 9, 5, 1, 2 ) = -1.709200;
621  rsigma_dot( 9, 6, 1, 2 ) = -3.906460;
622  rsigma_dot( 9, 7, 1, 2 ) = -3.906460;
623  rsigma_dot( 9, 8, 1, 2 ) = -5.355550;
624  rsigma_dot( 9, 9, 1, 2 ) = -6.070210;
625  rsigma_dot( 9, 10, 1, 2 ) = -5.971910;
626  rsigma_dot( 9, 11, 1, 2 ) = -6.070210;
627  rsigma_dot( 9, 12, 1, 2 ) = -5.852370;
628  rsigma_dot( 9, 13, 1, 2 ) = -4.901080;
629  rsigma_dot( 9, 14, 1, 2 ) = -2.807800;
630  rsigma_dot( 9, 15, 1, 2 ) = -2.220100;
631  rsigma_dot( 9, 16, 1, 2 ) = -0.610600;
632  rsigma_dot( 9, 17, 1, 2 ) = -1.709200;
633  rsigma_dot( 9, 18, 1, 2 ) = -0.610600;
634  rsigma_dot( 10, 1, 1, 2 ) = -0.568100;
635  rsigma_dot( 10, 2, 1, 2 ) = -0.568100;
636  rsigma_dot( 10, 3, 1, 2 ) = -0.568100;
637  rsigma_dot( 10, 4, 1, 2 ) = -2.177500;
638  rsigma_dot( 10, 5, 1, 2 ) = -2.765300;
639  rsigma_dot( 10, 6, 1, 2 ) = -4.002050;
640  rsigma_dot( 10, 7, 1, 2 ) = -4.802170;
641  rsigma_dot( 10, 8, 1, 2 ) = -5.598500;
642  rsigma_dot( 10, 9, 1, 2 ) = -5.809810;
643  rsigma_dot( 10, 10, 1, 2 ) = -5.809810;
644  rsigma_dot( 10, 11, 1, 2 ) = -5.732850;
645  rsigma_dot( 10, 12, 1, 2 ) = -5.572010;
646  rsigma_dot( 10, 13, 1, 2 ) = -4.858520;
647  rsigma_dot( 10, 14, 1, 2 ) = -3.935360;
648  rsigma_dot( 10, 15, 1, 2 ) = -1.666700;
649  rsigma_dot( 10, 16, 1, 2 ) = -1.666700;
650  rsigma_dot( 10, 17, 1, 2 ) = -0.568100;
651  rsigma_dot( 10, 18, 1, 2 ) = -1.666700;
652  rsigma_dot( 11, 1, 1, 2 ) = -0.527200;
653  rsigma_dot( 11, 2, 1, 2 ) = -0.527200;
654  rsigma_dot( 11, 3, 1, 2 ) = -0.527200;
655  rsigma_dot( 11, 4, 1, 2 ) = -1.625900;
656  rsigma_dot( 11, 5, 1, 2 ) = -3.746110;
657  rsigma_dot( 11, 6, 1, 2 ) = -4.190800;
658  rsigma_dot( 11, 7, 1, 2 ) = -4.761350;
659  rsigma_dot( 11, 8, 1, 2 ) = -5.897880;
660  rsigma_dot( 11, 9, 1, 2 ) = -5.850250;
661  rsigma_dot( 11, 10, 1, 2 ) = -5.789930;
662  rsigma_dot( 11, 11, 1, 2 ) = -5.371430;
663  rsigma_dot( 11, 12, 1, 2 ) = -5.608640;
664  rsigma_dot( 11, 13, 1, 2 ) = -5.306360;
665  rsigma_dot( 11, 14, 1, 2 ) = -4.288440;
666  rsigma_dot( 11, 15, 1, 2 ) = -3.092200;
667  rsigma_dot( 11, 16, 1, 2 ) = -0.527200;
668  rsigma_dot( 11, 17, 1, 2 ) = -0.527200;
669  rsigma_dot( 11, 18, 1, 2 ) = -0.527200;
670  rsigma_dot( 12, 1, 1, 2 ) = -0.488000;
671  rsigma_dot( 12, 2, 1, 2 ) = -0.488000;
672  rsigma_dot( 12, 3, 1, 2 ) = -0.488000;
673  rsigma_dot( 12, 4, 1, 2 ) = -3.623500;
674  rsigma_dot( 12, 5, 1, 2 ) = -4.722120;
675  rsigma_dot( 12, 6, 1, 2 ) = -4.692710;
676  rsigma_dot( 12, 7, 1, 2 ) = -5.197550;
677  rsigma_dot( 12, 8, 1, 2 ) = -5.544260;
678  rsigma_dot( 12, 9, 1, 2 ) = -5.740290;
679  rsigma_dot( 12, 10, 1, 2 ) = -5.478450;
680  rsigma_dot( 12, 11, 1, 2 ) = -5.347830;
681  rsigma_dot( 12, 12, 1, 2 ) = -5.363220;
682  rsigma_dot( 12, 13, 1, 2 ) = -5.103140;
683  rsigma_dot( 12, 14, 1, 2 ) = -4.722120;
684  rsigma_dot( 12, 15, 1, 2 ) = -3.984530;
685  rsigma_dot( 12, 16, 1, 2 ) = -3.196100;
686  rsigma_dot( 12, 17, 1, 2 ) = -0.488000;
687  rsigma_dot( 12, 18, 1, 2 ) = -0.488000;
688 
689  // section 11
690  rsigma_dot( 1, 1, 1, 1 ) = -0.552300;
691  rsigma_dot( 1, 2, 1, 1 ) = -0.552300;
692  rsigma_dot( 1, 3, 1, 1 ) = -0.552300;
693  rsigma_dot( 1, 4, 1, 1 ) = -0.552300;
694  rsigma_dot( 1, 5, 1, 1 ) = -0.552300;
695  rsigma_dot( 1, 6, 1, 1 ) = -0.552300;
696  rsigma_dot( 1, 7, 1, 1 ) = -0.552300;
697  rsigma_dot( 1, 8, 1, 1 ) = -0.552300;
698  rsigma_dot( 1, 9, 1, 1 ) = -0.552300;
699  rsigma_dot( 1, 10, 1, 1 ) = -1.651000;
700  rsigma_dot( 1, 11, 1, 1 ) = -0.552300;
701  rsigma_dot( 1, 12, 1, 1 ) = -0.552300;
702  rsigma_dot( 1, 13, 1, 1 ) = -0.552300;
703  rsigma_dot( 1, 14, 1, 1 ) = -0.552300;
704  rsigma_dot( 1, 15, 1, 1 ) = -0.552300;
705  rsigma_dot( 1, 16, 1, 1 ) = -0.552300;
706  rsigma_dot( 1, 17, 1, 1 ) = -0.552300;
707  rsigma_dot( 1, 18, 1, 1 ) = -0.552300;
708  rsigma_dot( 2, 1, 1, 1 ) = -0.487800;
709  rsigma_dot( 2, 2, 1, 1 ) = -0.487800;
710  rsigma_dot( 2, 3, 1, 1 ) = -0.487800;
711  rsigma_dot( 2, 4, 1, 1 ) = -0.487800;
712  rsigma_dot( 2, 5, 1, 1 ) = -0.487800;
713  rsigma_dot( 2, 6, 1, 1 ) = -0.487800;
714  rsigma_dot( 2, 7, 1, 1 ) = -0.487800;
715  rsigma_dot( 2, 8, 1, 1 ) = -0.487800;
716  rsigma_dot( 2, 9, 1, 1 ) = -0.487800;
717  rsigma_dot( 2, 10, 1, 1 ) = -1.586400;
718  rsigma_dot( 2, 11, 1, 1 ) = -2.097200;
719  rsigma_dot( 2, 12, 1, 1 ) = -1.586400;
720  rsigma_dot( 2, 13, 1, 1 ) = -0.487800;
721  rsigma_dot( 2, 14, 1, 1 ) = -0.487800;
722  rsigma_dot( 2, 15, 1, 1 ) = -0.487800;
723  rsigma_dot( 2, 16, 1, 1 ) = -0.487800;
724  rsigma_dot( 2, 17, 1, 1 ) = -0.487800;
725  rsigma_dot( 2, 18, 1, 1 ) = -0.487800;
726  rsigma_dot( 3, 1, 1, 1 ) = -0.427200;
727  rsigma_dot( 3, 2, 1, 1 ) = -0.427200;
728  rsigma_dot( 3, 3, 1, 1 ) = -0.427200;
729  rsigma_dot( 3, 4, 1, 1 ) = -0.427200;
730  rsigma_dot( 3, 5, 1, 1 ) = -0.427200;
731  rsigma_dot( 3, 6, 1, 1 ) = -1.525800;
732  rsigma_dot( 3, 7, 1, 1 ) = -2.036600;
733  rsigma_dot( 3, 8, 1, 1 ) = -2.036600;
734  rsigma_dot( 3, 9, 1, 1 ) = -0.427200;
735  rsigma_dot( 3, 10, 1, 1 ) = -2.373100;
736  rsigma_dot( 3, 11, 1, 1 ) = -1.525800;
737  rsigma_dot( 3, 12, 1, 1 ) = -2.992100;
738  rsigma_dot( 3, 13, 1, 1 ) = -2.373100;
739  rsigma_dot( 3, 14, 1, 1 ) = -0.427200;
740  rsigma_dot( 3, 15, 1, 1 ) = -0.427200;
741  rsigma_dot( 3, 16, 1, 1 ) = -0.427200;
742  rsigma_dot( 3, 17, 1, 1 ) = -0.427200;
743  rsigma_dot( 3, 18, 1, 1 ) = -0.427200;
744  rsigma_dot( 4, 1, 1, 1 ) = -0.370000;
745  rsigma_dot( 4, 2, 1, 1 ) = -0.370000;
746  rsigma_dot( 4, 3, 1, 1 ) = -0.370000;
747  rsigma_dot( 4, 4, 1, 1 ) = -0.370000;
748  rsigma_dot( 4, 5, 1, 1 ) = -0.370000;
749  rsigma_dot( 4, 6, 1, 1 ) = -0.370000;
750  rsigma_dot( 4, 7, 1, 1 ) = -2.315900;
751  rsigma_dot( 4, 8, 1, 1 ) = -2.567200;
752  rsigma_dot( 4, 9, 1, 1 ) = -2.767900;
753  rsigma_dot( 4, 10, 1, 1 ) = -3.588900;
754  rsigma_dot( 4, 11, 1, 1 ) = -3.414500;
755  rsigma_dot( 4, 12, 1, 1 ) = -4.033580;
756  rsigma_dot( 4, 13, 1, 1 ) = -3.804010;
757  rsigma_dot( 4, 14, 1, 1 ) = -0.370000;
758  rsigma_dot( 4, 15, 1, 1 ) = -0.370000;
759  rsigma_dot( 4, 16, 1, 1 ) = -0.370000;
760  rsigma_dot( 4, 17, 1, 1 ) = -0.370000;
761  rsigma_dot( 4, 18, 1, 1 ) = -0.370000;
762  rsigma_dot( 5, 1, 1, 1 ) = -0.316000;
763  rsigma_dot( 5, 2, 1, 1 ) = -0.316000;
764  rsigma_dot( 5, 3, 1, 1 ) = -0.316000;
765  rsigma_dot( 5, 4, 1, 1 ) = -0.316000;
766  rsigma_dot( 5, 5, 1, 1 ) = -0.316000;
767  rsigma_dot( 5, 6, 1, 1 ) = -1.925400;
768  rsigma_dot( 5, 7, 1, 1 ) = -3.683200;
769  rsigma_dot( 5, 8, 1, 1 ) = -4.122610;
770  rsigma_dot( 5, 9, 1, 1 ) = -3.812460;
771  rsigma_dot( 5, 10, 1, 1 ) = -3.683200;
772  rsigma_dot( 5, 11, 1, 1 ) = -3.360500;
773  rsigma_dot( 5, 12, 1, 1 ) = -4.969910;
774  rsigma_dot( 5, 13, 1, 1 ) = -4.869830;
775  rsigma_dot( 5, 14, 1, 1 ) = -3.149200;
776  rsigma_dot( 5, 15, 1, 1 ) = -0.316000;
777  rsigma_dot( 5, 16, 1, 1 ) = -0.316000;
778  rsigma_dot( 5, 17, 1, 1 ) = -0.316000;
779  rsigma_dot( 5, 18, 1, 1 ) = -0.316000;
780  rsigma_dot( 6, 1, 1, 1 ) = -0.264700;
781  rsigma_dot( 6, 2, 1, 1 ) = -0.264700;
782  rsigma_dot( 6, 3, 1, 1 ) = -0.264700;
783  rsigma_dot( 6, 4, 1, 1 ) = -0.264700;
784  rsigma_dot( 6, 5, 1, 1 ) = -0.264700;
785  rsigma_dot( 6, 6, 1, 1 ) = -2.829600;
786  rsigma_dot( 6, 7, 1, 1 ) = -5.680760;
787  rsigma_dot( 6, 8, 1, 1 ) = -5.076840;
788  rsigma_dot( 6, 9, 1, 1 ) = -4.271990;
789  rsigma_dot( 6, 10, 1, 1 ) = -4.234950;
790  rsigma_dot( 6, 11, 1, 1 ) = -3.761170;
791  rsigma_dot( 6, 12, 1, 1 ) = -5.567960;
792  rsigma_dot( 6, 13, 1, 1 ) = -6.108200;
793  rsigma_dot( 6, 14, 1, 1 ) = -4.956010;
794  rsigma_dot( 6, 15, 1, 1 ) = -1.874100;
795  rsigma_dot( 6, 16, 1, 1 ) = -0.264700;
796  rsigma_dot( 6, 17, 1, 1 ) = -0.264700;
797  rsigma_dot( 6, 18, 1, 1 ) = -0.264700;
798  rsigma_dot( 7, 1, 1, 1 ) = -0.215900;
799  rsigma_dot( 7, 2, 1, 1 ) = -0.215900;
800  rsigma_dot( 7, 3, 1, 1 ) = -0.215900;
801  rsigma_dot( 7, 4, 1, 1 ) = -0.215900;
802  rsigma_dot( 7, 5, 1, 1 ) = -0.215900;
803  rsigma_dot( 7, 6, 1, 1 ) = -5.011660;
804  rsigma_dot( 7, 7, 1, 1 ) = -7.313420;
805  rsigma_dot( 7, 8, 1, 1 ) = -6.580620;
806  rsigma_dot( 7, 9, 1, 1 ) = -4.726730;
807  rsigma_dot( 7, 10, 1, 1 ) = -4.748470;
808  rsigma_dot( 7, 11, 1, 1 ) = -2.923900;
809  rsigma_dot( 7, 12, 1, 1 ) = -5.700670;
810  rsigma_dot( 7, 13, 1, 1 ) = -6.799280;
811  rsigma_dot( 7, 14, 1, 1 ) = -6.059410;
812  rsigma_dot( 7, 15, 1, 1 ) = -2.161800;
813  rsigma_dot( 7, 16, 1, 1 ) = -0.215900;
814  rsigma_dot( 7, 17, 1, 1 ) = -0.215900;
815  rsigma_dot( 7, 18, 1, 1 ) = -0.215900;
816  rsigma_dot( 8, 1, 1, 1 ) = -0.169300;
817  rsigma_dot( 8, 2, 1, 1 ) = -0.169300;
818  rsigma_dot( 8, 3, 1, 1 ) = -0.169300;
819  rsigma_dot( 8, 4, 1, 1 ) = -0.169300;
820  rsigma_dot( 8, 5, 1, 1 ) = -0.169300;
821  rsigma_dot( 8, 6, 1, 1 ) = -6.132930;
822  rsigma_dot( 8, 7, 1, 1 ) = -8.735140;
823  rsigma_dot( 8, 8, 1, 1 ) = -7.789560;
824  rsigma_dot( 8, 9, 1, 1 ) = -5.199790;
825  rsigma_dot( 8, 10, 1, 1 ) = -4.981530;
826  rsigma_dot( 8, 11, 1, 1 ) = -3.882920;
827  rsigma_dot( 8, 12, 1, 1 ) = -5.059700;
828  rsigma_dot( 8, 13, 1, 1 ) = -7.329420;
829  rsigma_dot( 8, 14, 1, 1 ) = -7.217740;
830  rsigma_dot( 8, 15, 1, 1 ) = -4.723230;
831  rsigma_dot( 8, 16, 1, 1 ) = -0.169300;
832  rsigma_dot( 8, 17, 1, 1 ) = -0.169300;
833  rsigma_dot( 8, 18, 1, 1 ) = -0.169300;
834  rsigma_dot( 9, 1, 1, 1 ) = -0.124900;
835  rsigma_dot( 9, 2, 1, 1 ) = -0.124900;
836  rsigma_dot( 9, 3, 1, 1 ) = -0.124900;
837  rsigma_dot( 9, 4, 1, 1 ) = -1.223500;
838  rsigma_dot( 9, 5, 1, 1 ) = -3.069300;
839  rsigma_dot( 9, 6, 1, 1 ) = -6.756900;
840  rsigma_dot( 9, 7, 1, 1 ) = -8.997660;
841  rsigma_dot( 9, 8, 1, 1 ) = -7.292710;
842  rsigma_dot( 9, 9, 1, 1 ) = -5.704630;
843  rsigma_dot( 9, 10, 1, 1 ) = -4.442390;
844  rsigma_dot( 9, 11, 1, 1 ) = -4.202430;
845  rsigma_dot( 9, 12, 1, 1 ) = -4.567550;
846  rsigma_dot( 9, 13, 1, 1 ) = -7.550850;
847  rsigma_dot( 9, 14, 1, 1 ) = -8.291970;
848  rsigma_dot( 9, 15, 1, 1 ) = -5.617960;
849  rsigma_dot( 9, 16, 1, 1 ) = -1.223500;
850  rsigma_dot( 9, 17, 1, 1 ) = -1.223500;
851  rsigma_dot( 9, 18, 1, 1 ) = -0.124900;
852  rsigma_dot( 10, 1, 1, 1 ) = -0.082300;
853  rsigma_dot( 10, 2, 1, 1 ) = -0.082300;
854  rsigma_dot( 10, 3, 1, 1 ) = -0.082300;
855  rsigma_dot( 10, 4, 1, 1 ) = -1.180900;
856  rsigma_dot( 10, 5, 1, 1 ) = -3.932490;
857  rsigma_dot( 10, 6, 1, 1 ) = -6.858840;
858  rsigma_dot( 10, 7, 1, 1 ) = -8.186740;
859  rsigma_dot( 10, 8, 1, 1 ) = -6.228670;
860  rsigma_dot( 10, 9, 1, 1 ) = -5.720690;
861  rsigma_dot( 10, 10, 1, 1 ) = -4.548250;
862  rsigma_dot( 10, 11, 1, 1 ) = -4.451790;
863  rsigma_dot( 10, 12, 1, 1 ) = -4.548250;
864  rsigma_dot( 10, 13, 1, 1 ) = -7.242410;
865  rsigma_dot( 10, 14, 1, 1 ) = -8.619920;
866  rsigma_dot( 10, 15, 1, 1 ) = -6.349540;
867  rsigma_dot( 10, 16, 1, 1 ) = -3.026800;
868  rsigma_dot( 10, 17, 1, 1 ) = -1.180900;
869  rsigma_dot( 10, 18, 1, 1 ) = -1.180900;
870  rsigma_dot( 11, 1, 1, 1 ) = -0.041500;
871  rsigma_dot( 11, 2, 1, 1 ) = -0.041500;
872  rsigma_dot( 11, 3, 1, 1 ) = -0.041500;
873  rsigma_dot( 11, 4, 1, 1 ) = -1.987400;
874  rsigma_dot( 11, 5, 1, 1 ) = -4.275620;
875  rsigma_dot( 11, 6, 1, 1 ) = -6.582550;
876  rsigma_dot( 11, 7, 1, 1 ) = -6.915710;
877  rsigma_dot( 11, 8, 1, 1 ) = -6.166200;
878  rsigma_dot( 11, 9, 1, 1 ) = -5.483930;
879  rsigma_dot( 11, 10, 1, 1 ) = -4.595390;
880  rsigma_dot( 11, 11, 1, 1 ) = -4.616230;
881  rsigma_dot( 11, 12, 1, 1 ) = -4.931860;
882  rsigma_dot( 11, 13, 1, 1 ) = -6.246070;
883  rsigma_dot( 11, 14, 1, 1 ) = -8.296560;
884  rsigma_dot( 11, 15, 1, 1 ) = -6.734840;
885  rsigma_dot( 11, 16, 1, 1 ) = -3.337400;
886  rsigma_dot( 11, 17, 1, 1 ) = -1.140100;
887  rsigma_dot( 11, 18, 1, 1 ) = -0.041500;
888  rsigma_dot( 12, 1, 1, 1 ) = -0.002300;
889  rsigma_dot( 12, 2, 1, 1 ) = -0.002300;
890  rsigma_dot( 12, 3, 1, 1 ) = -0.002300;
891  rsigma_dot( 12, 4, 1, 1 ) = -2.199500;
892  rsigma_dot( 12, 5, 1, 1 ) = -4.729680;
893  rsigma_dot( 12, 6, 1, 1 ) = -6.131350;
894  rsigma_dot( 12, 7, 1, 1 ) = -5.729140;
895  rsigma_dot( 12, 8, 1, 1 ) = -5.604410;
896  rsigma_dot( 12, 9, 1, 1 ) = -5.211780;
897  rsigma_dot( 12, 10, 1, 1 ) = -4.892640;
898  rsigma_dot( 12, 11, 1, 1 ) = -4.992730;
899  rsigma_dot( 12, 12, 1, 1 ) = -5.275290;
900  rsigma_dot( 12, 13, 1, 1 ) = -5.945090;
901  rsigma_dot( 12, 14, 1, 1 ) = -7.716970;
902  rsigma_dot( 12, 15, 1, 1 ) = -6.693140;
903  rsigma_dot( 12, 16, 1, 1 ) = -3.972590;
904  rsigma_dot( 12, 17, 1, 1 ) = -1.948200;
905  rsigma_dot( 12, 18, 1, 1 ) = -0.002300;
906 
907  // section 22
908  rsigma_dot( 1, 1, 2, 2 ) = -0.550000;
909  rsigma_dot( 1, 2, 2, 2 ) = -0.550000;
910  rsigma_dot( 1, 3, 2, 2 ) = -0.550000;
911  rsigma_dot( 1, 4, 2, 2 ) = -0.550000;
912  rsigma_dot( 1, 5, 2, 2 ) = -0.550000;
913  rsigma_dot( 1, 6, 2, 2 ) = -0.550000;
914  rsigma_dot( 1, 7, 2, 2 ) = -0.550000;
915  rsigma_dot( 1, 8, 2, 2 ) = -2.159400;
916  rsigma_dot( 1, 9, 2, 2 ) = -0.550000;
917  rsigma_dot( 1, 10, 2, 2 ) = -2.159400;
918  rsigma_dot( 1, 11, 2, 2 ) = -0.550000;
919  rsigma_dot( 1, 12, 2, 2 ) = -0.550000;
920  rsigma_dot( 1, 13, 2, 2 ) = -0.550000;
921  rsigma_dot( 1, 14, 2, 2 ) = -0.550000;
922  rsigma_dot( 1, 15, 2, 2 ) = -0.550000;
923  rsigma_dot( 1, 16, 2, 2 ) = -0.550000;
924  rsigma_dot( 1, 17, 2, 2 ) = -0.550000;
925  rsigma_dot( 1, 18, 2, 2 ) = -0.550000;
926  rsigma_dot( 2, 1, 2, 2 ) = -0.485500;
927  rsigma_dot( 2, 2, 2, 2 ) = -0.485500;
928  rsigma_dot( 2, 3, 2, 2 ) = -0.485500;
929  rsigma_dot( 2, 4, 2, 2 ) = -0.485500;
930  rsigma_dot( 2, 5, 2, 2 ) = -0.485500;
931  rsigma_dot( 2, 6, 2, 2 ) = -1.584100;
932  rsigma_dot( 2, 7, 2, 2 ) = -3.050400;
933  rsigma_dot( 2, 8, 2, 2 ) = -2.094900;
934  rsigma_dot( 2, 9, 2, 2 ) = -2.431400;
935  rsigma_dot( 2, 10, 2, 2 ) = -1.584100;
936  rsigma_dot( 2, 11, 2, 2 ) = -0.485500;
937  rsigma_dot( 2, 12, 2, 2 ) = -0.485500;
938  rsigma_dot( 2, 13, 2, 2 ) = -0.485500;
939  rsigma_dot( 2, 14, 2, 2 ) = -0.485500;
940  rsigma_dot( 2, 15, 2, 2 ) = -0.485500;
941  rsigma_dot( 2, 16, 2, 2 ) = -0.485500;
942  rsigma_dot( 2, 17, 2, 2 ) = -0.485500;
943  rsigma_dot( 2, 18, 2, 2 ) = -0.485500;
944  rsigma_dot( 3, 1, 2, 2 ) = -0.424800;
945  rsigma_dot( 3, 2, 2, 2 ) = -0.424800;
946  rsigma_dot( 3, 3, 2, 2 ) = -0.424800;
947  rsigma_dot( 3, 4, 2, 2 ) = -0.424800;
948  rsigma_dot( 3, 5, 2, 2 ) = -0.424800;
949  rsigma_dot( 3, 6, 2, 2 ) = -3.258000;
950  rsigma_dot( 3, 7, 2, 2 ) = -4.231500;
951  rsigma_dot( 3, 8, 2, 2 ) = -3.980180;
952  rsigma_dot( 3, 9, 2, 2 ) = -2.370700;
953  rsigma_dot( 3, 10, 2, 2 ) = -1.523400;
954  rsigma_dot( 3, 11, 2, 2 ) = -2.822700;
955  rsigma_dot( 3, 12, 2, 2 ) = -0.424800;
956  rsigma_dot( 3, 13, 2, 2 ) = -0.424800;
957  rsigma_dot( 3, 14, 2, 2 ) = -0.424800;
958  rsigma_dot( 3, 15, 2, 2 ) = -0.424800;
959  rsigma_dot( 3, 16, 2, 2 ) = -0.424800;
960  rsigma_dot( 3, 17, 2, 2 ) = -0.424800;
961  rsigma_dot( 3, 18, 2, 2 ) = -0.424800;
962  rsigma_dot( 4, 1, 2, 2 ) = -0.367700;
963  rsigma_dot( 4, 2, 2, 2 ) = -0.367700;
964  rsigma_dot( 4, 3, 2, 2 ) = -0.367700;
965  rsigma_dot( 4, 4, 2, 2 ) = -0.367700;
966  rsigma_dot( 4, 5, 2, 2 ) = -3.312100;
967  rsigma_dot( 4, 6, 2, 2 ) = -3.412200;
968  rsigma_dot( 4, 7, 2, 2 ) = -5.509340;
969  rsigma_dot( 4, 8, 2, 2 ) = -4.128880;
970  rsigma_dot( 4, 9, 2, 2 ) = -3.075700;
971  rsigma_dot( 4, 10, 2, 2 ) = -2.765600;
972  rsigma_dot( 4, 11, 2, 2 ) = -2.765600;
973  rsigma_dot( 4, 12, 2, 2 ) = -2.765600;
974  rsigma_dot( 4, 13, 2, 2 ) = -0.367700;
975  rsigma_dot( 4, 14, 2, 2 ) = -0.367700;
976  rsigma_dot( 4, 15, 2, 2 ) = -0.367700;
977  rsigma_dot( 4, 16, 2, 2 ) = -0.367700;
978  rsigma_dot( 4, 17, 2, 2 ) = -0.367700;
979  rsigma_dot( 4, 18, 2, 2 ) = -0.367700;
980  rsigma_dot( 5, 1, 2, 2 ) = -0.313600;
981  rsigma_dot( 5, 2, 2, 2 ) = -0.313600;
982  rsigma_dot( 5, 3, 2, 2 ) = -0.313600;
983  rsigma_dot( 5, 4, 2, 2 ) = -1.412200;
984  rsigma_dot( 5, 5, 2, 2 ) = -3.021700;
985  rsigma_dot( 5, 6, 2, 2 ) = -4.779520;
986  rsigma_dot( 5, 7, 2, 2 ) = -6.825360;
987  rsigma_dot( 5, 8, 2, 2 ) = -4.657420;
988  rsigma_dot( 5, 9, 2, 2 ) = -3.358100;
989  rsigma_dot( 5, 10, 2, 2 ) = -3.609400;
990  rsigma_dot( 5, 11, 2, 2 ) = -4.245440;
991  rsigma_dot( 5, 12, 2, 2 ) = -3.021700;
992  rsigma_dot( 5, 13, 2, 2 ) = -1.412200;
993  rsigma_dot( 5, 14, 2, 2 ) = -0.313600;
994  rsigma_dot( 5, 15, 2, 2 ) = -0.313600;
995  rsigma_dot( 5, 16, 2, 2 ) = -0.313600;
996  rsigma_dot( 5, 17, 2, 2 ) = -0.313600;
997  rsigma_dot( 5, 18, 2, 2 ) = -0.313600;
998  rsigma_dot( 6, 1, 2, 2 ) = -0.262300;
999  rsigma_dot( 6, 2, 2, 2 ) = -0.262300;
1000  rsigma_dot( 6, 3, 2, 2 ) = -0.262300;
1001  rsigma_dot( 6, 4, 2, 2 ) = -1.871800;
1002  rsigma_dot( 6, 5, 2, 2 ) = -2.208200;
1003  rsigma_dot( 6, 6, 2, 2 ) = -6.862190;
1004  rsigma_dot( 6, 7, 2, 2 ) = -7.463490;
1005  rsigma_dot( 6, 8, 2, 2 ) = -4.837030;
1006  rsigma_dot( 6, 9, 2, 2 ) = -4.436700;
1007  rsigma_dot( 6, 10, 2, 2 ) = -4.916280;
1008  rsigma_dot( 6, 11, 2, 2 ) = -5.687270;
1009  rsigma_dot( 6, 12, 2, 2 ) = -4.405450;
1010  rsigma_dot( 6, 13, 2, 2 ) = -2.208200;
1011  rsigma_dot( 6, 14, 2, 2 ) = -0.262300;
1012  rsigma_dot( 6, 15, 2, 2 ) = -1.360900;
1013  rsigma_dot( 6, 16, 2, 2 ) = -0.262300;
1014  rsigma_dot( 6, 17, 2, 2 ) = -0.262300;
1015  rsigma_dot( 6, 18, 2, 2 ) = -0.262300;
1016  rsigma_dot( 7, 1, 2, 2 ) = -0.213500;
1017  rsigma_dot( 7, 2, 2, 2 ) = -0.213500;
1018  rsigma_dot( 7, 3, 2, 2 ) = -0.213500;
1019  rsigma_dot( 7, 4, 2, 2 ) = -0.213500;
1020  rsigma_dot( 7, 5, 2, 2 ) = -4.256580;
1021  rsigma_dot( 7, 6, 2, 2 ) = -8.696750;
1022  rsigma_dot( 7, 7, 2, 2 ) = -7.275720;
1023  rsigma_dot( 7, 8, 2, 2 ) = -4.746130;
1024  rsigma_dot( 7, 9, 2, 2 ) = -4.145350;
1025  rsigma_dot( 7, 10, 2, 2 ) = -4.923060;
1026  rsigma_dot( 7, 11, 2, 2 ) = -6.725270;
1027  rsigma_dot( 7, 12, 2, 2 ) = -5.620700;
1028  rsigma_dot( 7, 13, 2, 2 ) = -2.611400;
1029  rsigma_dot( 7, 14, 2, 2 ) = -0.213500;
1030  rsigma_dot( 7, 15, 2, 2 ) = -0.213500;
1031  rsigma_dot( 7, 16, 2, 2 ) = -0.213500;
1032  rsigma_dot( 7, 17, 2, 2 ) = -0.213500;
1033  rsigma_dot( 7, 18, 2, 2 ) = -0.213500;
1034  rsigma_dot( 8, 1, 2, 2 ) = -0.167000;
1035  rsigma_dot( 8, 2, 2, 2 ) = -0.167000;
1036  rsigma_dot( 8, 3, 2, 2 ) = -0.167000;
1037  rsigma_dot( 8, 4, 2, 2 ) = -1.265600;
1038  rsigma_dot( 8, 5, 2, 2 ) = -6.383610;
1039  rsigma_dot( 8, 6, 2, 2 ) = -9.098560;
1040  rsigma_dot( 8, 7, 2, 2 ) = -6.418910;
1041  rsigma_dot( 8, 8, 2, 2 ) = -4.561460;
1042  rsigma_dot( 8, 9, 2, 2 ) = -4.174340;
1043  rsigma_dot( 8, 10, 2, 2 ) = -4.979190;
1044  rsigma_dot( 8, 11, 2, 2 ) = -7.462060;
1045  rsigma_dot( 8, 12, 2, 2 ) = -7.337890;
1046  rsigma_dot( 8, 13, 2, 2 ) = -3.722400;
1047  rsigma_dot( 8, 14, 2, 2 ) = -0.167000;
1048  rsigma_dot( 8, 15, 2, 2 ) = -0.167000;
1049  rsigma_dot( 8, 16, 2, 2 ) = -0.167000;
1050  rsigma_dot( 8, 17, 2, 2 ) = -0.167000;
1051  rsigma_dot( 8, 18, 2, 2 ) = -0.167000;
1052  rsigma_dot( 9, 1, 2, 2 ) = -0.122600;
1053  rsigma_dot( 9, 2, 2, 2 ) = -0.122600;
1054  rsigma_dot( 9, 3, 2, 2 ) = -0.122600;
1055  rsigma_dot( 9, 4, 2, 2 ) = -1.221200;
1056  rsigma_dot( 9, 5, 2, 2 ) = -7.827370;
1057  rsigma_dot( 9, 6, 2, 2 ) = -8.067400;
1058  rsigma_dot( 9, 7, 2, 2 ) = -5.809530;
1059  rsigma_dot( 9, 8, 2, 2 ) = -4.966740;
1060  rsigma_dot( 9, 9, 2, 2 ) = -4.327250;
1061  rsigma_dot( 9, 10, 2, 2 ) = -4.884730;
1062  rsigma_dot( 9, 11, 2, 2 ) = -7.072410;
1063  rsigma_dot( 9, 12, 2, 2 ) = -8.543900;
1064  rsigma_dot( 9, 13, 2, 2 ) = -5.252450;
1065  rsigma_dot( 9, 14, 2, 2 ) = -1.732000;
1066  rsigma_dot( 9, 15, 2, 2 ) = -0.122600;
1067  rsigma_dot( 9, 16, 2, 2 ) = -0.122600;
1068  rsigma_dot( 9, 17, 2, 2 ) = -0.122600;
1069  rsigma_dot( 9, 18, 2, 2 ) = -0.122600;
1070  rsigma_dot( 10, 1, 2, 2 ) = -0.080000;
1071  rsigma_dot( 10, 2, 2, 2 ) = -0.080000;
1072  rsigma_dot( 10, 3, 2, 2 ) = -1.689400;
1073  rsigma_dot( 10, 4, 2, 2 ) = -2.913200;
1074  rsigma_dot( 10, 5, 2, 2 ) = -7.579420;
1075  rsigma_dot( 10, 6, 2, 2 ) = -6.907620;
1076  rsigma_dot( 10, 7, 2, 2 ) = -5.422330;
1077  rsigma_dot( 10, 8, 2, 2 ) = -5.161400;
1078  rsigma_dot( 10, 9, 2, 2 ) = -4.590850;
1079  rsigma_dot( 10, 10, 2, 2 ) = -4.970340;
1080  rsigma_dot( 10, 11, 2, 2 ) = -6.331900;
1081  rsigma_dot( 10, 12, 2, 2 ) = -8.504850;
1082  rsigma_dot( 10, 13, 2, 2 ) = -6.881280;
1083  rsigma_dot( 10, 14, 2, 2 ) = -2.644900;
1084  rsigma_dot( 10, 15, 2, 2 ) = -0.080000;
1085  rsigma_dot( 10, 16, 2, 2 ) = -1.178600;
1086  rsigma_dot( 10, 17, 2, 2 ) = -0.080000;
1087  rsigma_dot( 10, 18, 2, 2 ) = -0.080000;
1088  rsigma_dot( 11, 1, 2, 2 ) = -0.039200;
1089  rsigma_dot( 11, 2, 2, 2 ) = -1.137800;
1090  rsigma_dot( 11, 3, 2, 2 ) = -2.437100;
1091  rsigma_dot( 11, 4, 2, 2 ) = -4.356660;
1092  rsigma_dot( 11, 5, 2, 2 ) = -6.159470;
1093  rsigma_dot( 11, 6, 2, 2 ) = -6.109910;
1094  rsigma_dot( 11, 7, 2, 2 ) = -5.082600;
1095  rsigma_dot( 11, 8, 2, 2 ) = -5.203960;
1096  rsigma_dot( 11, 9, 2, 2 ) = -5.029610;
1097  rsigma_dot( 11, 10, 2, 2 ) = -5.180840;
1098  rsigma_dot( 11, 11, 2, 2 ) = -5.633880;
1099  rsigma_dot( 11, 12, 2, 2 ) = -7.850740;
1100  rsigma_dot( 11, 13, 2, 2 ) = -8.236990;
1101  rsigma_dot( 11, 14, 2, 2 ) = -3.258000;
1102  rsigma_dot( 11, 15, 2, 2 ) = -2.604100;
1103  rsigma_dot( 11, 16, 2, 2 ) = -1.648600;
1104  rsigma_dot( 11, 17, 2, 2 ) = -0.039200;
1105  rsigma_dot( 11, 18, 2, 2 ) = -0.039200;
1106  rsigma_dot( 12, 1, 2, 2 ) = 0.000000;
1107  rsigma_dot( 12, 2, 2, 2 ) = -2.197200;
1108  rsigma_dot( 12, 3, 2, 2 ) = -2.708000;
1109  rsigma_dot( 12, 4, 2, 2 ) = -5.093700;
1110  rsigma_dot( 12, 5, 2, 2 ) = -5.802070;
1111  rsigma_dot( 12, 6, 2, 2 ) = -5.587200;
1112  rsigma_dot( 12, 7, 2, 2 ) = -5.198450;
1113  rsigma_dot( 12, 8, 2, 2 ) = -5.459540;
1114  rsigma_dot( 12, 9, 2, 2 ) = -5.068860;
1115  rsigma_dot( 12, 10, 2, 2 ) = -5.272950;
1116  rsigma_dot( 12, 11, 2, 2 ) = -5.370590;
1117  rsigma_dot( 12, 12, 2, 2 ) = -6.597100;
1118  rsigma_dot( 12, 13, 2, 2 ) = -7.631870;
1119  rsigma_dot( 12, 14, 2, 2 ) = -5.043380;
1120  rsigma_dot( 12, 15, 2, 2 ) = -3.218800;
1121  rsigma_dot( 12, 16, 2, 2 ) = -2.833200;
1122  rsigma_dot( 12, 17, 2, 2 ) = -2.197200;
1123  rsigma_dot( 12, 18, 2, 2 ) = 0.000000;
1124 
1125  // section 21
1126  rsigma_dot( 1, 1, 2, 1 ) = -0.974100;
1127  rsigma_dot( 1, 2, 2, 1 ) = -0.974100;
1128  rsigma_dot( 1, 3, 2, 1 ) = -0.974100;
1129  rsigma_dot( 1, 4, 2, 1 ) = -0.974100;
1130  rsigma_dot( 1, 5, 2, 1 ) = -0.974100;
1131  rsigma_dot( 1, 6, 2, 1 ) = -0.974100;
1132  rsigma_dot( 1, 7, 2, 1 ) = -0.974100;
1133  rsigma_dot( 1, 8, 2, 1 ) = -2.920000;
1134  rsigma_dot( 1, 9, 2, 1 ) = -3.171400;
1135  rsigma_dot( 1, 10, 2, 1 ) = -0.974100;
1136  rsigma_dot( 1, 11, 2, 1 ) = -2.583600;
1137  rsigma_dot( 1, 12, 2, 1 ) = -0.974100;
1138  rsigma_dot( 1, 13, 2, 1 ) = -0.974100;
1139  rsigma_dot( 1, 14, 2, 1 ) = -0.974100;
1140  rsigma_dot( 1, 15, 2, 1 ) = -0.974100;
1141  rsigma_dot( 1, 16, 2, 1 ) = -0.974100;
1142  rsigma_dot( 1, 17, 2, 1 ) = -0.974100;
1143  rsigma_dot( 1, 18, 2, 1 ) = -0.974100;
1144  rsigma_dot( 2, 1, 2, 1 ) = -0.909600;
1145  rsigma_dot( 2, 2, 2, 1 ) = -0.909600;
1146  rsigma_dot( 2, 3, 2, 1 ) = -0.909600;
1147  rsigma_dot( 2, 4, 2, 1 ) = -0.909600;
1148  rsigma_dot( 2, 5, 2, 1 ) = -0.909600;
1149  rsigma_dot( 2, 6, 2, 1 ) = -0.909600;
1150  rsigma_dot( 2, 7, 2, 1 ) = -0.909600;
1151  rsigma_dot( 2, 8, 2, 1 ) = -4.343580;
1152  rsigma_dot( 2, 9, 2, 1 ) = -6.015540;
1153  rsigma_dot( 2, 10, 2, 1 ) = -4.205430;
1154  rsigma_dot( 2, 11, 2, 1 ) = -2.008200;
1155  rsigma_dot( 2, 12, 2, 1 ) = -2.519000;
1156  rsigma_dot( 2, 13, 2, 1 ) = -0.909600;
1157  rsigma_dot( 2, 14, 2, 1 ) = -0.909600;
1158  rsigma_dot( 2, 15, 2, 1 ) = -0.909600;
1159  rsigma_dot( 2, 16, 2, 1 ) = -0.909600;
1160  rsigma_dot( 2, 17, 2, 1 ) = -0.909600;
1161  rsigma_dot( 2, 18, 2, 1 ) = -0.909600;
1162  rsigma_dot( 3, 1, 2, 1 ) = -0.849000;
1163  rsigma_dot( 3, 2, 2, 1 ) = -0.849000;
1164  rsigma_dot( 3, 3, 2, 1 ) = -0.849000;
1165  rsigma_dot( 3, 4, 2, 1 ) = -0.849000;
1166  rsigma_dot( 3, 5, 2, 1 ) = -0.849000;
1167  rsigma_dot( 3, 6, 2, 1 ) = -0.849000;
1168  rsigma_dot( 3, 7, 2, 1 ) = -3.046200;
1169  rsigma_dot( 3, 8, 2, 1 ) = -5.990630;
1170  rsigma_dot( 3, 9, 2, 1 ) = -7.627760;
1171  rsigma_dot( 3, 10, 2, 1 ) = -6.013760;
1172  rsigma_dot( 3, 11, 2, 1 ) = -2.794900;
1173  rsigma_dot( 3, 12, 2, 1 ) = -2.458400;
1174  rsigma_dot( 3, 13, 2, 1 ) = -0.849000;
1175  rsigma_dot( 3, 14, 2, 1 ) = -0.849000;
1176  rsigma_dot( 3, 15, 2, 1 ) = -0.849000;
1177  rsigma_dot( 3, 16, 2, 1 ) = -0.849000;
1178  rsigma_dot( 3, 17, 2, 1 ) = -0.849000;
1179  rsigma_dot( 3, 18, 2, 1 ) = -0.849000;
1180  rsigma_dot( 4, 1, 2, 1 ) = -0.791800;
1181  rsigma_dot( 4, 2, 2, 1 ) = -0.791800;
1182  rsigma_dot( 4, 3, 2, 1 ) = -0.791800;
1183  rsigma_dot( 4, 4, 2, 1 ) = -0.791800;
1184  rsigma_dot( 4, 5, 2, 1 ) = -0.791800;
1185  rsigma_dot( 4, 6, 2, 1 ) = -2.401300;
1186  rsigma_dot( 4, 7, 2, 1 ) = -3.736250;
1187  rsigma_dot( 4, 8, 2, 1 ) = -6.512120;
1188  rsigma_dot( 4, 9, 2, 1 ) = -8.856760;
1189  rsigma_dot( 4, 10, 2, 1 ) = -7.235940;
1190  rsigma_dot( 4, 11, 2, 1 ) = -4.010690;
1191  rsigma_dot( 4, 12, 2, 1 ) = -2.401300;
1192  rsigma_dot( 4, 13, 2, 1 ) = -1.890400;
1193  rsigma_dot( 4, 14, 2, 1 ) = -0.791800;
1194  rsigma_dot( 4, 15, 2, 1 ) = -0.791800;
1195  rsigma_dot( 4, 16, 2, 1 ) = -0.791800;
1196  rsigma_dot( 4, 17, 2, 1 ) = -0.791800;
1197  rsigma_dot( 4, 18, 2, 1 ) = -0.791800;
1198  rsigma_dot( 5, 1, 2, 1 ) = -0.737700;
1199  rsigma_dot( 5, 2, 2, 1 ) = -0.737700;
1200  rsigma_dot( 5, 3, 2, 1 ) = -0.737700;
1201  rsigma_dot( 5, 4, 2, 1 ) = -0.737700;
1202  rsigma_dot( 5, 5, 2, 1 ) = -0.737700;
1203  rsigma_dot( 5, 6, 2, 1 ) = -2.683700;
1204  rsigma_dot( 5, 7, 2, 1 ) = -4.171730;
1205  rsigma_dot( 5, 8, 2, 1 ) = -7.436010;
1206  rsigma_dot( 5, 9, 2, 1 ) = -9.773380;
1207  rsigma_dot( 5, 10, 2, 1 ) = -8.510500;
1208  rsigma_dot( 5, 11, 2, 1 ) = -4.669570;
1209  rsigma_dot( 5, 12, 2, 1 ) = -3.445800;
1210  rsigma_dot( 5, 13, 2, 1 ) = -1.836400;
1211  rsigma_dot( 5, 14, 2, 1 ) = -0.737700;
1212  rsigma_dot( 5, 15, 2, 1 ) = -0.737700;
1213  rsigma_dot( 5, 16, 2, 1 ) = -0.737700;
1214  rsigma_dot( 5, 17, 2, 1 ) = -0.737700;
1215  rsigma_dot( 5, 18, 2, 1 ) = -0.737700;
1216  rsigma_dot( 6, 1, 2, 1 ) = -0.686500;
1217  rsigma_dot( 6, 2, 2, 1 ) = -0.686500;
1218  rsigma_dot( 6, 3, 2, 1 ) = -0.686500;
1219  rsigma_dot( 6, 4, 2, 1 ) = -0.686500;
1220  rsigma_dot( 6, 5, 2, 1 ) = -0.686500;
1221  rsigma_dot( 6, 6, 2, 1 ) = -0.686500;
1222  rsigma_dot( 6, 7, 2, 1 ) = -4.297370;
1223  rsigma_dot( 6, 8, 2, 1 ) = -7.670240;
1224  rsigma_dot( 6, 9, 2, 1 ) = -9.274480;
1225  rsigma_dot( 6, 10, 2, 1 ) = -8.780530;
1226  rsigma_dot( 6, 11, 2, 1 ) = -5.395980;
1227  rsigma_dot( 6, 12, 2, 1 ) = -2.632400;
1228  rsigma_dot( 6, 13, 2, 1 ) = -3.084300;
1229  rsigma_dot( 6, 14, 2, 1 ) = -0.686500;
1230  rsigma_dot( 6, 15, 2, 1 ) = -0.686500;
1231  rsigma_dot( 6, 16, 2, 1 ) = -0.686500;
1232  rsigma_dot( 6, 17, 2, 1 ) = -0.686500;
1233  rsigma_dot( 6, 18, 2, 1 ) = -0.686500;
1234  rsigma_dot( 7, 1, 2, 1 ) = -0.637700;
1235  rsigma_dot( 7, 2, 2, 1 ) = -0.637700;
1236  rsigma_dot( 7, 3, 2, 1 ) = -0.637700;
1237  rsigma_dot( 7, 4, 2, 1 ) = -0.637700;
1238  rsigma_dot( 7, 5, 2, 1 ) = -1.736300;
1239  rsigma_dot( 7, 6, 2, 1 ) = -3.682200;
1240  rsigma_dot( 7, 7, 2, 1 ) = -4.645000;
1241  rsigma_dot( 7, 8, 2, 1 ) = -6.967380;
1242  rsigma_dot( 7, 9, 2, 1 ) = -7.766160;
1243  rsigma_dot( 7, 10, 2, 1 ) = -7.608390;
1244  rsigma_dot( 7, 11, 2, 1 ) = -5.416790;
1245  rsigma_dot( 7, 12, 2, 1 ) = -3.773160;
1246  rsigma_dot( 7, 13, 2, 1 ) = -3.345700;
1247  rsigma_dot( 7, 14, 2, 1 ) = -2.834900;
1248  rsigma_dot( 7, 15, 2, 1 ) = -2.247100;
1249  rsigma_dot( 7, 16, 2, 1 ) = -0.637700;
1250  rsigma_dot( 7, 17, 2, 1 ) = -0.637700;
1251  rsigma_dot( 7, 18, 2, 1 ) = -0.637700;
1252  rsigma_dot( 8, 1, 2, 1 ) = -0.591100;
1253  rsigma_dot( 8, 2, 2, 1 ) = -0.591100;
1254  rsigma_dot( 8, 3, 2, 1 ) = -0.591100;
1255  rsigma_dot( 8, 4, 2, 1 ) = -0.591100;
1256  rsigma_dot( 8, 5, 2, 1 ) = -2.788400;
1257  rsigma_dot( 8, 6, 2, 1 ) = -4.087650;
1258  rsigma_dot( 8, 7, 2, 1 ) = -4.702020;
1259  rsigma_dot( 8, 8, 2, 1 ) = -6.084200;
1260  rsigma_dot( 8, 9, 2, 1 ) = -6.877140;
1261  rsigma_dot( 8, 10, 2, 1 ) = -6.827510;
1262  rsigma_dot( 8, 11, 2, 1 ) = -5.980210;
1263  rsigma_dot( 8, 12, 2, 1 ) = -4.087650;
1264  rsigma_dot( 8, 13, 2, 1 ) = -3.635700;
1265  rsigma_dot( 8, 14, 2, 1 ) = -4.202060;
1266  rsigma_dot( 8, 15, 2, 1 ) = -2.537100;
1267  rsigma_dot( 8, 16, 2, 1 ) = -0.591100;
1268  rsigma_dot( 8, 17, 2, 1 ) = -0.591100;
1269  rsigma_dot( 8, 18, 2, 1 ) = -0.591100;
1270  rsigma_dot( 9, 1, 2, 1 ) = -0.546700;
1271  rsigma_dot( 9, 2, 2, 1 ) = -0.546700;
1272  rsigma_dot( 9, 3, 2, 1 ) = -0.546700;
1273  rsigma_dot( 9, 4, 2, 1 ) = -0.546700;
1274  rsigma_dot( 9, 5, 2, 1 ) = -3.254700;
1275  rsigma_dot( 9, 6, 2, 1 ) = -4.102040;
1276  rsigma_dot( 9, 7, 2, 1 ) = -4.916140;
1277  rsigma_dot( 9, 8, 2, 1 ) = -5.788440;
1278  rsigma_dot( 9, 9, 2, 1 ) = -5.777800;
1279  rsigma_dot( 9, 10, 2, 1 ) = -6.126420;
1280  rsigma_dot( 9, 11, 2, 1 ) = -5.997730;
1281  rsigma_dot( 9, 12, 2, 1 ) = -4.721080;
1282  rsigma_dot( 9, 13, 2, 1 ) = -4.102040;
1283  rsigma_dot( 9, 14, 2, 1 ) = -4.157610;
1284  rsigma_dot( 9, 15, 2, 1 ) = -3.111600;
1285  rsigma_dot( 9, 16, 2, 1 ) = -1.645300;
1286  rsigma_dot( 9, 17, 2, 1 ) = -0.546700;
1287  rsigma_dot( 9, 18, 2, 1 ) = -0.546700;
1288  rsigma_dot( 10, 1, 2, 1 ) = -0.504100;
1289  rsigma_dot( 10, 2, 2, 1 ) = -0.504100;
1290  rsigma_dot( 10, 3, 2, 1 ) = -0.504100;
1291  rsigma_dot( 10, 4, 2, 1 ) = -0.504100;
1292  rsigma_dot( 10, 5, 2, 1 ) = -3.548700;
1293  rsigma_dot( 10, 6, 2, 1 ) = -4.873580;
1294  rsigma_dot( 10, 7, 2, 1 ) = -5.176960;
1295  rsigma_dot( 10, 8, 2, 1 ) = -5.756400;
1296  rsigma_dot( 10, 9, 2, 1 ) = -5.573030;
1297  rsigma_dot( 10, 10, 2, 1 ) = -5.702630;
1298  rsigma_dot( 10, 11, 2, 1 ) = -5.668920;
1299  rsigma_dot( 10, 12, 2, 1 ) = -5.078840;
1300  rsigma_dot( 10, 13, 2, 1 ) = -4.678520;
1301  rsigma_dot( 10, 14, 2, 1 ) = -4.310790;
1302  rsigma_dot( 10, 15, 2, 1 ) = -4.217700;
1303  rsigma_dot( 10, 16, 2, 1 ) = -2.113600;
1304  rsigma_dot( 10, 17, 2, 1 ) = -0.504100;
1305  rsigma_dot( 10, 18, 2, 1 ) = -0.504100;
1306  rsigma_dot( 11, 1, 2, 1 ) = -0.463300;
1307  rsigma_dot( 11, 2, 2, 1 ) = -0.463300;
1308  rsigma_dot( 11, 3, 2, 1 ) = -0.463300;
1309  rsigma_dot( 11, 4, 2, 1 ) = -2.409200;
1310  rsigma_dot( 11, 5, 2, 1 ) = -4.313460;
1311  rsigma_dot( 11, 6, 2, 1 ) = -5.225480;
1312  rsigma_dot( 11, 7, 2, 1 ) = -5.493750;
1313  rsigma_dot( 11, 8, 2, 1 ) = -5.353660;
1314  rsigma_dot( 11, 9, 2, 1 ) = -5.154660;
1315  rsigma_dot( 11, 10, 2, 1 ) = -5.397780;
1316  rsigma_dot( 11, 11, 2, 1 ) = -5.726000;
1317  rsigma_dot( 11, 12, 2, 1 ) = -5.480590;
1318  rsigma_dot( 11, 13, 2, 1 ) = -4.018660;
1319  rsigma_dot( 11, 14, 2, 1 ) = -4.313460;
1320  rsigma_dot( 11, 15, 2, 1 ) = -4.574180;
1321  rsigma_dot( 11, 16, 2, 1 ) = -2.861200;
1322  rsigma_dot( 11, 17, 2, 1 ) = -0.463300;
1323  rsigma_dot( 11, 18, 2, 1 ) = -1.561900;
1324  rsigma_dot( 12, 1, 2, 1 ) = -0.424100;
1325  rsigma_dot( 12, 2, 2, 1 ) = -0.424100;
1326  rsigma_dot( 12, 3, 2, 1 ) = -2.822000;
1327  rsigma_dot( 12, 4, 2, 1 ) = -3.132100;
1328  rsigma_dot( 12, 5, 2, 1 ) = -5.115440;
1329  rsigma_dot( 12, 6, 2, 1 ) = -5.697090;
1330  rsigma_dot( 12, 7, 2, 1 ) = -5.747100;
1331  rsigma_dot( 12, 8, 2, 1 ) = -5.386930;
1332  rsigma_dot( 12, 9, 2, 1 ) = -5.019210;
1333  rsigma_dot( 12, 10, 2, 1 ) = -5.186260;
1334  rsigma_dot( 12, 11, 2, 1 ) = -5.775950;
1335  rsigma_dot( 12, 12, 2, 1 ) = -5.600240;
1336  rsigma_dot( 12, 13, 2, 1 ) = -4.394380;
1337  rsigma_dot( 12, 14, 2, 1 ) = -3.257300;
1338  rsigma_dot( 12, 15, 2, 1 ) = -3.920600;
1339  rsigma_dot( 12, 16, 2, 1 ) = -3.643000;
1340  rsigma_dot( 12, 17, 2, 1 ) = -1.522700;
1341  rsigma_dot( 12, 18, 2, 1 ) = -0.424100;
1342 }
1343 
1344 
1345 } // ns sspot
1346 } // ns potentials
1347 } // ns fldsgn
1348 } // ns protocols
Generated on Sat Jun 1 2013 11:50:06 for Rosetta 3.5 by   doxygen 1.8.4