#include <vector>
#include <iostream>
#include <string>
using namespace std;

int maxN = 30000;

double TestSTL(void);
double TestArray(void);


int main(int narg, char *argv[])
{

	string what(narg>1 ? argv[1]:"");

	double sumXsq = 0.0;
	if(what=="array"){
		sumXsq = TestArray();
	}else if(what=="stl"){
		sumXsq = TestSTL();
	}else{
		cout<<"Must pass either \"array\" or \"stl\" as argument!"<<endl;
	}

	cout<<"sumXsq = "<<sumXsq<<endl;
}

double TestSTL(void)
{
  vector<double> xxx(maxN);

  for( int i = 0; i < maxN; i++ )xxx[i] = i;

  double sumXsq =0;
  for( int i = 0; i < maxN; i++)
    for( int j =0; j < maxN; j++)
      sumXsq += xxx[i] * xxx[i];

  return sumXsq;
}

double TestArray(void)
{
  vector<double> xxx;
  xxx.resize(maxN);
  double *zzz = new double  [maxN];

  for( int i = 0; i < maxN; i++ )zzz[i] = i;

  double sumXsq =0;
  for( int i = 0; i < maxN; i++)
    for( int j =0; j < maxN; j++)
      sumXsq += zzz[i] * zzz[i];

  return sumXsq;
}