Create a constraints matrix for triangulation

Hi everybody,

I would like to create a 2D triangulation with a lot of free borders. I think that the easiest way to do that is to use the DelaunayTri function.

I believe that the way to create free border is to use a constraint function and the inOutstatus of the triangulation :

    dt0=DelaunayTri(x0,C);
    inside0 = dt0.inOutStatus();
    tr0 = TriRep(dt0(inside0, :), dt0.X);

But I don't know how to create the Constraint matrix !

I have a n*n matrix of the edges I would like to use. How can I convert this matrix in a constraint matrix?

I hope you understand my problem and hopefully you will help me ! Thanks

[Olivier]

I tried to create the constraint matrix by looking the closer point.

My algorithm is very long to compute and is full of bugs. I am not sure to be able to create the constraint matrix by this way.

Is there any existing hints to know when you create the constraint matrix from a complex geometry ?