[ADMB Users] Fwd: > dvariable in separable function

[Weihai Liu]

---------- Forwarded message ----------
From: Weihai Liu <liuweih at msu.edu>
Date: Tue, Nov 1, 2011 at 9:12 AM
Subject: Re: [ADMB Users] > dvariable in separable function
To: Mollie Brooks <mbrooks at ufl.edu>


for your second try, use these inside your separable function,
fpen=dvariable(0.0);
if(value(fpen)>0.0)

if it works, let me know.

weihai

On Mon, Oct 31, 2011 at 11:42 PM, Mollie Brooks <mbrooks at ufl.edu> wrote:

> Hello,
> I need to use posfun in a separable function. Then I'd like to output the
> value of fpen if it's more than 0. I declared fpen as a dvariable within
> the separable function. There doesn't seem to be the right operator defined.
> sepvar_shts_nb.cpp:243: error: no match for ‘operator>’ in ‘fpen > 0.0’
> /home/mbrooks/admb/include/fvar.hpp:7923: note: candidates are: int
> operator>(double, const prevariable&)
> /home/mbrooks/admb/include/df32fun.h:290: note:                 int
> operator>(const df3_two_variable&, double)
> /home/mbrooks/admb/include/df32fun.h:296: note:                 int
> operator>(const df3_two_variable&, const df3_two_variable&)
>
>
> Instead, I tried declaring it as a number in the parameter section and
> sending it to the separable function, but I have to send it as a dvariable
> so I have the same problem.
>
> SEPARABLE_FUNCTION void obs(dvariable fpen,int i,const dvar_vector&
> fcoeffs,const dvariable& indiv_std,const dvariable& indiv_slope_std,const
> dvariable& year_std,const dvariable& plot_std,const dvariable&
> rcoeff_indiv,const dvariable& rcoeff_indiv_slope,const dvariable&
> rcoeff_plot,const dvariable& rcoeff_year,const dvariable& k)
>        dvariable prediction=fixmm(i)*fcoeffs+
>                                indiv_std*rcoeff_indiv+
>
>  indiv_slope_std*rcoeff_indiv_slope*fixmm(i,1)+ //The 1st column should be
> shts
>                                year_std*rcoeff_year+
>                                plot_std*rcoeff_plot;
>        fpen = 0.0;
>        prediction=posfun(prediction, 0.0000001, fpen);
>        if(fpen>0.0)cout<<endl<<"fpen="<<fpen<<endl;
>        jnll+=100000*fpen;
>
>  jnll-=gammln(k+shts_next(i))-gammln(k)-gammln(shts_next(i)+1)+k*log(k)-k*log(prediction+k)+shts_next(i)*log(prediction)-shts_next(i)*log(prediction+k);
>
> sepvar_shts_nb.cpp:245: error: no match for ‘operator=’ in ‘fpen = 0.0’
> /home/mbrooks/admb/include/df1b2fnl.h:114: note: candidates are:
> funnel_init_df1b2variable& funnel_init_df1b2variable::operator=(const
> funnel_init_df1b2variable&)
> sepvar_shts_nb.cpp:247: error: no match for ‘operator>’ in ‘fpen > 0.0’
> /home/mbrooks/admb/include/fvar.hpp:7923: note: candidates are: int
> operator>(double, const prevariable&)
> /home/mbrooks/admb/include/df32fun.h:290: note:                 int
> operator>(const df3_two_variable&, double)
> /home/mbrooks/admb/include/df32fun.h:296: note:                 int
> operator>(const df3_two_variable&, const df3_two_variable&)
>
> Has anyone had this problem with using the > operator? Any ideas for
> workarounds?
> thanks,
> Mollie
>
> Mollie Brooks
> Ph.D. Candidate
> NSF IGERT Fellow
> Biology Department
> University of Florida
> mbrooks at ufl.edu
> http://people.biology.ufl.edu/mbrooks
>
>
>
>
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20111101/50f962cd/attachment.html>