Skip to main content

Matlab, nested functions, and variables

I rejoiced that MATLAB lets you have locally declared functions. Ahh, where are Python's namespaces? All hail Python. Fie promiscuous MATLAB! poor shadow, false pretender.. ahem..

Ok, so I go and nest a function and load some data from a file and Boom! MATLAB cries Uncle, with some wimpy excuse that it can't load variables into a static workspace. I ask you, is that grown up behavior?

So I have to do this round about thing: Suppose I have a mat file segmenting.mat with 6 variables called segment, dive_trials, dive_segment, trials, interaction, and capture. I have to load them into an explicit variable (thisisbogus) and then I have them as structure fields. Ewww.

thisisbogus = load('segmenting')

thisisbogus =

segment: [43x2 double]
dive_trials: [2 3 4 5 7 8 9 10 12 16 19 20 26 29 32 38 39 41 43]
dive_segment: [19x2 double]
trials: [1x38 double]
interaction: [43x1 double]
capture: [43x1 double]

I didn't want to type the bazillion lines of code to assign each field back to a variable, so I automated it

thisisbogus = load('segmenting')
bog = fieldnames(thisisbogus)
for n=1:length(bog),
disp([bog{n} ' = getfield(thisisbogus,bog{' num2str(n) '});'])
end

which prints

segment = getfield(thisisbogus,bog{1});
dive_trials = getfield(thisisbogus,bog{2});
dive_segment = getfield(thisisbogus,bog{3});
trials = getfield(thisisbogus,bog{4});
interaction = getfield(thisisbogus,bog{5});
capture = getfield(thisisbogus,bog{6});

Which can be pasted into the MATLAB function.

I get this nagging feeling that MathWorks is getting sloppy, I wonder why. This makes me upset, just like Microsoft Word's half-baked implementation of a cross-referencing system.

Then why don't I just go use Python and Latex then, you say? I do, I do, but JEB prefers Word files, and my co-authors use MATLAB and Word...

Comments

  1. FYI, you don't have to copy and paste that output to do your reassignment. You can execute the commands directly in your loop if you just change the "disp" to an "eval".

    ReplyDelete
  2. Unfortunately, no. The evals will require MATLAB to generate variables on the fly, which it can not do in nested functions (aka static workspaces)

    ReplyDelete
  3. Good point. How about an "evalin('caller',..." then? I got it to work as follows:

    function output=nestTest()
    nest()
    function []=nest()
    for i = 0:25
    evalin('caller',['output.','a'+i,'=',num2str(i),';'])
    end
    end
    end


    Yielding:

    >> nestTest()
    ans =
    a: 0
    b: 1
    c: 2 ...

    ReplyDelete

Post a Comment

Popular posts from this blog

A note on Python's __exit__() and errors

Python's context managers are a very neat way of handling code that needs a teardown once you are done. Python objects have do have a destructor method ( __del__ ) called right before the last instance of the object is about to be destroyed. You can do a teardown there. However there is a lot of fine print to the __del__ method. A cleaner way of doing tear-downs is through Python's context manager , manifested as the with keyword. class CrushMe: def __init__(self): self.f = open('test.txt', 'w') def foo(self, a, b): self.f.write(str(a - b)) def __enter__(self): return self def __exit__(self, exc_type, exc_val, exc_tb): self.f.close() return True with CrushMe() as c: c.foo(2, 3) One thing that is important, and that got me just now, is error handling. I made the mistake of ignoring all those 'junk' arguments ( exc_type, exc_val, exc_tb ). I just skimmed the docs and what popped out is that you need to return True or...

Store numpy arrays in sqlite

Use numpy.getbuffer (or sqlite3.Binary ) in combination with numpy.frombuffer to lug numpy data in and out of the sqlite3 database: import sqlite3, numpy r1d = numpy.random.randn(10) con = sqlite3.connect(':memory:') con.execute("CREATE TABLE eye(id INTEGER PRIMARY KEY, desc TEXT, data BLOB)") con.execute("INSERT INTO eye(desc,data) VALUES(?,?)", ("1d", sqlite3.Binary(r1d))) con.execute("INSERT INTO eye(desc,data) VALUES(?,?)", ("1d", numpy.getbuffer(r1d))) res = con.execute("SELECT * FROM eye").fetchall() con.close() #res -> #[(1, u'1d', <read-write buffer ptr 0x10371b220, size 80 at 0x10371b1e0>), # (2, u'1d', <read-write buffer ptr 0x10371b190, size 80 at 0x10371b150>)] print r1d - numpy.frombuffer(res[0][2]) #->[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] print r1d - numpy.frombuffer(res[1][2]) #->[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] Note that for work where data ty...