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

    TRAVIS - Trajectory Analyzer and Visualizer

    http://www.travis-analyzer.de/

    Copyright (c) 2009-2022 Martin Brehm
                  2012-2022 Martin Thomas
                  2016-2022 Sascha Gehrke

    Please cite:  J. Chem. Phys. 2020, 152 (16), 164105.         (DOI 10.1063/5.0005078 )
                  J. Chem. Inf. Model. 2011, 51 (8), 2007-2023.  (DOI 10.1021/ci200217w )

    This file was written by Martin Brehm.

    ---------------------------------------------------------------------------

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.

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


// This must always be the first include directive
#include "config.h"

#include "xdoublearray.h"


const char *GetRevisionInfo_xdoublearray(unsigned int len) {
	static char buf[256];
	GET_REVISION_INFO( buf, len );
	return buf;
}


const char *GetSourceVersion_xdoublearray() {
	static char buf[256];
	GET_SOURCE_VERSION( buf );
	return buf;
}


#ifndef DEBUG_ARRAYS
#define m_sName (NULL)
#endif


CxDoubleArray::CxDoubleArray()
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::CxDoubleArray()\n");
#endif
	m_pData = NULL;
#ifdef DEBUG_ARRAYS
	m_sName = NULL;
#endif
	m_iSize = 0;
	m_iMaxSize = 0;
	m_iGrow = 16;
	BXOUT;
}


CxDoubleArray::CxDoubleArray(const char *name)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::CxDoubleArray(const char *)\n");
#endif
	m_pData = NULL;
#ifdef DEBUG_ARRAYS
	m_sName = NULL;
#endif
	m_iSize = 0;
	m_iMaxSize = 0;
	m_iGrow = 10;
	SetName(name);
	BXOUT;
}

	
CxDoubleArray::~CxDoubleArray()
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::~CxDoubleArray()\n");
#endif
	RemoveAll();
#ifdef DEBUG_ARRAYS
	if (m_sName != NULL)
	{
		delete[] m_sName;
		m_sName = NULL;
	}
#endif
	BXOUT;
}

	
CxDoubleArray::CxDoubleArray(const CxDoubleArray &o) : CxObject()
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::CxDoubleArray(CxDoubleArray &)...");
#endif
	unsigned long z;
	m_iGrow = o.m_iGrow;

	if ((o.m_iMaxSize != 0) && (o.m_pData != NULL)) {

		m_iSize = o.m_iSize;
		m_iMaxSize = o.m_iMaxSize;

		try { m_pData = new double[m_iMaxSize]; } catch(...) { m_pData = NULL; }
		if (m_pData == NULL) NewException((double)m_iMaxSize*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
		
		for (z=0;z<m_iSize;z++)
			m_pData[z] = o.m_pData[z];

	} else {
		m_iSize = 0;
		m_iMaxSize = 0;
		m_pData = NULL;
	}

#ifdef DEBUG_ARRAYS
	if (o.m_sName != NULL)
		SetName(o.m_sName);
	else
		m_sName = NULL;
#endif

#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDoubleArray::SetAt(unsigned long pos, double f)
{
	BXIN;
	unsigned long z;
#ifdef DEBUG_CDOUBLEARRAY
	bool s = false;
	mprintf("@ CxDoubleArray::Add(double)...");
#endif
	if (pos >= m_iMaxSize)
	{
#ifdef DEBUG_CDOUBLEARRAY
		s = true;
		mprintf("\n");
#endif
		SetMaxSize(pos+1);
	}
	m_pData[pos] = f;
	if (pos >= m_iSize)
	{
		for (z=m_iSize;z<pos;z++)
			m_pData[z] = 0.0;
		m_iSize = pos+1;
	}
#ifdef DEBUG_CDOUBLEARRAY
	if (s)
		mprintf("@ done.\n");
	else mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDoubleArray::Add(double f)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	bool s = false;
	mprintf("@ CxDoubleArray::Add(double)...");
#endif
	if (m_iSize+1 > m_iMaxSize)
	{
#ifdef DEBUG_CDOUBLEARRAY
		s = true;
		mprintf("\n");
#endif
//		SetMaxSize(m_iMaxSize + m_iGrow);
		if (m_iMaxSize == 0)
			SetMaxSize(m_iMaxSize + m_iGrow);
				else SetMaxSize(m_iMaxSize*2);
	}
	m_pData[m_iSize] = f;
	m_iSize++;
#ifdef DEBUG_CDOUBLEARRAY
	if (s)
		mprintf("@ done.\n");
			else mprintf("done.\n");
#endif
	BXOUT;
}
	

void CxDoubleArray::SetSize(unsigned long i)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::SetSize(int): %d...",i);
#endif
	if (m_iSize == i)
		return;
	double *temp;

	try { temp = new double[i]; } catch(...) { temp = NULL; }
	if (temp == NULL) NewException((double)i*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	
	if (m_pData != NULL) {
		if (m_iSize > i)
			memcpy(temp,m_pData,i*sizeof(double));
		else
			memcpy(temp,m_pData,m_iSize*sizeof(double));
		delete[] m_pData;
	}
	m_pData = temp;
	m_iMaxSize = i;
	m_iSize = i;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}


void CxDoubleArray::SetMaxSize(unsigned long i)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::SetMaxSize(int): %d...",i);
#endif
	double *temp;

	try { temp = new double[i]; } catch(...) { temp = NULL; }
	if (temp == NULL) NewException((double)i*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	
	if (m_pData != NULL) {
		memcpy(temp,m_pData,m_iSize*sizeof(double));
		delete[] m_pData;
	}
	m_pData = temp;
	m_iMaxSize = i;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDoubleArray::SetGrow(unsigned long i)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::SetGrow(int): %d\n",i);
#endif
	m_iGrow = i;
	BXOUT;
}
		
	
void CxDoubleArray::RemoveAll()
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::RemoveAll():...");
#endif
	if (m_pData != NULL)
		delete[] m_pData;
	m_pData = NULL;
	m_iSize = 0;
	m_iMaxSize = 0;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}


void CxDoubleArray::RemoveAll_KeepSize()
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::RemoveAll_KeepSize():...");
#endif
	m_iSize = 0;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDoubleArray::RemoveAt(unsigned long pos, unsigned long count)
{
	BXIN;
	double *temp;
		
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::RemoveAt(int, int): %d, %d...",pos,count);
#endif

	try { temp = new double[m_iSize-count]; } catch(...) { temp = NULL; }
	if (temp == NULL) NewException((double)(m_iSize-count)*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	
	if (m_pData != NULL) {
		memcpy(temp,m_pData,pos*sizeof(double));
		memcpy(&temp[pos],&m_pData[pos+count],(m_iSize-pos-count)*sizeof(double));
		delete[] m_pData;
	}
	m_pData = temp;
	m_iSize-=count;
	m_iMaxSize = m_iSize;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDoubleArray::InsertAt(double f, unsigned long pos)
{
	BXIN;
	double *temp;
		
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::InsertAt(double, int): %d...");
#endif

	try { temp = new double[m_iSize+1]; } catch(...) { temp = NULL; }
	if (temp == NULL) NewException((double)(m_iSize+1)*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	
	if (m_pData != NULL) {
		memcpy(temp,m_pData,pos*sizeof(double));
		memcpy(&temp[pos+1],&m_pData[pos],(m_iSize-pos)*sizeof(double));
		delete[] m_pData;
	}
	temp[pos] = f;
	m_pData = temp;
	m_iSize++;
	m_iMaxSize = m_iSize;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}


void CxDoubleArray::CopyFrom(CxDoubleArray *d)
{
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::CopyFrom(CxDoubleArray *d): %X...",d);
#endif
	if (m_iMaxSize != d->m_iMaxSize)
	{
		m_iMaxSize = d->m_iMaxSize;
		if (m_pData != NULL)
			delete[] m_pData;

		try { m_pData = new double[d->m_iMaxSize]; } catch(...) { m_pData = NULL; }
		if (m_pData == NULL) NewException((double)d->m_iMaxSize*sizeof(double),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	}
	m_iSize = d->m_iSize;
	m_iGrow = d->m_iGrow;
	if (d->m_pData != NULL)
		memcpy(m_pData,d->m_pData,m_iSize*sizeof(double));
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}


void CxDoubleArray::SetName(const char *name)
{
	UNUSED(name);
#ifdef DEBUG_ARRAYS
	BXIN;
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("@ CxDoubleArray::SetName(const char *): \"%s\"...",name);
#endif
	if (m_sName != NULL)
		delete[] m_sName;
	try { m_sName = new char[strlen(name)+1]; } catch(...) { m_sName = NULL; }
	if (m_sName == NULL) NewException((double)(strlen(name)+1)*sizeof(char),__FILE__,__LINE__,__PRETTY_FUNCTION__);
	strcpy(m_sName,name);
#ifdef DEBUG_CDOUBLEARRAY
	mprintf("done.\n");
#endif
	BXOUT;
#endif
}
