src/IntRank2/SolidSys.cc
// sosiz - 1.) still need to deal with when is coord vector 0 and when are the
//             coordinate flags necessary?
//         2.) need to see if the linked list of interactions is working here
1.) still need to deal with when is coord vector 0 and when are the
coordinate flags necessary.
2.) need to see if the linked list of interactions is working here
1 . ) still need to deal with when is coord vector 0 and when are the coordinate flags necessary . 2 . ) need to see if the linked list of interactions is working here
void solid_sys::setCoord(int i, coord& pt)
  {
  string I, S;
  if(i<0 || i>spins()) 			
    ssys_fatal(105, Gdec(i));
  if(!SCoords.size()) 			
    ssys_fatal(107, Gdec(i));
  SCoords.put(pt, i);			
  int j,k,l,ns=spins();			
  for(j=0,l=0; j<ns-1; j++)		
    {
    I = symbol(j);
    for(k=j+1; k<ns; k++,l++)		
      {
      S = symbol(k);
      if(k==i || j==i)			
        {				
        if(Dvec.DCC(l))			
          Dvec(l) = IntDip(I,S,SCoords(j),SCoords(k));
        }
      }
    }
  }
void solid_sys :: setCoord ( int i , coord & pt ) { string I , S ; if ( i < 0 || i > spins ( ) ) ssys_fatal ( 105 , Gdec ( i ) ) ; if ( ! SCoords . size ( ) ) ssys_fatal ( 107 , Gdec ( i ) ) ; SCoords . put ( pt , i ) ; int j , k , l , ns = spins ( ) ; for ( j = 0 , l = 0 ; j < ns - 1 ; j ++ ) { I = symbol ( j ) ; for ( k = j + 1 ; k < ns ; k ++ , l ++ ) { S = symbol ( k ) ; if ( k == i || j == i ) { if ( Dvec . DCC ( l ) ) Dvec ( l ) = IntDip ( I , S , SCoords ( j ) , SCoords ( k ) ) ; } } } }
sosiz - 1.)
sosiz - 1 . )
// Input                ssys    : Solid spin system //                      i       : spin index  //			pt	: spin coordinate // Output               none    : Spin i coordinate set to pt // Note				: Coordinate can only be set if there //				  already exists a complete set! // Insure valid spin index // Insure coordinates exist // Set the spin coordinate // Get the number of spins // Loop over all spin pairs // and see which involve the // new spin.  Check if the spin // is involved in any dipolar // interaction and, if so,