(*1*) # Generated by macrovsa_experiments.C do NOT edit
(*2*) with(Statistics):with(plots):
(*3*) X := Vector([100, 400, 1024, 2500, 4096, 10000, 1000000, 25000000], datatype=float):
(*4*) Y := Vector([0.004445, 0.0157525, 0.0639475, 0.236397, 0.495985, 1.87314, 2089.46, 433039], datatype=float):
(*5*) Xu := Vector([0.004445, 6.95311e-310, 0.0639475, 6.95311e-310, 0.495985, 1.87314, 6.95311e-310], datatype=float):
(*6*) Yu := Vector([0.00095, 0.00079, 0.001645, 0.003435, 0.00507, 0.02295, 1.78942], datatype=float):
(*7*) fit1 := Fit(a+b*d^c, X, Y, d, initialvalues = [a = 0.05, b = 0.00003, c = 1.35], output = [leastsquaresfunction, residualstandarddeviation]);
                                                   -6  1.65710837903622
fit1 := [0.219028870290985 + 0.238417656900174 x 10   d                ,

    0.3289928865]

(*8*) f1 := unapply(fit1[1], d):
(*9*) fit2 := Fit(a+b*d^1.5, X, Y, d, output=[leastsquaresfunction, residualstandarddeviation], summarize=true);
Summary:
----------------
Model: -195.59430+.34658012e-5*d^1.5
----------------
Coefficients:
    Estimate    Std. Error    t-value    P(>|t|)
a   -195.5943    196.7979     -0.9939     0.3587
b    0.0000      0.0000        778.3269   0.0000
----------------
R-squared: 1.0000, Adjusted R-squared: 1.0000
                                                    -5  1.5
 fit2 := [-195.594301310103 + 0.346580120438532 x 10   d   , 520.082179395601]

(*10*) f2 := unapply(fit2[1], d):
(*11*) plotsetup(jpeg, plotoutput="../public/macrovsa_experiments/durations_interpolation.jpg", plotoptions="width=600,height=600"):
(*12*) display(plot(f1(d), d=100..10000), pointplot(X, Y), title="binding computation time in msec");
(*13*) fit3 := Fit(a+b*d, Xu, Yu, d, initialvalues = [a = 0.00, b = 0.00003], output = [leastsquaresfunction, residualstandarddeviation]);
      fit3 := [0.331164898251176 - 0.202621842820917 d, 0.722177763734505]

(*14*) f3 := unapply(fit3[1], d):
(*15*) plotsetup(jpeg, plotoutput="../public/macrovsa_experiments/durations_interpolation_2.jpg", plotoptions="width=600,height=600"):
(*16*) display(plot(f3(d), d=100..10000), pointplot(Xu, Yu), title="bundling computation time in msec");
(*17*) quit: