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

    TRAVIS - Trajectory Analyzer and Visualizer

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

    Copyright (c) 2009-2026 Martin Brehm
                  2012-2017 Martin Thomas
                  2016-2020 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 "xdvec3array.h"


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


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


#ifndef DEBUG_ARRAYS
#define m_sName (NULL)
#endif


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


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

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

	
CxDVec3Array::CxDVec3Array(const CxDVec3Array &o) : CxObject()
{
	BXIN;
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::CxDVec3Array(CxDVec3Array &)...");
#endif
//	unsigned long z;
	int z;
	m_iGrow = o.m_iGrow;
//	m_pData = new CxDVector3*[m_iMaxSize];

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

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

		try { m_pData = new CxDVector3[m_iMaxSize]; } catch(...) { m_pData = NULL; }
		if (m_pData == NULL) NewException((double)m_iMaxSize*sizeof(CxDVector3),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
	
		if (o.m_pData != NULL) {
			//memcpy(m_pData,o.m_pData,sizeof(CxDVector3)*m_iSize);
			for (z=0;z<(int)m_iSize;z++)
				m_pData[z] = CxDVector3(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_CDVEC3ARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDVec3Array::SetAt(unsigned long pos, CxDVector3 o)
{
	BXIN;
#ifdef DEBUG_CDVEC3ARRAY
	bool s = false;
	mprintf("@ CxDVec3Array::Add(CxDVector3)...");
#endif
	if (pos >= m_iMaxSize)
	{
#ifdef DEBUG_CDVEC3ARRAY
		s = true;
		mprintf("\n");
#endif
		SetMaxSize(pos+1);
	}
//	*m_pData[pos] = *o;
	m_pData[pos] = o;
	if (pos >= m_iSize)
		m_iSize = pos+1;
#ifdef DEBUG_CDVEC3ARRAY
	if (s)
		mprintf("@ done.\n");
	else mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDVec3Array::Add(CxDVector3 o)
{
	BXIN;
#ifdef DEBUG_CDVEC3ARRAY
	bool s = false;
	mprintf("@ CxDVec3Array::Add(CxDVector3)...");
#endif
	if (m_iSize+1 > m_iMaxSize)
	{
#ifdef DEBUG_CDVEC3ARRAY
		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] = *o;
	m_pData[m_iSize] = o;
	m_iSize++;
#ifdef DEBUG_CDVEC3ARRAY
	if (s)
		mprintf("@ done.\n");
			else mprintf("done.\n");
#endif
	BXOUT;
}

	
void CxDVec3Array::SetSize(unsigned long i)
{
	BXIN;
	int z;
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::SetSize(int): %d...",i);
#endif
	if (m_iSize == i)
		return;
//	CxDVector3 **temp;
	CxDVector3 *temp;
//	temp = new CxDVector3*[i];

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


void CxDVec3Array::SetMaxSize(unsigned long i)
{
	BXIN;
	int z;
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::SetMaxSize(int): %d...",i);
#endif
//	CxDVector3 **temp;
	CxDVector3 *temp;
	if (i <= m_iSize)
	{
		eprintf("CxDVec3Array::SetMaxSize(%lu): Size %lu > %lu.\n",i,m_iSize,i);
		return;
	}
//	temp = new CxDVector3*[i];

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

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


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

	
void CxDVec3Array::RemoveAt(unsigned long pos, unsigned long count)
{
	BXIN;
//	CxDVector3 **temp;
	CxDVector3 *temp;
	int z;
		
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::RemoveAt(int, int): %d, %d...",pos,count);
#endif
//	temp = new CxDVector3*[m_iSize-count];

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

	
void CxDVec3Array::InsertAt(CxDVector3 o, unsigned long pos)
{
	BXIN;
	int z;
//	CxDVector3 **temp;
	CxDVector3 *temp;
		
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::InsertAt(CxDVector3, int): %d...");
#endif
//	temp = new CxDVector3*[m_iSize+1];

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


void CxDVec3Array::CopyFrom(CxDVec3Array *a)
{
	BXIN;
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::CopyFrom(CxDVec3Array*)...");
#endif
	unsigned long z;
	if (m_iMaxSize != a->m_iMaxSize)
	{
		if (m_pData != NULL)
		{
/*			for (z=0;z<m_iSize;z++)
				delete m_pData[z];*/
			delete[] m_pData;
		}
		m_iMaxSize = a->m_iMaxSize;
//		m_pData = new CxDVector3*[m_iMaxSize];

		try { m_pData = new CxDVector3[m_iMaxSize]; } catch(...) { m_pData = NULL; }
		if (m_pData == NULL) NewException((double)m_iMaxSize*sizeof(CxDVector3),__FILE__,__LINE__,__PRETTY_FUNCTION__,m_sName);
		
/*		for (z=0;z<m_iMaxSize;z++)
			m_pData[z] = new CxDVector3();*/
	}
	m_iSize = a->m_iSize;
	m_iGrow = a->m_iGrow;
	for (z=0;z<m_iSize;z++)
		m_pData[z] = a->m_pData[z];
//		*m_pData[z] = *a->m_pData[z];
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("done.\n");
#endif
	BXOUT;
}


void CxDVec3Array::SetName(const char *name)
{
	UNUSED(name);
#ifdef DEBUG_ARRAYS
	BXIN;
#ifdef DEBUG_CDVEC3ARRAY
	mprintf("@ CxDVec3Array::SetName(const char *s): \"%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_CDVEC3ARRAY
	mprintf("done.\n");
#endif
	BXOUT;
#endif
}
