continuous.schmidt.decompose: Perform a continuous Schmidt decomposition

In RQEntangle: Quantum Entanglement of Bipartite System

Description

Perform a continuous Schmidt decomposition

Usage

continuous.schmidt.decompose(bifunc, x1lo, x1hi, x2lo, x2hi, nbx1 = 100,
  nbx2 = 100, keep = min(10, nbx1, nbx2))

Arguments

bifunc	bipartitite continuous wavefunction
x1lo	lower limit of x1
x1hi	upper limit of x1
x2lo	lower limit of x2
x2hi	upper limit of x2
nbx1	number of discretized x1 (default: 100)
nbx2	number of discretized x2 (default: 100)
keep	number of Schmidt modes to keep (default: minimum of 10, nbx1, and nbx2)

Value

Schmidt modes, including the eigenvalues, and the lambda interpolated function of the Schmidt modes

Examples

coupled.harm.fcn<- function(x1,x2) exp(-((0.5*(x1+x2))**2))*exp(-(x1-x2)**2)*sqrt(2./pi)
continuous.schmidt.decompose(coupled.harm.fcn, -10, 10, -10, 10)

RQEntangle documentation built on May 2, 2019, 3:01 a.m.