Ask Your Question
2

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

asked 2011-12-15 19:59:22 +0100

updated 2015-01-14 14:19:08 +0100

FrédéricC gravatar image

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/ti... .

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!

edit retag flag offensive close merge delete

Comments

I don't think it's the fact that `s` is a list, it's that it contains the function `f`. If you had "s = [sin(x), 12, 'hello']", you could do "save(s, 'file_name')" and then "load('file_name.sobj')" works fine.

John Palmieri gravatar imageJohn Palmieri ( 2011-12-16 00:29:37 +0100 )edit

Hi, John: thanks a lot for your reply! I agree with you. But my actual goal is really want to save a list with a lot of functions like f. In my project, I start from an abstract form of a differential equation, and simplify it to a simpler form (here I want to save my work because the former step is slow, and I want to use it when I start Sage next time), and solve the differential equation with different additional conditions. It's a differential equation thus I have to define a general function instead of put a "sin(x)" in the list.

tririver gravatar imagetririver ( 2011-12-16 06:52:44 +0100 )edit

2 Answers

Sort by » oldest newest most voted
0

answered 2011-12-17 22:02:59 +0100

DSM gravatar image

updated 2011-12-17 22:39:05 +0100

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.

edit flag offensive delete link more

Comments

Thanks a lot! Seems a best answer before the bug is fixed.

tririver gravatar imagetririver ( 2011-12-19 09:00:23 +0100 )edit
0

answered 2011-12-17 19:53:37 +0100

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.

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 2011-12-15 19:59:22 +0100

Seen: 1,800 times

Last updated: Dec 17 '11