test2d_g.c

//  gcc -Wall -Wextra test2d_g.c -Xlinker -rpath=../lib
//  -L ../lib -ljigsaw -o test2d_g

//  Uses TRIPOD to build an (unconstrained) DT in E^2,
//  comparable to more conventional Delaunay tessellation
//  approaches.

#   include "../inc/lib_jigsaw.h"

#   include "print.h"
#   include "stdio.h"

    int test2d_g (int _verb)
    {
        int _retv = 0;

    /*-------------------------------- setup JIGSAW types */
        jigsaw_jig_t _jjig ;
        jigsaw_init_jig_t(&_jjig) ;

        jigsaw_msh_t _init ;
        jigsaw_init_msh_t(&_init) ;

        jigsaw_msh_t _tria ;
        jigsaw_init_msh_t(&_tria) ;

    /*-------------------------------- pts to triangulate */

        jigsaw_VERT2_t _init_vert2[9] = {
            { {0., 0.}, +0 } ,
            { {2., 0.}, +0 } ,
            { {1., 2.}, +0 } ,
            { {0., 1.}, +0 } ,
            { {.5, .0}, +0 } ,
            { {1., .7}, +0 } ,
            { {.6, 1.}, +0 } ,
            { {.0, .4}, +0 } ,
            { {.3, .3}, +0 }
            } ;

        _init._flags
            = JIGSAW_EUCLIDEAN_MESH;

        _init._vert2._data = &_init_vert2[0] ;
        _init._vert2._size = +9 ;

    /*-------------------------------- build TRIPOD r-DT. */

        _jjig._verbosity = _verb;

        _jjig._mesh_dims =   +2 ;

        _retv = tripod (
            & _jjig ,               // the config. opts
            & _init ,               // init. data
               NULL ,               // empty geom. data
            & _tria ) ;

    /*-------------------------------- print TRIPOD r-DT. */

        if (_verb > 0 )
        output_msh_data_2(&_tria);

        jigsaw_free_msh_t(&_tria);

        printf (
       "[2d_g] TRIPOD returned code : %d \n", _retv) ;


        return _retv ;
    }

#   ifndef __SKIP_MAIN__
    int main () { return test2d_g(1) ; }
#   endif//__SKIP_MAIN__