src/IntRank2/SolidSys.cc
/* For this spin system we (should) know the number of spins and their isotope
   types.  The dipolar interactions are all to be placed into the dipolar
   interaction vector, Dvec.  This must be read in using spin (NOT interaction)
   indices.  For example, { Iso(i), coord(i), Iso(j), coord(j) } for all spin
   pairs.  However, not everyone likes to use spin coordinates for designating
   a dipolar interaction, so we must allow for their setting without the use
   of coordinates, but using { Isotopes, spin indices, + other info }

	   Input		ssys	: Solid spin system (this)
	   			pset	: A parameter set
	  			ccount  : Count of specified coordinates
	   Output		none	: System dipolar interactions are
	  			  	  set from parameters in pset
	   Note				: Assumes that
					    1.) System isotopes are set
	  				    2.) Any spin coordinates have
	  					been set (setCoords)
	  				    3.) All dipolar interactions are
	  				        specified by spin indices
	  				        (input formats can vary)     */
For this spin system we (should) know the number of spins and their isotope
types.  The dipolar interactions are all to be placed into the dipolar
interaction vector, Dvec.  This must be read in using spin (NOT interaction)
indices.


For this spin system we ( should ) know the number of spins and their isotope types . The dipolar interactions are all to be placed into the dipolar interaction vector Dvec . This must be read in using spin ( NOT interaction ) indices .
void solid_sys::setDs(const ParameterSet& pset, int ccount)
  {
  vector<Isotope> Isos = IsoVec();		
  if(ccount && cflags.size())			
    Dvec = IntDipVec(Isos, SCoords, 0);		
  else						
    Dvec = IntDipVec(pset,Isos,-1,0);		
  }
void solid_sys :: setDs ( const ParameterSet & pset , int ccount ) { vector < Isotope > Isos = IsoVec ( ) ; if ( ccount && cflags . size ( ) ) Dvec = IntDipVec ( Isos , SCoords , 0 ) ; else Dvec = IntDipVec ( pset , Isos , -1 , 0 ) ; }
For this spin system we (should) know the number of spins and their isotope
types.
For this spin system we ( should ) know the number of spins and their isotope types .
// Get vector system isotopes // If we know spin coordinates // set interactions using them // If no spin coordinates // set interactions using Isos