如何在函數參數中使用元編程？

Question

它是我第二天與朱莉婭學習和實驗。雖然我仔細閱讀了關於元編程的文檔（但也許不夠小心）和幾個相同的線程。我仍然無法弄清楚如何在函數內部使用它。 我tryed做一些數據更靈活的模擬如下功能：

using Distributions 
function gendata(N,NLATENT,NITEMS) 
    latent = repeat(rand(Normal(6,2),N,NLATENT), inner=(1,NITEMS)) 
    errors = rand(Normal(0,1),N,NLATENT*NITEMS) 
    x = latent+errors 
end 

通過這樣做：

using Distributions 
function gendata(N,NLATENT,NITEMS,LATENT_DIST="Normal(0,1)",ERRORS_DIST="Normal(0,1)") 
    to_eval_latent = parse("latent = repeat(rand($LATENT_DIST,N,NLATENT), inner=(1,NITEMS))") 
    eval(to_eval_latent) 
    to_eval_errors = parse("error = rand($ERRORS_DIST,N,NLATENT*NITEMS)") 
    eval(to_eval_errors) 
    x = latent+errors 
end 

但由於EVAL不要在局部範圍內工作，它不工作。我能做些什麼來圍繞這個工作？

也是原本的功能，似乎不是那麼快，我在性能上做了什麼大的失誤？

我真的appriciate任何recommandation。 在此先感謝。

Answer 1

不需要在那裏使用eval，通過將分佈類型作爲關鍵字args（或帶有默認值的命名參數）傳遞，您可以保持相同的靈活性。解析和評估「字符串型」參數通常會使優化失敗，應該避免。

function gendata(N,NLATENT,NITEMS; LATENT_DIST=Normal(0,1),ERRORS_DIST=Normal(0,1)) 
     latent = repeat(rand(LATENT_DIST,N,NLATENT), inner=(1,NITEMS)) 
     errors = rand(ERRORS_DIST,N,NLATENT*NITEMS) 
     x = latent+errors 
end 


julia> gendata(10,2,3, LATENT_DIST=Pareto(.3)) 
... 


julia> gendata(10,2,3, ERRORS_DIST=Gamma(.6)) 
... 

等

Answer 2

你不是真的要在這裏使用eval（慢，不會產生類型信息，將與編譯干擾等），但如果你想了解什麼出了問題，這裏是你會怎麼做：

無論它從代碼的其餘部分分開：

function gendata(N,NLATENT,NITEMS,LDIST_EX="Normal(0,1)",EDIST_EX="Normal(0,1)") 

    # Eval your expressions separately 

    LATENT_DIST = eval(parse(LDIST_EX)) 
    ERRORS_DIST = eval(parse(EDIST_EX)) 

    # Do your thing 

    latent = repeat(rand(LATENT_DIST,N,NLATENT), inner=(1,NITEMS)) 
    errors = rand(ERROR_DIST,N,NLATENT*NITEMS) 
    x = latent+errors  
end 

或者使用插補引用表達式：

function gendata(N,NLATENT,NITEMS,LDIST_EX="Normal(0,1)",EDIST_EX="Normal(0,1)") 

    # Obtain expression objects 

    LATENT_DIST = parse(LDIST_EX) 
    ERRORS_DIST = parse(EDIST_EX) 

    # Eval but interpolate in everything that's local to the function 
    # And you can't introduce local variables with eval so keep them 
    # out of it. 

    latent = eval(:(repeat(rand($LATENT_DIST,$N,$NLATENT), inner=(1,$NITEMS)))) 
    errors = eval(:(rand($ERRORS_DIST, $N, $NLATENT*$NITEMS))) 
    x = latent+errors 
end 

你也可以使用一個單一的eval使用let塊引入一個自包含的範圍：

function gendata(N,NLATENT,NITEMS,LDIST_EX="Normal(0,1)",EDIST_EX="Normal(0,1)") 

    LATENT_DIST = parse(LDIST_EX) 
    ERRORS_DIST = parse(EDIST_EX) 
    x = 
    @eval let 
    latent = repeat(rand($LATENT_DIST,$N,$NLATENT), inner=(1,$NITEMS)) 
    errors = (rand($ERRORS_DIST, $N, $NLATENT*$NITEMS)) 
    latent+errors 
    end 
end 

((@eval x) == eval(:(x)))

好吧，希望你能理解eval事情更好一點。第二天，我的意思是，你應該試驗;）