4

I would like to some programming that is very generic. Particularly I am interested in the following:

Let's say I want to write a function 

fuu[t_, f_] := Integrate[Exp[a*f[q]], {q, 0, 2*Pi}] /. a -> t

where t is a real value and f some function.

Now if I run fuu multiple times for different values of t, say 1,2,3 etc. and keep the function f as Cos, it will re-run the integration multiple times, which is undesirable.

Storing values is also possible if I write 

fuu[t_, f_] := fuu[t, f] = Integrate[Exp[a*f[q]], {q, 0, 2*Pi}] /. a -> t

The problem is that the results from the integration are only stored for specific values of t then. How do I store it for generic t but f fixed as Cos?

Sure one may argue why not work around this by some other techniques. But the goal in Mathematica is to have effortless code that is still fast - so does anyone have suggestions?

If someone can give me a solution. How can this solution be combined with Let's say that I would want to compile the result from Integrate to get something even better?

Roman
  • 47,322
  • 2
  • 55
  • 121
Michael
  • 123
  • 5
  • Have a look at this discussion. – Leonid Shifrin Jun 25 '19 at 21:52
  • Thank you for the link. This is very generous of you but this does not really solve my problem because different instances in the solutions provided have to be put by hand and are not dynamically determined based on what type of problem I keep throwing at mathematica. Nevertheless I very much appreciate your comment Leonid. – Michael Jun 25 '19 at 21:58

1 Answers1

9

partial memoization:

Clear[fuu];
fuu[f_] := fuu[f] =
  Function[t, Evaluate[Integrate[Exp[t*f[q]], {q, 0, 2*Pi}]]]

test:

fuu[Cos]
(*    Function[t$, 2 π BesselI[0, t$]]    *)

use this function for different values of t:

fuu[Cos][0.2]
(*    6.34617    *)

fuu[Cos][0.4]
(*    6.53704    *)
Roman
  • 47,322
  • 2
  • 55
  • 121