A search of SymPy *.py files for "reduced density" gives a hit in method `_reduced_density` which is called from the private method `Qubit._eval_trace` which is called by `Trace(qubit).doit()`. Does that get you any closer to a solution?
Chris /c On Sunday, July 31, 2022 at 1:36:47 AM UTC-5 hmaso...@ku.edu.tr wrote: > Hi, > I'm new to SymPy, and I'm in the middle of a calculation where I use the > groundstates of a Heisenberg model hamiltonian to write my density matrix. > I found the function Density() which rewrites the matrix as a superposition > of qubits, but then I am stuck. > Is there any function to calculate the reduced density matrix? > Or any function for partial trace? > Or any way to decompose, let's say, |000> to |0>|0>|0>? > > Any help is appreciated > Hamideh > -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/df23ae67-5b97-4211-8ab0-6d5f34c76113n%40googlegroups.com.