anvaka · April 22, 2013 21:16
diff --git a/nbody_plain_objects.js b/nbody_plain_objects.js
 function runtIt(n) {
   'use strict';
   /* The Computer Language Benchmarks Game
      http://benchmarksgame.alioth.debian.org/
      contributed by Isaac Gouy 
      Optimized by Roy Williams
      "Deoptimized" to plain JS objects by Anvaka*/

   /**
    * @type {number}
    */
   var PI = 3.141592653589793;

   /**
    * @type {number}
    */
   var SOLAR_MASS = 4 * PI * PI;

   /**
    * @type {number}
    */
   var DAYS_PER_YEAR = 365.24;

   /**
    * @param {number} x
    * @param {number} y
    * @param {number} z
    * @param {number} vx
    * @param {number} vy
    * @param {number} vz
    * @param {number} mass
    * @constructor
    */
   function Body(x,y,z,vx,vy,vz,mass) {
      this.X = x;
      this.Y = y;
      this.Z = z;
      this.VX = vx;
      this.VY = vy;
      this.VZ = vz;
      this.MASS = mass;
   }



   /**
    * @param {number} px
    * @param {number} py
    * @param {number} pz
    */
   Body.prototype.offsetMomentum = function(px,py,pz) {
      this.VX = -px / SOLAR_MASS;
      this.VY = -py / SOLAR_MASS;
      this.VZ = -pz / SOLAR_MASS;
   };

   /**
    * @param {Array.<Body>} bodies
    * @constructor
    */
   function NBodySystem(bodies){
      /**
       * @type {Array.<Body>}
       */
      this.bodies = bodies;
      var px = 0.0;
      var py = 0.0;
      var pz = 0.0;
      var size = this.bodies.length;
      for (var i=0; i<size; i++){
         var b = this.bodies[i];
         var m = b.MASS;
         px += b.VX * m;
         py += b.VY * m;
         pz += b.VZ * m;
      }
      this.bodies[0].offsetMomentum(px,py,pz);
   }

   /**
    * @param {number} dt
    */
   NBodySystem.prototype.advance = function(dt){
      var dx, dy, dz, distance, mag;
      var size = this.bodies.length;

      for (var i=0; i<size; i++) {
         var bodyi = this.bodies[i];
         var imass = bodyi.MASS;
         for (var j=i+1; j<size; j++) {
            var bodyj = this.bodies[j];
            var jmass = bodyj.MASS;
            dx = bodyi.X - bodyj.X;
            dy = bodyi.Y - bodyj.Y;
            dz = bodyi.Z - bodyj.Z;

            distance = Math.sqrt(dx*dx + dy*dy + dz*dz);
            mag = dt / (distance * distance * distance);

            bodyi.VX -= dx * jmass * mag;
            bodyi.VY -= dy * jmass * mag;
            bodyi.VZ -= dz * jmass * mag;

            bodyj.VX += dx * imass * mag;
            bodyj.VY += dy * imass * mag;
            bodyj.VZ += dz * imass * mag;
         }
         bodyi.X += dt * bodyi.VX;
         bodyi.Y += dt * bodyi.VY;
         bodyi.Z += dt * bodyi.VZ;
      }
   };

   /**
    * @return {number}
    */
   NBodySystem.prototype.energy = function(){
      var dx, dy, dz, distance;
      var e = 0.0;
      var size = this.bodies.length;

      for (var i=0; i<size; i++) {
         var bodyi = this.bodies[i];

         e += 0.5 * bodyi.MASS *
            ( bodyi.VX * bodyi.VX + bodyi.VY * bodyi.VY + bodyi.VZ * bodyi.VZ );

         for (var j=i+1; j<size; j++) {
            var bodyj = this.bodies[j];
            dx = bodyi.X - bodyj.X;
            dy = bodyi.Y - bodyj.Y;
            dz = bodyi.Z - bodyj.Z;

            distance = Math.sqrt(dx*dx + dy*dy + dz*dz);
            e -= (bodyi.MASS *
                  bodyj.MASS) / distance;
         }
      }
      return e;
   };


    var sun = new Body(0.0, 0.0, 0.0, 0.0, 0.0, 0.0, SOLAR_MASS),
        jupiter = new Body(
            4.84143144246472090e+00,
            -1.16032004402742839e+00,
            -1.03622044471123109e-01,
            1.66007664274403694e-03 * DAYS_PER_YEAR,
            7.69901118419740425e-03 * DAYS_PER_YEAR,
            -6.90460016972063023e-05 * DAYS_PER_YEAR,
            9.54791938424326609e-04 * SOLAR_MASS
         ),
         saturn = new Body(
            8.34336671824457987e+00,
            4.12479856412430479e+00,
            -4.03523417114321381e-01,
            -2.76742510726862411e-03 * DAYS_PER_YEAR,
            4.99852801234917238e-03 * DAYS_PER_YEAR,
            2.30417297573763929e-05 * DAYS_PER_YEAR,
            2.85885980666130812e-04 * SOLAR_MASS
         ),

         uranus = new Body(
            1.28943695621391310e+01,
            -1.51111514016986312e+01,
            -2.23307578892655734e-01,
            2.96460137564761618e-03 * DAYS_PER_YEAR,
            2.37847173959480950e-03 * DAYS_PER_YEAR,
            -2.96589568540237556e-05 * DAYS_PER_YEAR,
            4.36624404335156298e-05 * SOLAR_MASS
         ),

         neptune = new Body(
            1.53796971148509165e+01,
            -2.59193146099879641e+01,
            1.79258772950371181e-01,
            2.68067772490389322e-03 * DAYS_PER_YEAR,
            1.62824170038242295e-03 * DAYS_PER_YEAR,
            -9.51592254519715870e-05 * DAYS_PER_YEAR,
            5.15138902046611451e-05 * SOLAR_MASS
         );

     var bodies = new NBodySystem(new Array(sun, jupiter, saturn, uranus, neptune));
     print(bodies.energy().toFixed(9));
     for (var i=0; i<n; i++){ bodies.advance(0.01); }
     print(bodies.energy().toFixed(9));
 }

 /**
 * @type {number}
 */
 var n = arguments[0];
 runtIt(n);
	function runtIt(n) {
	'use strict';
	/* The Computer Language Benchmarks Game
	http://benchmarksgame.alioth.debian.org/
	contributed by Isaac Gouy
	Optimized by Roy Williams
	"Deoptimized" to plain JS objects by Anvaka*/

	/**
	* @type {number}
	*/
	var PI = 3.141592653589793;

	/**
	* @type {number}
	*/
	var SOLAR_MASS = 4 * PI * PI;

	/**
	* @type {number}
	*/
	var DAYS_PER_YEAR = 365.24;

	/**
	* @param {number} x
	* @param {number} y
	* @param {number} z
	* @param {number} vx
	* @param {number} vy
	* @param {number} vz
	* @param {number} mass
	* @constructor
	*/
	function Body(x,y,z,vx,vy,vz,mass) {
	this.X = x;
	this.Y = y;
	this.Z = z;
	this.VX = vx;
	this.VY = vy;
	this.VZ = vz;
	this.MASS = mass;
	}



	/**
	* @param {number} px
	* @param {number} py
	* @param {number} pz
	*/
	Body.prototype.offsetMomentum = function(px,py,pz) {
	this.VX = -px / SOLAR_MASS;
	this.VY = -py / SOLAR_MASS;
	this.VZ = -pz / SOLAR_MASS;
	};

	/**
	* @param {Array.<Body>} bodies
	* @constructor
	*/
	function NBodySystem(bodies){
	/**
	* @type {Array.<Body>}
	*/
	this.bodies = bodies;
	var px = 0.0;
	var py = 0.0;
	var pz = 0.0;
	var size = this.bodies.length;
	for (var i=0; i<size; i++){
	var b = this.bodies[i];
	var m = b.MASS;
	px += b.VX * m;
	py += b.VY * m;
	pz += b.VZ * m;
	}
	this.bodies[0].offsetMomentum(px,py,pz);
	}

	/**
	* @param {number} dt
	*/
	NBodySystem.prototype.advance = function(dt){
	var dx, dy, dz, distance, mag;
	var size = this.bodies.length;

	for (var i=0; i<size; i++) {
	var bodyi = this.bodies[i];
	var imass = bodyi.MASS;
	for (var j=i+1; j<size; j++) {
	var bodyj = this.bodies[j];
	var jmass = bodyj.MASS;
	dx = bodyi.X - bodyj.X;
	dy = bodyi.Y - bodyj.Y;
	dz = bodyi.Z - bodyj.Z;

	distance = Math.sqrt(dxdx + dydy + dz*dz);
	mag = dt / (distance * distance * distance);

	bodyi.VX -= dx * jmass * mag;
	bodyi.VY -= dy * jmass * mag;
	bodyi.VZ -= dz * jmass * mag;

	bodyj.VX += dx * imass * mag;
	bodyj.VY += dy * imass * mag;
	bodyj.VZ += dz * imass * mag;
	}
	bodyi.X += dt * bodyi.VX;
	bodyi.Y += dt * bodyi.VY;
	bodyi.Z += dt * bodyi.VZ;
	}
	};

	/**
	* @return {number}
	*/
	NBodySystem.prototype.energy = function(){
	var dx, dy, dz, distance;
	var e = 0.0;
	var size = this.bodies.length;

	for (var i=0; i<size; i++) {
	var bodyi = this.bodies[i];

	e += 0.5 * bodyi.MASS *
	( bodyi.VX * bodyi.VX + bodyi.VY * bodyi.VY + bodyi.VZ * bodyi.VZ );

	for (var j=i+1; j<size; j++) {
	var bodyj = this.bodies[j];
	dx = bodyi.X - bodyj.X;
	dy = bodyi.Y - bodyj.Y;
	dz = bodyi.Z - bodyj.Z;

	distance = Math.sqrt(dxdx + dydy + dz*dz);
	e -= (bodyi.MASS *
	bodyj.MASS) / distance;
	}
	}
	return e;
	};


	var sun = new Body(0.0, 0.0, 0.0, 0.0, 0.0, 0.0, SOLAR_MASS),
	jupiter = new Body(
	4.84143144246472090e+00,
	-1.16032004402742839e+00,
	-1.03622044471123109e-01,
	1.66007664274403694e-03 * DAYS_PER_YEAR,
	7.69901118419740425e-03 * DAYS_PER_YEAR,
	-6.90460016972063023e-05 * DAYS_PER_YEAR,
	9.54791938424326609e-04 * SOLAR_MASS
	),
	saturn = new Body(
	8.34336671824457987e+00,
	4.12479856412430479e+00,
	-4.03523417114321381e-01,
	-2.76742510726862411e-03 * DAYS_PER_YEAR,
	4.99852801234917238e-03 * DAYS_PER_YEAR,
	2.30417297573763929e-05 * DAYS_PER_YEAR,
	2.85885980666130812e-04 * SOLAR_MASS
	),

	uranus = new Body(
	1.28943695621391310e+01,
	-1.51111514016986312e+01,
	-2.23307578892655734e-01,
	2.96460137564761618e-03 * DAYS_PER_YEAR,
	2.37847173959480950e-03 * DAYS_PER_YEAR,
	-2.96589568540237556e-05 * DAYS_PER_YEAR,
	4.36624404335156298e-05 * SOLAR_MASS
	),

	neptune = new Body(
	1.53796971148509165e+01,
	-2.59193146099879641e+01,
	1.79258772950371181e-01,
	2.68067772490389322e-03 * DAYS_PER_YEAR,
	1.62824170038242295e-03 * DAYS_PER_YEAR,
	-9.51592254519715870e-05 * DAYS_PER_YEAR,
	5.15138902046611451e-05 * SOLAR_MASS
	);

	var bodies = new NBodySystem(new Array(sun, jupiter, saturn, uranus, neptune));
	print(bodies.energy().toFixed(9));
	for (var i=0; i<n; i++){ bodies.advance(0.01); }
	print(bodies.energy().toFixed(9));
	}

	/**
	* @type {number}
	*/
	var n = arguments[0];
	runtIt(n);