Merge 675bbd9 into cad6f8a

src/propmachinery.jl

@@ -173,6 +173,28 @@ QuEvolutionOp{QE<:QuEquation}(eq::QE, dt::Float64) = QuEvolutionOp(operator(eq),
 
 QuEvolutionOp{QE<:QuEquation}(eq::QE, tf::Float64, ti::Float64) = QuEvolutionOp(eq, tf-ti)
 
+function Base.show(io::IO, qprop::QuPropagator)
+    field_params = fieldnames(qprop.eq)
+    println(io, "Summarizing the system :")
+    if :lindblad in field_params && :hamiltonian in field_params
+        println(io, "Equation type : $(coeffs(qprop.eq.hamiltonian))")
+        println(io, "Lindblad operator of the system : $(coeffs(qprop.eq.lindblad))")
+        println(io, "Hamiltonian of the system : $(coeffs(qprop.eq.hamiltonian))")
+        println(io, "Collapse operators of the system : $(qprop.eq.collapse_ops)")
+        println(io, "Density matrix : $(coeffs(qprop.init_state))")
+    elseif :hamiltonian in field_params
+        println(io, "Equation type : $(typeof(qprop.eq))")
+        println(io, "Hamiltonian of the system : $(coeffs(qprop.eq.hamiltonian))")
+        println(io, "Initial state : $(coeffs(qprop.init_state))")
+    elseif :liouvillian in field_params
+        println(io, "Equation type : $(typeof(qprop.eq))")
+        println(io, "Liouvillian of the system : $(coeffs(qprop.eq.liouvillian))")
+        println(io, "Density matrix : $(coeffs(qprop.init_state))")
+    end
+    println(io, "Time steps used : $(qprop.tlist)")
+    println(io, "Solver used : $(qprop.method)")
+end
+
 export  QuStateEvolution,
       QuPropagator,
       QuEvolutionOp,