CGAL users discussion list

[Kwok Jasper <>]

Excuse me, after some further investigation, I found that I get the exception when I insert a Point object into  a list
(i.e.) the 2 lines of code below 

            co_input.push_back(Point(tempX2 - tempX1 , tempY2 - tempY1, tempZ2 - tempZ1));
            co_input.push_back(p);

inside the following code.


    Point_3 * co_mean = NULL;
    list<Point> co_input;
    list<Point>::iterator co_itr;


vector<Vertex_handle>::iterator nv_itr;

    double temp_fifth_nearest_dist;
    double tempX1;
    double tempY1;
    double tempZ1;
    double tempX2;
    double tempY2;
    double tempZ2;
    double dist;


    for(curr = start; curr != end; ++curr)
    {
        if (co_mean != NULL)
        {
           delete co_mean;
           co_mean = NULL;
        }

        if (co_input.size() != 0)
        {
           co_input.clear();
        }

        x_total = 0;
        y_total = 0;
        z_total = 0;

        x_variance = 0;
        y_variance = 0;
        z_variance = 0;
 
        x_y_covariance = 0;
        y_z_covariance = 0;
        x_z_covariance = 0;
        temp_fifth_nearest_dist = curr -> fifth_nearest_dist;

        cout << temp_fifth_nearest_dist << endl;



        nv_itr = (curr -> neighbor_vertex_handle).begin();


        while(1)
        {

            if (nv_itr == (curr -> neighbor_vertex_handle).end())
            {
               break;
            }
       
            tempX1 = (curr -> point()).x();
            tempY1 = (curr -> point()).y();
            tempZ1 = (curr -> point()).z();

            cout << "b0" << endl;

            tempX2 = ((*nv_itr) -> get_point()).x();
            tempY2 = ((*nv_itr) -> get_point()).y();
            tempZ2 = ((*nv_itr) -> get_point()).z();
           

            cout << "b1" << endl;

            dist = pow((tempX1 - tempX2),2);
          
            dist += pow((tempY1 - tempY2),2);
            dist += pow((tempZ1 - tempZ2),2);
            dist = sqrt(dist);
           
            co_input.push_back(Point(tempX2 - tempX1 , tempY2 - tempY1, tempZ2 - tempZ1));
            co_input.push_back(p);

           
            ++nv_itr; 
   }


May I ask if is possible the problem actually relates to the use of

    list<Point> co_input;
    list<Point>::iterator co_itr;

i.e. the std::list object with the Point object in CGAL?
And even I change the std::list into std::vector, the same exception still exist..


Thank you


           
         

---

From:
To:
Date: Sun, 8 Feb 2009 16:06:38 +0800
Subject: RE: [cgal-discuss] [HELP] about some advice related to the usage of

Thanks.

Actually, I used the following to declare the triangulation

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;


typedef CGAL::Triangulation_data_structure_3<Vertex_base_for_ssr<K> > Tds;
typedef CGAL::Delaunay_triangulation_3<K, Tds> Delaunay_triangulation;

typedef Delaunay_triangulation::Vertex_handle Vertex_handle;
typedef Delaunay_triangulation::Point Point;
typedef Delaunay_triangulation::Finite_vertices_iterator
Finite_vertices_iterator;
typedef Delaunay_triangulation::Finite_facets_iterator
Finite_facets_iterator;

typedef Delaunay_triangulation::Triangle Triangle;
typedef Delaunay_triangulation::Segment Segment;


I am using the exact kernel.

Is it possibly, I need to reinstall CGAL to solve the problem?
Or is there still any other possible reason for the probelm to occur?

Thank you very much
From: laurent.rineau__
To: 
Date: Sat, 7 Feb 2009 18:33:53 +0100
Subject: Re: [cgal-discuss] [HELP] about some advice related to the usage of

On Saturday 07 February 2009 17:58:49 Kwok Jasper wrote:
> Excuse me, for the code in the question that I have asked before,
>
>
>
> may I ask if I have understood the CGAL manual properly an used the
> structure and function there correctly.
>
> Although the my whole code can be compiled, it get exception


Perhaps you do not use an filtered or exact kernel.

 

-- 
Laurent Rineau, PhD
Engineer at GeometryFactory
http://www.geometryfactory.com/

---

Windows Live Gallery 讓才華綻放，讓全世界看見你的創造。

---

內在氣質才是致勝之道－您是哪一 最原始的溝通是用心的交流。