CGAL users discussion list

["Christopher O'Donnell" <>]

I'm trying to input a Polyhedron_3 using the built-in .off input stream. All of the points seem to be input correctly, but all of the Plane_3 objects are read as (0.5, 0.5, 0.5, 0.5). Also, I get lots of strange errors when I try to manipulate them (inf > sup assertion failure and seg faults).

Here is the code I'm using, and it fails on the corner.off from examples/Polyhedron:

#include <CGAL/Simple_cartesian.h>
#include <CGAL/Polyhedron_3.h>
#include <CGAL/Exact_predicates_exact_constructions_kernel.h>

#include <CGAL/IO/Polyhedron_iostream.h>
#include <iostream>
#include <fstream>

typedef CGAL::Exact_predicates_exact_constructions_kernel Kernel;
typedef CGAL::Polyhedron_3<Kernel, CGAL::Polyhedron_items_3> Polyhedron;

typedef Polyhedron::Point_iterator Point_iterator;
typedef Polyhedron::Plane_iterator Plane_iterator;

int main(int argc, char** argv)
{
    Polyhedron P;
    std::ifstream input(argv[1]);
    input >> P;

    for (Point_iterator i = P.points_begin(); i != P.points_end(); ++i)
        std::cout << *i << std::endl;
    for (Plane_iterator i = P.planes_begin(); i != P.planes_end(); ++i)
        std::cout << *i << std::endl;

    fflush(stdout);
    return 0;
}

This is the output I get from corner.off:

0 0 0
1 0 0
2 0 0
0 0 1
1 0 1
2 0 1
0 0 2
1 0 2
0 1 0
1 1 0
2 1 0
0 1 1
1 1 1

2 1 1
0 1 2
1 1 2
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5

0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5
0.5 0.5 0.5 0.5

Any help would be appreciated.

Chris