I've been trying to use a vector

interface.ccl:
real cm[num_wave_components] TYPE=GF Timelevels=1 tags='Prolongation="None"'
real sm[num_wave_components] TYPE=GF Timelevels=1 tags='Prolongation="None"'

storage.ccl:
STORAGE: cm
STORAGE: sm

When the code runs, I find that cm is not allocated (assert(cm != 0) fails).

I found this syntax elsewhere, e.g. CarpetCG

carpet/CarpetDev/CarpetCG/interface.ccl:
real dirvars[cg_maxsolvevars] TYPE=GF Timelevels=1 
tags='Prolongation="None"'

carpet/CarpetDev/CarpetCG/schedule.ccl:
STORAGE: dirvars

I've tried a number of variations of the syntax. I've discovered that 
removing the square brackets
allows the storage allocation to happen (see below), but I really do 
want the vector.

real cm TYPE=GF Timelevels=1 tags='Prolongation="None"'
real sm TYPE=GF Timelevels=1 tags='Prolongation="None"'

Note, this is an old Carpet Nov 2012 (newer carpets sync differently and 
my ghost zones don't update).

Any idea what am I doing wrong?

Cheers,
Steve
_______________________________________________
Users mailing list
[email protected]
http://cactuscode.org/mailman/listinfo/users

Reply via email to