src/Bloch/BlochM.cc
/* These are just direct functions which return 3Nx1 column vectors (as MagVec)
   Them may contain multiple sub-vectors. They are formally generated by the
   formula
                                         R
                               |M  > =   - |M  >
                                 ss      G   eq

   but of course this will have troubles whenever the matrix R is zero or the
   matrix G cannot be easily inverted. When R is zero (no-relaxation) then
   there simply is no steady-state defined. Alternativly, if R=G (i.e. there
   are no rf-fields appled) the the steady state is one and the same as the
   equilibrium state.                                                        */
These are just direct functions which return 3Nx1 column vectors (as MagVec)
Them may contain multiple sub-vectors. They are formally generated by the
formula
R
|M  > =   - |M  >
ss      G   eq

but of course this will have troubles whenever the matrix R is zero or the
matrix G cannot be easily inverted. When R is zero (no-relaxation) then
there simply is no steady-state defined. Alternativly, if R=G  the the steady state is one and the same as the
equilibrium state.
These are just direct functions which return 3Nx1 column vectors ( as MagVec ) Them may contain multiple sub - vectors . They are formally generated by the formula R |M > = - |M > ss G eq but of course this will have troubles whenever the matrix R is zero or the matrix G cannot be easily inverted . When R is zero ( no - relaxation ) then there simply is no steady - state defined . Alternativly if R = G the the steady state is one and the same as the equilibrium state .
MagVec BlochMss(const matrix& G, const matrix& R, const MagVec& Meq)
  {
  if(R.is_zero()) return Meq;
  return (inv(G)*R)*Meq;
  }
MagVec BlochMss ( const matrix & G , const matrix & R , const MagVec & Meq ) { if ( R . is_zero ( ) ) return Meq ; return ( inv ( G ) * R ) * Meq ; }
These are just direct functions which return 3Nx1 column vectors (as MagVec)
Them may contain multiple sub-vectors.
These are just direct functions which return 3Nx1 column vectors ( as MagVec ) Them may contain multiple sub - vectors .