How to save a list to a file, with sage objects in it?

Question

Hi, I am trying to save intermediate step of my calculation into a file, so that next time I can read the file to resume the calculation. Say, I need to save

f = function('a',var('t'))
s = [f] # together with other things in the list

How to save s to a file?

I tried two ways:

(1) Standard python module pickle cannot save sage objects, thus fails

(2) I converted s to a python object:

s = Sequence([f])

Now I can use s.save('file_name'). However, there is a bug that I cannot load using load('file_name')

RuntimeError: unknown function 'a' in archive

It is a known bug: http://trac.sagemath.org/sage_trac/ticket/11919 .

Unfortunately I am not expert enough to fix that bug, and I really want to get my work done before waiting for the fix. Is there any other way to save the list s currently?

In the bug report it is mentioned that an older version of Pynac doesn't have this bug. If there is no other work around, is it possible to downgrade Pynac within sage or I have to downgrade the whole sage?

Thanks!

Accepted Answer

Okay, this really isn't a solution, but might serve as a temporary workaround-- you can't save formal functions, but you can convert them to maxima objects, and those can be saved.

For example:

sage: f = function('a',var('t'))
sage: parent(f)
Symbolic Ring
sage: # convert it to maxima (maxima(f) does the same thing)
sage: f._maxima_()
'a(t)
sage: # note the prefix ' now!
sage: parent(_)
Maxima_lib
sage: # convert forward and back
sage: (f._maxima_())._sage_()
a(t)
sage: # do we have equality?
sage: bool((f._maxima_())._sage_() == f)
True

So while we have

sage: loads(dumps(f))
---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
[...]
RuntimeError: unknown function 'a' in archive

we're lucky enough to have

sage: bool(loads(dumps(f._maxima_()))._sage_() == f)
True

It's straightforward to bundle the above technique for these guys. Not everything can be converted seamlessly between Sage and maxima, so you might want to only convert some things, but you get the idea.

Thanks a lot! Seems a best answer before the bug is fixed.
tririver (Dec 19 '11)

Answer 2

0	A partial answer: in Sage 4.7.0, the above attempt (2) works. Thus it is indeed a bug brought in from either 4.7.1 or 4.7.2. I also tried 4.8 alpha 5, the bug is not fixed --looking forward to a fix. link posted Dec 17 '11 tririver 282 ● 2 ● 4 ● 14 http://www.physics.mcgill...

How to save a list to a file, with sage objects in it?

2 Answers:

Your answer

Question tools

Tags:

Stats:

Related questions