algorithm - Can a Monte Carlo pi calculation be used for a world record? -


i have random function calculate pi monte carlo style:

enter image description here

max=10000000; format long;  in = 0; tic k=1:max     x = rand();     y = rand();     if sqrt(x^2 + y^2) < 1         in = in + 1;     end end  toc calc_pi = 4*in/max  epsilon = abs(pi - calc_pi)  calcpi(100000000);

if iterate 10e100 times, algorithm compete the world record? if so, how can find number of iteration give nth digit?

this nice exercise calculating pi, inefficient one. remarks:

  • my statistics rusty before had coffee, guess error scales 1 / sqrt(n_guess). n digits, need error of 10^(-n), need (10^n)^2 random guesses. if 1e100 guesses, proposed, on order of 50 digits of pi! number of iteration exponentional function of number of required digits, horribly slow. algorithm maybe linear in number of digits want.

  • with large number of guesses required, have start questioning quality of random number generator.

  • your algorithm limited floating-point errors 1e-16 or so. calculating digits of pi requires sort of arbitrary precision number format.

  • to speed algorithm, can leave out sqrt().

  • don't use variable called max, overwrites existing function. use n_guess or so.

quick , dirty test prove theory (after coffee):

pie = @(n) 4 * nnz(rand(n,1).^2 + rand(n, 1).^2 < 1) / n; ntrial = round(logspace(1, 8, 100)); pies = arrayfun(pie, ntrial);  loglog(ntrial, abs(pies - pi), '.k', ntrial, ntrial.^-.5, '--r') xlabel('ntrials') ylabel('epsilon') legend('monte carlo', '1 / sqrt(ntrial)')

monte carlo result


-

Comments

Post a Comment

Popular posts from this blog

scala - 'wrong top statement declaration' when using slick in IntelliJ -

Image
i'm porting project scala 2.9 2.10, therefore have use slick instead of scalaquery. i'm using slick 2.1.0 since supports ms access. according this tutorial , upgrade guide changed robs object class , added val robs : why error message wrong top statement declaration , how rid of it? edit: i'm new scala... according this question seems can't put val outside of methods or classes, right? code above directly in packge. right approach slick then? should move code in class or trait? according this changed val robs = tablequery[robs] to object robs extends tablequery(new robs(_)) {} . no warnings or errors anymore. :)
Read more

c# - DevExpress.Wpf.Grid.InfiniteGridSizeException was unhandled -

Image
when using devexpress, see error: devexpress.wpf.grid.infinitegridsizeexception unhandled message="by default, infinite grid height not allowed since grid rows rendered , hence grid work slowly. fix issue, should place grid container give finite height grid, or should manually specify grid's height or maxheight. note can avoid exception setting gridcontrol.allowinfinitegridsize static property true, in case grid run slowly." the problem dxgrid has infinite height. to fix, set height non-infinite. snoop absolutely invaluable this: if "height" xaml element infinite (i.e. 0 or nan ), can set using one of following: option 1: height="{binding path=actualheight, relativesource={relativesource mode=findancestor, ancestortype=uielement}}" option 2: verticalalignment="stretch" option 3: height="auto" hint: use verticalalignment="stretch" if child of grid <rowdefinition height=
Read more

PySide and Qt Properties: Connecting signals from Python to QML -

i trying run simple pyside appliccation design made in qml. call function in qml python in separated thread. i have following qml view: import qt 4.7 rectangle { visible: true width: 100; height: 100 color: "lightgrey" text { id: datestr objectname: "datestr" x: 10 y: 10 text: "init text" font.pixelsize: 12 function updatedata(text) { datestr.text = qstr(text) console.log("you said: " + text) } } } here application: some imports.... #!/usr/bin/python import sys import time pyside import qtdeclarative, qtcore pyside.qtgui import qapplication pyside.qtdeclarative import qdeclarativeview the qthread worker.... #worker thread class worker(qtcore.qthread): updateprogress = qtcore.signal(int) def __init__(self, child): qtcore.qthread.__init__(self) self.child = child def run(self):
Read more