/*
 * Test_Sparse_jac.cpp
 *
 *  Created on: May 17, 2013
 *      Author: ad007804
 */

/*----------------------------------------------------------------------------
 ADOL-C -- Automatic Differentiation by Overloading in C++
 File:     sparse_jacobian.cpp
 Revision: $Id: sparse_jacobian.cpp 299 2012-03-21 16:08:40Z kulshres $
 Contents: example for computation of sparse jacobians

 Copyright (c) Andrea Walther, Andreas Griewank, Andreas Kowarz,
               Hristo Mitev, Sebastian Schlenkrich, Jean Utke, Olaf Vogel

 This file is part of ADOL-C. This software is provided as open source.
 Any use, reproduction, or distribution of the software constitutes
 recipient's acceptance of the terms of the accompanying license file.

---------------------------------------------------------------------------*/

#include <math.h>
#include <cstdlib>
#include <cstdio>

#include <adolc/adolc.h>
#include <adolc/adolc_sparse.h>

#include "sys/types.h"
#include "sys/sysinfo.h"

#define tag 1

void   ceval_ad(adouble *x, adouble *c);
void   ceval(double *x, double *c);

#define m 10000
int main() {

	adouble factive[m];
	double f[m];
	double xpassive[m];

	adouble xactive[m];
	adouble a[m], b[m];
	for (int q=0;q<m;q++)
		xpassive[q] = 1.+q;

	while(1)
	{
		trace_on(0);

		for (int q=0;q<m;q++)
			xactive[q]<<=xpassive[q];
		for (int q=0;q<m;q++)
		{
			a[q] = 2*xactive[q];
			b[q] = 4*xactive[q];
		}
		for (int q=0;q<m;q++)
			factive[q] = a[q]*b[q];

		for (int q=0;q<m;q++)
			factive[q]>>=f[q];
		trace_off();

		double **Jac;
		Jac = myalloc2(m,m);
		jacobian(0,m,m,xpassive,Jac);

		//Check RAM usage
		struct sysinfo memInfo;
		sysinfo (&memInfo);

		long long physMemUsed = memInfo.totalram - memInfo.freeram;
		physMemUsed *= memInfo.mem_unit;

		std::cout<<"physMemUsed: "<< physMemUsed<<std::endl;
	}


    printf("\n");

}


/***************************************************************************/

void ceval(double *x, double *c) {
    c[0] = 2*x[0]+x[1]-2.0;
    c[1] = x[2]*x[2]+x[3]*x[3]-2.0;
    c[2] = 3*x[4]*x[5] - 3.0;
}

/***************************************************************************/

void ceval_ad(adouble *x, adouble *c) {
   // x[0] = x[1];
	c[0] = 2*x[0]+x[1]-2.0;
    c[0] += cos(x[3])*sin(x[4]);
    c[1] = x[2]*x[2]+x[3]*x[3]-2.0;
    c[2] = 3*x[4]*x[5] - 3.0+sin(x[4]*x[5]);
}

/***************************************************************************/