sympy.geometry Point class is working slow

Question

I have a code which reads unstructured mesh. I wrote wrappers around geometric entities of sympy.geometry such as:

class Point:
    def __init__(self, x, y, parent_mesh):
        self.shape = sympy.geometry.Point(x,y)
        self.parent_mesh = parent_mesh
        self.parent_cell = list()

Everything works fine but initialization of sympy.geometry.Point takes a lot of time for each Point. Actually, the code did not finish execution for thousands of points. Similar code written in C++ finished in a few seconds. Without it the code is fast enough (I removed it and timed). I read that a possible reason could be that sympy.geometry converts floating point numbers to rationals for precision. Is there a way (flag) to speed up sympy.geometry as I do not need exact precision?

score 4 · Accepted Answer · answered Aug 08 '16 at 13:20

4

Take a look at the Point class documentation, specifically, in one of the first examples:

Floats are automatically converted to Rational unless the evaluate flag is False.

So, you could pass a flag named evaluate during initialization of your Point classes:

self.shape = sympy.geometry.Point(x,y, evaluate=False)

which apparently signals what you're after.

answered Aug 08 '16 at 13:20

Dimitris Fasarakis Hilliard

150,925
31
268
253

1

Yes this solved the problem. I should have read the documentation better. Thanks. – Shibli Aug 08 '16 at 15:36
Although now sympy.geometry.Triangle is working very slow and there is nothing about `evaluate` in this case. – Shibli Aug 08 '16 at 16:08
`Triangle` can be initialized with 3 `Point`s so the best you could do is make those Points not convert to Rationals. – Dimitris Fasarakis Hilliard Aug 08 '16 at 16:19
1

I did but no improvement. Also as the number of digits of points increase execution time increases dramatically. – Shibli Aug 08 '16 at 16:34

sympy.geometry Point class is working slow

1 Answers1

Linked