topic Re: How to use call dynamic_array right? in New SAS User

How to use call dynamic_array right?

whymath — Tue, 10 Dec 2019 14:59:28 GMT

I have try to write a user-defined function MinDis() and apply it in data step, this function is used to compute the minimum distance from one point to each element of an (numeric)array. Code fllowing:

proc fcmp outlib = work.funcs.Math;
    function MinDis(Exp,Arr[*]);
        array dis[2] /nosymbols;
        call dynamic_array(dis,dim(Arr));
        do i = 1 to dim(Arr);
            dis[i] = abs((Exp - Arr[i]));
        end;
        return(min(of dis[*]));
    endsub;
quit;

option cmplib=work.funcs ; 

data MinDis;
    input LamdazLower LamdazUpper @;
    cards;
    2.50 10.0
    2.51 10.8
    2.49 9.97
    2.75 9.50
    ;
run;

data _null_;
    set ;

    array _PTN_ [14] _temporary_ (0.5,1,1.5,2,2.5,3,4,5,6,7,8,9,10,12);
    StdLamZLow = MinDis(LamdazLower,_PTN_);
    StdLamZUpp = MinDis(LamdazUpper,_PTN_);
    put _all_;
run;

It was compile rightly but gave wrong results. StdLamZLow just get the minimum distance from LamdazLower to the first two element of array _PTN_.

When I rewrite the dim of dis as 999 or something very big and get rid of call dynamic_array statement I would get it right. But I surely want to know why min(of dis[*]) just take dis as a 2-dim array.

By the way, how can I use implied DO-loops do over ... to instead of explicit DO loops? I have tried several times but haven`t success yet.

Thanks for any hints.

Re: How to use call dynamic_array right?

PeterClemmensen — Tue, 10 Dec 2019 13:05:00 GMT

A few things:

- What is your desired result given _PIN_? Probably makes it easier to help you.

- I assume that you misspelled nosymbols as symbols in your Array Statement?

- You should always use the Varargs Option when you want the function to accept a variable number of arguments (like an array)

- Do Over syntax is for implicit arrays only. PROC FCMP does not support implicit arrays.

Re: How to use call dynamic_array right?

whymath — Tue, 10 Dec 2019 14:58:18 GMT

1 The desired result of StdLamZLow should be the minimum distance from LamdazLower to each element of array _PTN_, that is minimum of (abs(LamdazLower-_PTN_[i])), am I more clear now?

2. Yes, I misspelled the option nosymbols. Now it is corrected.

3. I don't realize it is neccessary cause I just get an rig

4. Thank you for the fourth tip.

Re: How to use call dynamic_array right?

KachiM — Wed, 11 Dec 2019 08:50:45 GMT

@whymath

The issue here is that OF operator for a Dynamically created Array in FCMP which is NOT supported. There are two ways to solve your problem - one using Dynamic array and the other without it. In the former, we compare absolute difference of the previous and the current to choose the minimum. In the second we use MIN function to choose the minimum.

Let us see the first solution.


proc fcmp outlib = work.funcs.Math;
    function MinDis(Exp,Arr[*]);
        array dis[1] / nosym;
        call dynamic_array(dis,dim(Arr));
        mindis = constant('BIG');
        do i = 1 to dim(Arr);
            dis[i] = abs((Exp - Arr[i]));
            if mindis > dis[i] then mindis = dis[i]; * get minimum distance;
        end;
        return(mindis);
    endsub;
quit;

NOSYM is a short form for NOSYMBOLS. MINDIS is assigned to an impossible big positive value as Minimum Distance. Each absolute difference is compared to previous to get the current minimum. Finally, MINDIS has the final minimum value which is returned to the Data Step. The array DIS[ ] is a waste of memory and hence the second solution.

The second solution avoids the use of the array, DIS[ ] but uses the _temporary_ array passed by the Data Step. Here we compare the previous MINDIS with the i-th absolute difference to choose the minimum value.

/** No Dynamic_Array **/

proc fcmp outlib = work.funcs.Math;
    function MinDis(Exp,Arr[*]);
        mindis = Constant('BIG');
        do i = 1 to dim(Arr);
            mindis = min(mindis, abs(Exp - Arr[i]));
        end;
        return(mindis);
    endsub;
quit;

Hope this answers your question.

Re: How to use call dynamic_array right?

whymath — Thu, 05 Mar 2020 14:17:39 GMT

Very impressive. thanks a lot.