From: Joachim Schimpf <joachim.schimpf_at_infotech.monash.edu.au>

Date: Wed, 16 Sep 2009 11:43:01 +1000

Date: Wed, 16 Sep 2009 11:43:01 +1000

liu.yu.senior.china wrote: > Dear All, > > I have a constraint with exponentiation as follows: > > :- lib(ic). > dv(Fs):- > Fs::[18000,20000,22000,24000], > Fmax is 24000, > Fmin is 18000, > Lamda0 is 1e-6, > D is 4, > Si is 1, > Lamda = Lamda0*10^(D*(Fmax-Fs)/(Fmax-Fmin)), > Rfmax = exp(-Lamda*1000000/Fmax/1e5), > Q = (1-Rfmax)^(Si-1), > eval(Q) $>= -1000. > > Since Rfmax is close to 1 (0.99..) and Si is 1, Q will be 1. So the > constraint is actually redundent and does not have any effect on Fs. > However, I get a No as a result. Anything wrong? Yes, there seems to be a bug in the implementation of exponentiation that causes failure when the exponent is 0. I've recorded this at http://eclipse-clp.org/bugzilla/show_bug.cgi?id=674 > > My second question is aobut eval/1. By using $>=/2, the goals are > suspended before all variables are instantiated. Then why do I have to > use eval? I get errors without it. > > If I don't use Q to pass the expresson and use eval((1-Rfmax)^(Si-1)) > $>= -1000, I get the same error as I use Q and do not use eval. The > error is "some constants are expected". Eval is used to mark variables - which occur as arguments of constraints in your source code, and - which, at runtime, will not be numbers or variables, but symbolic expressions So if you write your constraint as (1-Rfmax)^(Si-1) $>= -1000 then you have to mark Rfmax, because its value is a symbolic expression. (1-eval(Rfmax))^(Si-1) $>= -1000 Si on the other hand is just a number and doesn't need marking. For efficiency, it is best however, to expand the constraint as much as possible in your source code, because then the system can preprocess it better at compile time, i.e. in your case (1-exp(-(Lamda0*10^(D*(Fmax-Fs)/(Fmax-Fmin)))*1000000/Fmax/1e5))^(Si-1) $>= -1000. -- JoachimReceived on Wed Sep 16 2009 - 02:43:20 CEST

*
This archive was generated by hypermail 2.2.0
: Thu Feb 02 2012 - 02:31:58 CET
*