[OctDev] Re: Changes for the geometry package

[Brad B. <br...@sh...>]

At 10:57 AM 2/3/2004, Pedro Tarrafeta wrote:
<clip>

>When I set Qt in the options qh_new_qhull won't give a triangulated hull. As 
>far as I understand your code, qh_triangulate() is called before producing 
>the output, so I should just make a call to qh_triangulate(), is this right?.

Yes, you can call qh_triangulate().    It may be better to fix the underlying problem and create a function (see below) called qh_prepare_output().  

qh_produce_output() includes several housekeeping calls before calling the print functions.   In the next version of qhull, these calls will move to a new function called qh_prepare_output().   Qhull library users should call qh_prepare_output() if qh_new_qhull() succeeds.  The variables qh VORONOI, etc. are set from the option string.  

void qh_prepare_output(void) { 
  if (qh VORONOI) {
    qh_clearcenters (qh_ASvoronoi);
    qh_vertexneighbors();
  }
  if (qh TRIangulate) {
    qh_triangulate(); 
    if (qh VERIFYoutput && !qh CHECKfrequently) 
      qh_checkpolygon (qh facet_list);
  }
  qh_findgood_all (qh facet_list); 
  if (qh GETarea)
    qh_getarea(qh facet_list);
  if (qh KEEParea || qh KEEPmerge || qh KEEPminArea < REALmax/2)
    qh_markkeep (qh facet_list);
}

                                                --Brad