aluque/grempy

Error using numpy rank

Opened this issue · 0 comments

Hi,
trying to use the library I ran into the following problem:

Traceback (most recent call last):
  File "./em.py", line 317, in <module>
     main()
  File "./em.py", line 283, in main
    sim.update_e()
  File "/home/kohlmann/Software/grempy/langevin.py", line 169, in update_e
    * avg(self.j_(R), axis=R))
  File "/home/kohlmann/Software/grempy/langevin.py", line 159, in j_
    slices = [s_[:]] * rank(self.j)
NameError: name 'rank' is not defined

I wanted to substitute it by linalg.matrix_rank(self.j), but it gave back:

Traceback (most recent call last):
  File "./em.py", line 317, in <module>
    main()
  File "./em.py", line 283, in main
    sim.update_e()
  File "/home/kohlmann/Software/grempy/langevin.py", line 169, in update_e
    * avg(self.j_(R), axis=R))
  File "/home/kohlmann/Software/grempy/langevin.py", line 159, in j_
    slices = [s_[:]] * linalg.matrix_rank(self.j)
TypeError: unsupported operand type(s) for *: 'slice' and 'int'