tests/test1.R

In Rdsdp: R Interface to DSDP Semidefinite Programming Library

# Test a simple example (Sanity check)
library(Rdsdp)

	# Sedumi format example
	K=c()
	K$s=c(2,3)
	K$l=2

	C=matrix(c(0,0,2,1,1,2,c(3,0,1,
                       0,2,0,
                       1,0,3)),1,15,byrow=TRUE)
	A=matrix(c(0,1,0,0,0,0,c(3,0,1,
                         0,4,0,
                         1,0,5),
          	1,0,3,1,1,3,rep(0,9)), 2,15,byrow=TRUE)
	b <- c(1,2)
        C=-C

    ret = dsdp(A,b,C,K)
    stopifnot(all.equal(ret$y,c(-1.0,-0.75),tolerance=1e-05))