Difference between pages "Dartboard PI" and "Asin"

From Mesham
(Difference between pages)
Jump to navigationJump to search
(Download)
 
 
Line 1: Line 1:
<metadesc>Mesham is a type oriented programming language allowing the writing of high performance parallel codes which are efficient yet simple to write and maintain</metadesc>
 
 
== Overview ==
 
== Overview ==
  
[[Image:dartboard.jpg|thumb|260px|right|Dartboard method to find PI]]
+
The asin(d) function will find the inverse sine of the value or variable ''d'' passed to it.
The dartboard method to find PI is a simple algorithm to find the value of PI. At this point it must be noted that there are much better methods out there to find PI, however the dartboard method is embarasingly parallel and as such quite simple to parallelise. The basic premise is that you can throw n darts randomly at a round dartboard on a square backing. As each dart is thrown randomly the ratio of darts hitting the board to those landing on the square is equal to the ratio between the two areas - which is PI / 4. Of course, the more darts you simulate throwing at the board, the better the approximation of PI - in our example each process will perform this simulated throwing a number of times, and then each process's approximation of PI is combined and averaged by one of the processes to obtain the result. Very ruffly, this means that with d darts, thrown over r rounds on n processes, the time taken parallely is the time it takes to simulate throwing d * r darts, yet a sequential algorithm would need to simulate throwing d * r * n darts. (We have excluded the consideration of communication costs from the parallel situation to simplify the concept.) Hopefully quite obviously, in the example by changing the number of processes, the number of rounds and the number of darts to throw in each round will directly change the accuracy of the result.
 
  
== Source Code ==
+
* '''Pass:''' A [[Double]] to find the inverse sine of
 +
 
 +
* '''Returns:''' A [[Double]] representing the inverse sine
 +
 
 +
== Example ==
  
 
  #include <maths>
 
  #include <maths>
#include <io>
 
#include <string>
 
 
var m:=64; // number of processes
 
 
   
 
   
 
  function void main() {
 
  function void main() {
     var calculatedPi:array[Double,m]:: allocated[single[on[0]]];
+
     var d:=asin(23);
    var mypi:Double;
+
     var y:=asin(d);
    var p;
 
    par p from 0 to m - 1 {
 
      var darts:=10000; // number of darts to simulate throwing each round
 
      var rounds:=100; // number of rounds of darts to throw
 
      var i;
 
      for i from 0 to rounds - 1 {
 
          mypi:=mypi + (4.0 * (throwdarts(darts) / darts));
 
      };
 
      mypi:=mypi / rounds;
 
      calculatedPi[p]:=mypi;
 
    };
 
    sync;
 
     proc 0 {
 
      var avepi:Double;
 
      var i;
 
      for i from 0 to m - 1 {   
 
          avepi:=avepi + calculatedPi[i];
 
      };
 
      avepi:=avepi / m;
 
      print(dtostring(avepi, "%.2f")+"\n");
 
    };
 
};
 
 
function Int throwdarts(var darts:Int)
 
{
 
    var score:=0;
 
    var n:=0;
 
    for n from 0 to darts - 1 {
 
      var xcoord:=randomnumber(0,1);
 
      var ycoord:=randomnumber(0,1);
 
      if ((pow(xcoord,2) + pow(ycoord,2)) < 1.0) {
 
          score++; // hit the dartboard!
 
      };
 
    };
 
    return score;
 
 
  };
 
  };
  
''This code requires at least Mesham version 1.0''
+
''Since: Version 1.0''
 
 
== Notes ==
 
 
 
An interesting aside is that we have used a function in this example, yet there is no main function. The throwdarts function will simulate throwing the darts for each round. As already noted in the language documentation, the main function is optional and without it the compiler will set the program entry point to be the start of the source code.
 
== Download ==
 
 
 
The dartboard method to compute PI source code is located [http://www.mesham.com/downloads/pi.mesh here] a legacy version for Mesham 0.5 can be downloaded [http://www.mesham.com/downloads/pi-0.5.mesh here]
 
  
[[Category:Example Codes]]
+
[[Category:Function Library]]
 +
[[Category:Maths Functions]]

Revision as of 14:32, 15 April 2013

Overview

The asin(d) function will find the inverse sine of the value or variable d passed to it.

  • Pass: A Double to find the inverse sine of
  • Returns: A Double representing the inverse sine

Example

#include <maths>

function void main() {
   var d:=asin(23);
   var y:=asin(d);
};

Since: Version 1.0