示例#1
0
文件: propagator.py 项目: foosa/qudy
 def __init__(self, *args, **keyword_args):
     
     # If there is one argument, then the user only gave a control
     # function.  Assume a qubit system, i.e. SU(2^n).
     if len(args) == 1:
         self.ideal_control = args[0]
         dim = self.ideal_control.number_controls
             
         # Convert Lie dimensonality into number of qubits.
         n = log( dim + 1 ) / log( 4 )
         self.number_qubits = n
         self.lie_algebra = "su(%i)" %(int(2**n))
         if n == 1:
             self.hamiltonians = [Hx,Hy,Hz]
         else:
             self.hamiltonians = routines.product_operator( n )
         
     # If there are two arguments, then the user gave both a set of
     # controls and also a set of Hamiltonians.
     elif len(args) == 2:
         self.ideal_control = args[0]
         self.hamiltonians = args[1]
         
     # If there are three arguments, then the user gave controls,
     # Hamiltonians and Lindblad operators.
     elif len(args) == 3:
         self.ideal_control = args[0]
         self.hamiltonians = args[1]
         self.lindblad = args[2]
         
         # Since Lindblad operators were specified, we should use
         # Master equation methods.  Calculate superoperator
         # representations of relevant operators.
         
     if not len( self.ideal_control ) == len( self.hamiltonians ):
         raise ValueError("Control dimension mismatch.  Controls and " + \
                          "Hamiltonians must be of the same length.")
         
     # Carry over several constants from controls
     self.dimension  =  self.ideal_control.dimension
     self.times      =  self.ideal_control.times
     self.timemin    =  self.ideal_control.timemin
     self.timemax    =  self.ideal_control.timemax
     
     # Create an control.  For propagator instances, these are
     # identical to ideal_control, however for imperfect instances
     # these differ.
     self.control = self.ideal_control.copy()
     
     # Parse through keyword arguments.  Sets default solution
     # method.  Other keywords that are not understood will be
     # quietly ignored.
     if keyword_args.has_key( 'solution' ):
         
         method = keyword_args['solution']
         valid_inputs = ['trotter', 'dyson', 'magnus', 'lindblad']
         
         if method in valid_inputs:
             self.solution_method = method
             
         else:
             # Keyword was not valid.  Check to see if already
             # defined.
             
             # Return to default and warn user
             self.solution_method = 'trotter'
             warn('Solution method not understood, defaulting to \'trotter\'.')
     else:
         # set default solution method
         self.solution_method = 'trotter'
         
     # Set order of pertubation theory
     default_order = 4
     if keyword_args.has_key( 'order' ):
         
         order = keyword_args['order']
         
         try:
             if (order % 1.0) == 0:
                 self.order = order
         
             else:
                 # Order was not an integer, revert to defaults
                 warn('Pertubation theory order must be integer valued.')
                 self.order = default_order
                 
         except TypeError:
             # Order was not numeric, revert to defaults
             warn('Pertubation theory order must be integer valued.')
             self.order = default_order
             
     else:
         # Set default order
         self.order = default_order
示例#2
0
文件: propagator.py 项目: foosa/qudy
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
# 
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <http://www.gnu.org/licenses/>.

from quantop import *
import control
import integration
import routines
import imperfect


__all__ = ['propagator','rotation','R']

[Hx,Hy,Hz] = routines.product_operator(1)

class propagator:
    """
    A class for quantum dynamical propagators.  The evolution is
    specified by a bilinear control system, and optionally, a set of
    Linblad decay channels.
    
    Inputs may be an instance of the control class and also lists of N
    x N matrices, where N is the dimensionality of the Hilbert space.
    If no Lindblad channels are specified, then the evolution is
    assumed to be unitary.  If no Hamiltonians are specified and the
    dimensionality of the control system is n(n-1), then we assume the
    dynamical Lie algebra is SU(2^n) and use a product-operator basis
    for the