Bug in compilefis function

[egormir2188]

The compilefis function from FuzzyLogic.jl does not always compile the fuzzy logic function correctly. This can be verified by plotting the response surface using the compiled and uncompiled function. The surfaces do not match in some cases.

[lucaferranti]

Hi, thank you! Could you provide an example?

[On-datra]

Hi, thank you! Could you provide an example?

Hello!
This can be verified by plotting the response surface using the compiled and uncompiled function. The surfaces do not match in some cases.

I tested this code on Julia 1.11.1

using Plots
using FuzzyLogic

S_fis = @sugfis function SugenoCtrl(E, delE)::U
    E := begin
        domain = -1.0:1.0
        N = TriangularMF(-2.0, -1.0, 0.0)
        Z = TriangularMF(-1.0, 0.0, 1.0)
        P = TriangularMF(0.0, 1.0, 2.0)
    end

    delE := begin
        domain = -1.0:1.0
        N = TriangularMF(-2.0, -1.0, 0.0)
        Z = TriangularMF(-1.0, 0.0, 1.0)
        P = TriangularMF(0.0, 1.0, 2.0)
    end    

    U := begin
        domain = -1.0:1.0
        NB = -1.0
        NM  = -0.5
        Z = 0.0
        PM  = 0.5
        PB = 1.0
    end

   E == N  && delE == N  -->  U == NB
   E == Z  && delE == N  -->  U == NM
   E == P  && delE == N  -->  U == Z
   E == N  && delE == Z  -->  U == NM
   E == Z  && delE == Z  -->  U == Z
   E == P  && delE == Z  -->  U == PM
   E == N  && delE == P  -->  U == Z
   E == Z  && delE == P  -->  U == PM
   E == P  && delE == P  -->  U == PB
end;
S_fis_ex = compilefis(S_fis)

##-------------------------------------------------------------##
## Next I get the SugenoCtrl function using compilefis
function SugenoCtrl(E, delE)
             N = max(min((E - -2.0) / 1.0, (0.0 - E) / 1.0), 0)
             Z = max(min((E - -1.0) / 1.0, (1.0 - E) / 1.0), 0)
             P = max(min((E - 0.0) / 1.0, (2.0 - E) / 1.0), 0)
             N = max(min((delE - -2.0) / 1.0, (0.0 - delE) / 1.0), 0)
             Z = max(min((delE - -1.0) / 1.0, (1.0 - delE) / 1.0), 0)
             P = max(min((delE - 0.0) / 1.0, (2.0 - delE) / 1.0), 0)
             ant1 = N * N
             ant2 = Z * N
             ant3 = P * N
             ant4 = N * Z
             ant5 = Z * Z
             ant6 = P * Z
             ant7 = N * P
             ant8 = Z * P
             ant9 = P * P
             tot_weight = ant1 + ant2 + ant3 + ant4 + ant5 + ant6 + ant7 + ant8 + ant9
             NB = -1.0
                   NM = -0.5
             Z = 0.0
             PM = 0.5
             PB = 1.0
             U = (ant1 * NB + ant2 * NM + ant3 * Z + ant4 * NM + ant5 * Z + ant6 * PM + ant7 * Z + ant8 * PM + ant9 * PB) / tot_weight
                   return U
                     end
E_range = -1.0:0.05:1.0 
delE_range = -1.0:0.05:1.0
fis_table = zeros(length(E_range), length(delE_range))
SugenoCtrl_table = zeros(length(E_range), length(delE_range))
i = 0;
for e in E_range
    i = i+1
    j = 0;
    for de in delE_range
        j = j+1
        # table by fuzzy logic function 
        fis_table[i, j] = S_fis([e, de])[:U]
        # compiled function table
        SugenoCtrl_table[i, j] = SugenoCtrl(e, de)       
    end
end
plot(layout=@layout([a b с]))
surface!(E_range, delE_range, fis_table, subplot=1,
       title="fuzzy function")
surface!(E_range, delE_range, SugenoCtrl_table, subplot=2,
       title="compiled function")   
gensurf!(S_fis, linewidth=3, subplot=3,
        title="with gensurf")