Andrea Cassioli andreacassioli

## test.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                andreacassioli
                / test.ipynb
            
            
              Created
              August 19, 2016 08:18
            
              
                test
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## quadratic_fail.mps
* ENCODING=ISO-8859-1
NAME
ROWS
 N  obj
 L  c1
 G  c2
 L  c3
 G  c4
 L  c5
 G  c6

## gist:9114130
M = Model("minimal sphere enclosing a set of points dual")

print('Declaring the variables...')
y= M.variable('y',NDSet(k,n+1),Domain.inQCone(k,n+1))

c=[0. for i in range(n+1)]
c[0]=1.

print('Defining the constraints...')
M.constraint('equalities', Expr.mul(DenseMatrix(1,k,1.0), y), Domain.equalsTo(c) )

## gist:9110599
def sphere_enclosing_dual(n,k,p):
    print('Creating the Fusion optimization model...')

    start=time.clock()
    M = Model("minimal sphere enclosing a set of points dual")


    print('Declaring the variables...')
    y= M.variable('y',NDSet(k,n+1),Domain.inQCone(k,n+1))
    print('done!\n')

## plot.py
import sys

import mosek
from  mosek.fusion import *
import numpy as np
import time
import plot


def sphere_enclosing_primal(n,k,p):
	M = Model("minimal sphere enclosing a set of points dual")

	print('Declaring the variables...')
	y= M.variable('y',NDSet(k,n+1),Domain.inQCone(k,n+1))

	c=[0. for i in range(n+1)]
	c[0]=1.

	print('Defining the constraints...')
	M.constraint('equalities', Expr.mul(DenseMatrix(1,k,1.0), y), Domain.equalsTo(c) )
	def sphere_enclosing_dual(n,k,p):
	print('Creating the Fusion optimization model...')

	start=time.clock()
	M = Model("minimal sphere enclosing a set of points dual")


	print('Declaring the variables...')
	y= M.variable('y',NDSet(k,n+1),Domain.inQCone(k,n+1))
	print('done!\n')
	import sys

	import mosek
	from mosek.fusion import *
	import numpy as np
	import time
	import plot


	def sphere_enclosing_primal(n,k,p):