Compilation without MainEvaluate

Question

I'm trying to avoid MainEvaluate in CompiledFunction. I've tried below but I couldn't get rid of it. I'd appreciate any help!!

Here f1 and f2 will be rank 1 list with real type.

decision[f1_, f2_] := And @@ Thread[LessEqual[f1, f2]] && Or @@ Thread[Less[f1, f2]];

dominates = Compile[{{f1, _Real, 1}, {f2, _Real, 1}}, 
decision[f1, f2], {{decision[_, _], True | False}}, 
CompilationOptions ->{"InlineCompiledFunctions" -> True}];
<< CompiledFunctionTools`
CompilePrint[dominates]

The expression in decision[] is not compilable. Try rewriting it using UnitStep[]. — Michael E2, Oct 12 '16 at 23:40
So you want all the elements of f1 less than or equal to the corresponding element of f2 and at least one the pairs unequal? — Michael E2, Oct 12 '16 at 23:42

score 1 · Accepted Answer · edited Apr 13 '17 at 12:55

1

Perhaps this?

dominates = Compile[{{f1, _Real, 1}, {f2, _Real, 1}}, 
  Length[Cases[UnitStep[f2 - f1], 0, 1]] == 0 && 
   Length[Cases[Unitize[f2 - f1], 1, 1]] > 0
  ]

For further discussion, see para. 2.3 of this answer: Performance tuning in Mathematica?

edited Apr 13 '17 at 12:55

Community

1

answered Oct 12 '16 at 23:50

Michael E2

235,386
17
334
747

Compilation without MainEvaluate

1 Answers1