[simonmcintyre]

Hi,

As the title says I'm trying to write a UDF to store the initial velocity components from a solution in a UDM variable (I need these intial velocities for a VG source term model I'm trying to implement). To do this I'm trying to use DEFINE_ON_DEMAND.

Currently I'm having trouble interpreting the UDF and receiving a segmentation error. Can anyone shed some light on why this might be happening?

I'm basing my code on the post here (Udm value on the walls).

I have already allocated 8 UDM locations and am running in serial, though would like to convert to parallel if I can get it working.


Code:

#include "udf.h"

int ID1 = 13;
int ID2 = 14;
int ID3 = 15;

DEFINE_ON_DEMAND(allocate_vel_initial)
{
Domain *d;
Thread *tc;
Thread *tf;
Thread *tA;
Thread *tB;
Thread *tC;
cell_t c;
face_t f;
d = Get_Domain(1);
tA = Lookup_Thread(d,ID1);
tB = Lookup_Thread(d,ID2);
tC = Lookup_Thread(d,ID3);
/* Loop through all the cell threads in the domain */
thread_loop_c(tc,d)
{
/* Loop through the cells in each cell thread */
begin_c_loop(c,tc)
{
C_UDMI(c,tc,0) = 1;
C_UDMI(c,tc,1) = 2;
C_UDMI(c,tc,2) = 3;
C_UDMI(c,tc,3) = 4;
}
end_c_loop(c,tc)
}

/* loops over all face threads in a domain*/
thread_loop_f(tf,d)
if (tf == tA) continue;
if (tf == tB) continue;
if (tf == tC) continue;
{
tc = THREAD_T0(tf); /*unique cell thread next to boundary face thread*/
/* Loop over faces in a face thread to get the information stored on faces. */
begin_f_loop(f,tf)
{
c = F_C0(f,tf); /*index of cell next to given face on boundary*/
F_UDMI(c,tf,4) = 5;
F_UDMI(c,tf,5) = 6;
F_UDMI(c,tf,6) = 7;
F_UDMI(c,tf,7) = 8;
}
end_f_loop(f,tc)
}
}


Error:

Error: received a fatal signal (Segmentation fault)
Error object: #f