a quick question about scheduling:

I have a scalar variable, which contains the elements of a rotation matrix that rotates the z-axis into a given vector...

I want to calculate some quantities (that are grid functions) in coordinates that have been rotated by this matrix...

I run into the following problem:

when scheduling the calculation of the quantities in local mode, the function doesn't seem to have access to the rotation matrix elements... 

when scheduling the calculation in global, loop-local mode, I can't output the data at the grid points (the output gives 1.96808407167164e+243 at each point for every refinement level apart from the innermost, where the values are output correctly)

I guess I could make the rotation matrix elements a grid function, but that seems to be a bit of an overkill with regard to the memory needed for the 9 CCTK_REALs i need to store the matrix as a SCALAR..

Any ideas on how to do this?

best wishes,

Vassili