Mathematica：編譯NSolve，NMinimize和其他數值函數

Question

有沒有辦法使用Compile []編譯這些函數？看來Mathematica總是失敗。

我在NSolve和NMinimize中有純粹的數值函數，但它仍然不起作用。

一個小小的例子是：

tempfunc = Compile[{}, NSolve[x^2 - 5 x + 6 == 0, x]] 
tempfunc[] 

這會給出一個錯誤......

不知道如何編寫這些功能，讓他們更快？

Answer 1

看來這是不可能的。評估：

Compile`CompilerFunctions[] // Sort 

這將給出「可編譯」功能列表。 Mathematica的10.3它會給：

{Abs, AddTo, And, Append, AppendTo, Apply, ArcCos, ArcCosh, ArcCot, ArcCoth, 
ArcCsc, ArcCsch, ArcSec, ArcSech, ArcSin, ArcSinh, ArcTan, ArcTanh, Arg, 
Array, ArrayDepth, Internal`Bag, Internal`BagPart, BitAnd, BitNot, BitOr, 
BitXor, Block, BlockRandom, Boole, Break, Cases, Catch, Ceiling, Chop, 
Internal`CompileError, System`Private`CompileSymbol, Complement, 
ComposeList, CompoundExpression, Conjugate, ConjugateTranspose, Continue, 
Cos, Cosh, Cot, Coth, Count, Csc, Csch, Decrement, Delete, DeleteCases, 
Dimensions, Divide, DivideBy, Do, Dot, Drop, Equal, Erf, Erfc, EvenQ, Exp, 
Fibonacci, First, FixedPoint, FixedPointList, Flatten, 
NDSolve`FEM`FlattenAll, Floor, Fold, FoldList, For, FractionalPart, FreeQ, 
Gamma, Compile`GetElement, Goto, Greater, GreaterEqual, Gudermannian, 
Haversine, If, Im, Implies, Increment, Indexed, Inequality, Compile`InnerDo, 
Insert, IntegerDigits, IntegerPart, Intersection, InverseGudermannian, 
InverseHaversine, Compile`IteratorCount, Join, Label, Last, Length, Less, 
LessEqual, List, Log, Log10, Log2, LogGamma, LogisticSigmoid, LucasL, Map, 
MapAll, MapAt, MapIndexed, MapThread, NDSolve`FEM`MapThreadDot, MatrixQ, 
Max, MemberQ, Min, Minus, Mod, Compile`Mod1, Module, Most, N, Negative, 
Nest, NestList, NonNegative, Not, OddQ, Or, OrderedQ, Out, Outer, Part, 
Partition, Piecewise, Plus, Position, Positive, Power, PreDecrement, 
PreIncrement, Prepend, PrependTo, Product, Quotient, Random, RandomChoice, 
RandomComplex, RandomInteger, RandomReal, RandomSample, RandomVariate, 
Range, Re, Internal`ReciprocalSqrt, ReplacePart, Rest, Return, Reverse, 
RotateLeft, RotateRight, Round, RuleCondition, SameQ, Scan, Sec, Sech, 
SeedRandom, Select, Set, SetDelayed, Compile`SetIterate, Sign, Sin, Sinc, 
Sinh, Sort, Sqrt, Internal`Square, Internal`StuffBag, Subtract, 
SubtractFrom, Sum, Switch, Table, Take, Tan, Tanh, TensorRank, Throw, Times, 
TimesBy, Tr, Transpose, Unequal, Union, Unitize, UnitStep, UnsameQ, VectorQ, 
Which, While, With, Xor} 

我試圖

cnm = Compile[{{a, _Real}}, 
    Block[{v, r}, 
    {v, r} = NMinimize[2 x^2 + x - a, x]; 
    {v, x /. r} 
    ] 
    ] 

和結果的非編數學代碼評價：

In[35]:= cnm[2] 

During evaluation of In[35]:= CompiledFunction::cfse: Compiled expression {-2.125,{x->-0.25}} should be a machine-size real number. >> 

During evaluation of In[35]:= CompiledFunction::cfex: Could not complete external evaluation at instruction 1; proceeding with uncompiled evaluation. >> 

Out[35]= {-2.125, -0.25}