/*
* Licensed to the Apache Software Foundation (ASF) under one
* or more contributor license agreements. See the NOTICE file
* distributed with this work for additional information
* regarding copyright ownership. The ASF licenses this file
* to you under the Apache License, Version 2.0 (the
* "License"); you may not use this file except in compliance
* with the License. You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing,
* software distributed under the License is distributed on an
* "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
* KIND, either express or implied. See the License for the
* specific language governing permissions and limitations
* under the License.
*/
package org.codehaus.groovy.benchmarks.vm5;
class psum {
static final double twothirds = 2.0/3.0;
public static void main(String[] args){
long start = System.currentTimeMillis ();
int n = Integer.parseInt(args[0]);
double a1 = 0.0, a2 = 0.0, a3 = 0.0, a4 = 0.0, a5 = 0.0;
double a6 = 0.0, a7 = 0.0, a8 = 0.0, a9 = 0.0, alt = -1.0;
for (int k=1; k<=n; k++){
double k2 = Math.pow(k,2), k3 = k2*k;
double sk = Math.sin(k), ck = Math.cos(k);
alt = -alt;
a1 += Math.pow(twothirds,k-1);
a2 += Math.pow(k,-0.5);
a3 += 1.0/(k*(k+1.0));
a4 += 1.0/(k3 * sk*sk);
a5 += 1.0/(k3 * ck*ck);
a6 += 1.0/k;
a7 += 1.0/k2;
a8 += alt/k;
a9 += alt/(2.0*k -1.0);
}
System.out.printf("%.9f\t(2/3)^k\n", a1);
System.out.printf("%.9f\tk^-0.5\n", a2);
System.out.printf("%.9f\t1/k(k+1)\n", a3);
System.out.printf("%.9f\tFlint Hills\n", a4);
System.out.printf("%.9f\tCookson Hills\n", a5);
System.out.printf("%.9f\tHarmonic\n", a6);
System.out.printf("%.9f\tRiemann Zeta\n", a7);
System.out.printf("%.9f\tAlternating Harmonic\n", a8);
System.out.printf("%.9f\tGregory\n", a9);
System.out.println((System.currentTimeMillis () - start) + "ms");
}
}