Solve instances with coordinates

[jkhamphousone]

Currently, the solver only solves a subset of TSPLIB adapted to RSP instances. We would like to make it solve instances with two matrices R and S as an input (ring costs) and (star costs):

R = [
     0  16   7  14
    16   0   3   5
     7   3   0  16
    14   5  16   0 
]
S = [
     0   6   7  14
    16   0   4   5
     7   3   0   1
    14   5   1   0 
]
RSP.rspoptimize(pars, R, S)```

[jkhamphousone]

closed by #16