[Mohammad74]

Hello,
I've written a DEFINE_EXECUTE_ON_LOADING macro to find a cell in the domain closest to a specified point. I did it by looping over all the cells and finding the cell with minimum distance from that point. Later I will use this cell's identifiers (hubthread_id, hubcell_id in the code below) in a DEFINE_ADJUST macro from another source file to address it. Here's a shortened version of my code:

#include "udf.h"

Thread *hubthread_id ;
cell_t hubcell_id ;

DEFIINE_EXECUTE_ON_LOADING(macro_name,libname)
{
    Domain *dom ;
    Thread *c_thread ;
    cell_t cell ;
    double d, d_new, c[3], temp[3];
    double x_h[3] = {0.0,3.0,0.0} ;

    dom = Get_Domain(1) ;
    d = 1.2 ;

    thread_loop_c(c_thread,dom)
    {
        begin_c_loop(cell,c_thread)
        {
            C_CENTROID(c, cell, c_thread);
            d_new = sqrt(pow(x[0]-c[0],2)+pow(x[1]-c[1],2)+pow(x[2]-c[2],2));
            if (d_new < d)
            {
                hubthread_id = c_thread ;
                hubcell_id = cell ;
                d = d_new ;
            }
                  
        }
        end_c_loop(cell,c_thread)
    }

    C_CENTROID(temp, hubcell_id, hubthread_id);               /*I'm trying to see if its working properly, but gave me an error*/
    Message("\n %d %d %d", temp[0],temp[1],temp[2]) ;

}

I can compile it with no errors but after loading the library the last two lines are causing my code to break with these errors:

SharedScreenshot.jpg

Any reasons why this is happening? or any better ways to do this whole process?

[AlexanderZ]

no idea how you could compile it without errors, as there are few typos

I recommend you to compile code always
I was able to compile following code, test it:

#include "udf.h"

Thread *hubthread_id ;
cell_t hubcell_id ;

DEFINE_EXECUTE_ON_LOADING(macro_name,libname)
{
    Domain *dom ;
    Thread *c_thread ;
    cell_t cell ;
    real d, d_new, c[ND_ND], temp[ND_ND];
    real x[3] = {0.0,3.0,0.0} ;

    dom = Get_Domain(1) ;
    d = 1.2 ;

    thread_loop_c(c_thread,dom)
    {
        begin_c_loop(cell,c_thread)
        {
            C_CENTROID(c, cell, c_thread);
            d_new = sqrt(pow(x[0]-c[0],2)+pow(x[1]-c[1],2)+pow(x[2]-c[2],2));
            if (d_new < d)
            {
                hubthread_id = c_thread ;
                hubcell_id = cell ;
                d = d_new ;
            }
                  
        }
        end_c_loop(cell,c_thread)
    }

    C_CENTROID(temp, hubcell_id, hubthread_id);               /*I'm trying to see if its working properly, but gave me an error*/
    Message0("\n %f %f %f", temp[0],temp[1],temp[2]) ;

}

[Mohammad74]

Originally Posted by AlexanderZ

no idea how you could compile it without errors, as there are few typos

Hi,
I'm sorry, Yeah you are correct about typos. But i made these mistakes while typing the code in this forum.
Anyway, I compiled the code you provided and again it was compiled with no errors, but after loading the library, it gives a lot of errors and I couldn't find the reason behind it. Here are images of my tui after loading the library made from your code:

1.jpg2.jpg