continuous.schmidt.decompose {RQEntangle}

Perform a continuous Schmidt decomposition

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)