h_constrained_gravity: h_constrained gravity
In CRC1266-A2/moin: Modelling Interactions

Constrained gravity model

1 2	h_constrained_gravity(E, Si = 1, Sj = 1, beta, C, c, gammas, g, margin = 1)

`E`	a matrix of random values (to be "optimized" by the Hamiltonian)
`Si`	a vectory of node attributes, e.g. carrying capacity or population
`Sj`	a vectory of node attributes, e.g. carrying capacity or population
`beta`	lagrange multiplier
`C`	distance matrix
`c`	distance constraining cost of the system
`gammas`	lagrange multiplier vector
`g`	constraining vector
`margin`	constraining dimension (either sum of rows or columns)