hDiaphony: function to do ...
In quaddisc: Quadratic discrepancies

Description Usage Arguments Details Value Note Author(s) References See Also Examples

~~ A concise (1-5 lines) description of what the function does. ~~

1	hDiaphony(vX, vY, lExtraArgsKernel = NULL)

`vX`	~~Describe `vX` here~~
`vY`	~~Describe `vY` here~~
`lExtraArgsKernel`	~~Describe `lExtraArgsKernel` here~~

~~ If necessary, more details than the description above ~~

~Describe the value returned If it is a LIST, use

`comp1`	Description of 'comp1'
`comp2`	Description of 'comp2'

...

~~further notes~~

~Make other sections like Warning with Warning .... ~

~~who you are~~

~put references to the literature/web site here ~

~~objects to See Also as help, ~~~

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (vX, vY, lExtraArgsKernel = NULL) 
{
    f <- function(vX) -1 + prod(1 + 2 * pi^2 * (vX^2 - vX + 1/6))
    iD <- length(vX)
    vZ <- vX - vY
    for (i in 1:iD) vZ[i] <- mod(vZ[i], 1)
    dResult <- (1/((1 + pi^2/3)^iD - 1)) * f(vZ)
    return(dResult)
  }