Simplex algorithms for 3D limit analysis of roman groin vaults