gill1109 wrote:...regarding formula (13), in good mathematical typography the function "sign" should be set in Roman text, not in italic math symbols, c.f. "log", "sin" and so on. Secondly, in view of the definition of the domain, the second line of (13) is superfluous and indeed misleading. The outcomes of measurements are +/-1. There never is an outcome "0" because there wouldn't have been a state if the outcome had been 0.
Joy Christian wrote:Lucien Hardy ... asked: What happens when? Well, nothing happens. There is a more detailed discussion about this ambiguity in notation in this chapter: http://libertesphilosophica.info/blog/w ... hapter.pdf [cf. eqs. (A.9.48) and (A.9.49) on page 254 and eqs. (A.9.53) and (A.9.54) on page 256]. The difference between these two sets of prescriptions is the difference between ontology and computability.
gill1109 wrote:My answer would have been: "it doesn't happen".
Joy Christian wrote:Fine with me. The bottom line is that there exist no states for whichfor any
, so the question of corresponding outcomes does not arise.
gill1109 wrote:Perhaps you can comment then on the two simulations I have just posted: one is a copy of Michel's, the other "converts" Michel's scheme fromto
. Both exhibit strong violations of CHSH but in fact the sample sizes in my simulation are large enough that we can be certain that neither model reproduces the singlet correlation (cosine difference of angles). They are both systematically off target, the "full" simulation
is further off than the restricted one
.
Conclusion: either your model is "wrong", or Michel's code does not simulate your model, but something else instead.
gill1109 wrote:They do give results systematically slightly different from the singlet correlation, provided one implements his algorithm with higher numerical precision and with larger sample sizes.
alpha <- 0 *2* pi / 360 ## CHSH: try 0 and 90 degrees
beta <- 45 *2* pi / 360 ## CHSH: try 45 and 135 degrees
a <- c(cos(alpha), sin(alpha)) ## S2 version: c(cos(alpha), sin(alpha), 0)
b <- c(cos(beta), sin(beta)) ## S2 version: c(cos(beta), sin(beta), 0)
sumprod <- 0
N <- 0
for (i in 1:10^6) {
e0 <- rnorm(2) ## S2 version: rnorm(3)
e0 <- e0/sqrt(sum(e0^2))
theta0 <- runif(1, 0, pi/2)
ca <- sum(e0*a)
cb <- sum(e0*b)
s <- (sin(theta0)^2)/2
if (abs(ca) > s & abs(cb) > s ) {
sumprod <- sumprod + sign(ca)*sign(cb)
N <- N+1
}
}
sum(a*b) ## Theoretical correlation
N ## Number pairs of particles
sumprod/N ## Observed correlation
> sum(a*b) ## Theoretical correlation
[1] 0.7071068
> N ## Number pairs of particles
[1] 690569
> sumprod/N ## Observed correlation
[1] 0.700198
> sum(a*b) ## Theoretical correlation
[1] 0.7071068
> N ## Number pairs of particles
[1] 601609
> sumprod/N ## Observed correlation
[1] 0.7071985
alpha <- 0 * 2* pi / 360
beta <- 45 * 2 * pi / 360
a <- c(cos(alpha), sin(alpha))
b <- c(cos(beta), sin(beta))
M <- 10^8
t <- runif(M, 0, 2*pi)
x <- cos(t)
y <- sin(t)
e <- rbind(x,y)
ca <- colSums(e*a)
cb <- colSums(e*b)
theta <- runif(M, 0, pi/2)
s <- (sin(theta)^2) / 2
good <- abs(ca) > s & abs(cb) > s
N <- sum(good)
corr <- sum(sign(ca[good])*sign(cb[good]))/N
> corr
[1] 0.6995283
> sqrt(1/N)
[1] 0.0001203852
> sum(a*b)
[1] 0.7071068
alpha <- 0 * 2 * pi / 360
beta <- 45 * 2 * pi / 360
a <- c(cos(alpha), sin(alpha), 0)
b <- c(cos(beta), sin(beta), 0)
M <- 10^8
z <- runif(M, -1, 1)
t <- runif(M, 0, 2*pi)
r <- sqrt(1 - z^2)
x <- r * cos(t)
y <- r * sin(t)
e <- rbind(x,y,z)
ca <- colSums(e*a)
cb <- colSums(e*b)
theta <- runif(M, 0, pi/2)
s <- (sin(theta)^2) / 2
good <- abs(ca) > s & abs(cb) > s
N <- sum(good)
corr <- sum(sign(ca[good])*sign(cb[good]))/N
> corr
[1] 0.7052408
> sum(a*b)
[1] 0.7071068
> sqrt(1/N)
[1] 0.0001289189
gill1109 wrote:OK, so it is now established that Michel Fodje (minkwe)'s simulation based on picking e_0 from S^1 does not exactly reproduce the cosine correlation we are after. Neither does an S^2 simulation....
Return to Sci.Physics.Foundations
Users browsing this forum: No registered users and 4 guests