From d87857fa398087026805cbffa90c6f15d2c6c1c1 Mon Sep 17 00:00:00 2001
From: Mubashshir <ahmubashshir@users.noreply.github.com>
Date: Mon, 26 Mar 2018 01:17:51 +0600
Subject: [PATCH] Update

* Added Auto indentation support.
* Cross Platform support.
---
 Makefile         |   7 +
 examples/mat.cpp |  78 ++++++++
 include/mat.h    | 476 +++++++++++++++++++++++++++++++++++++++++++++++
 3 files changed, 561 insertions(+)
 create mode 100644 Makefile
 create mode 100644 examples/mat.cpp
 create mode 100644 include/mat.h

diff --git a/Makefile b/Makefile
new file mode 100644
index 0000000..ae49152
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,7 @@
+CPP=g++
+CPPFLAGS=-Wall -fexceptions -g -fPIC
+.PHONY: all
+all: mat
+mat:
+	@echo "Building example"
+	$(CPP) -Wall -I include examples/mat.cpp -o mat
diff --git a/examples/mat.cpp b/examples/mat.cpp
new file mode 100644
index 0000000..95f427c
--- /dev/null
+++ b/examples/mat.cpp
@@ -0,0 +1,78 @@
+
+/*
+ * 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 2 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, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
+ * MA 02110-1301, USA.
+ * 
+ */
+#include<iostream>
+#define __debug__
+#define _limit_ 5
+#include <mat.h>
+using namespace std;
+int main()
+{
+    int line=7;
+	//~ try
+	{
+	   	unsigned int col,row;
+	   	line+=3;
+        //cout<<"\n1";
+        cout<<"Row:";
+	   	line+=3;
+        //cout<<"\n2";
+        cin>>row;
+	   	line+=3;
+        //	cout<<"\n3";
+        cout<<"Column:";
+	   	line+=3;
+        //cout<<"\n4";
+        cin>>col;
+	   	line+=3;
+        //	cout<<"\n5";
+        line+=4;
+        mat <int> mata(row,col,"MatA");
+        //cout<<"\n6";
+        mat <int> matb(row,col,"MatB");
+	   	line+=3;
+        //	cout<<"\n7";
+        mat <int> ans(row,col,"MatAns");
+	   	line+=3;
+        //cout<<"\n8";
+        mata.in();
+	   	line+=3;
+        //cout<<"\n9";
+        matb.in();
+	   	line+=3;
+        //cout<<"\n10";
+        mata.out();
+	   	line+=3;
+        //cout<<"\n11";
+        matb.out();
+	   	line+=3;
+        //cout<<"\n12";
+        ans.add(mata,matb);
+	   	line+=3;
+        //cout<<"\n13";
+        ans.out();
+	   	line+=3;
+        //cout<<"\n14";
+	}
+	//catch(...)
+	//{
+		//cout<<"Exception Occurred at Line "<<line<<endl;
+		//return -3;
+	//}
+	return 0;
+}
diff --git a/include/mat.h b/include/mat.h
new file mode 100644
index 0000000..a77cf5c
--- /dev/null
+++ b/include/mat.h
@@ -0,0 +1,476 @@
+
+/*
+ * matrics.hpp
+ * 
+ * Copyright 2018 Ahmad Hasan Mubashshir <ahmubashshir@gmail.com>
+ * 
+ * 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 2 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, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
+ * MA 02110-1301, USA.
+ * 
+ * 
+ */
+
+#ifndef __MATRICS__
+#define __MATRICS__
+#define __MATRICS_MAGOR__ 3
+#define __MATRICS_MINOR__ 2
+#ifndef devN
+#define devN "\x41\x68\x6D\x61\x64\x20\x48\x61\x73\x61\x6E\x20\x4D\x75\x62\x61\x73\x68\x73\x68\x69\x72"
+#endif // devN
+#ifndef devC
+#define devC "\x50\x6F\x6C\x79\x6D\x61\x74\x68\x20\x54\x65\x61\x6d"
+#endif // devC
+#if defined(linux) || defined(__linux__) || defined(__linux) || defined(unix) || defined(__unix) || defined(__unix__) || defined(__APPLE__) || defined(__MACH__) || defined(__FreeBSD__) || defined(__hpux) || defined(_AIX)
+	#define STARTL "\342\216\241"
+	#define MIDL "\342\216\242"
+	#define ENDL "\342\216\243"
+	#define STARTR "\342\216\244"
+	#define MIDR "\342\216\245"
+	#define ENDR "\342\216\246"
+#elif defined(_WIN32)||defined(__CYGWIN__)||defined(__MINGW32__)||defined(__MINGW64__)||defined(__WINNT__)
+	#define STARTL "\xDA"
+	#define MIDL "\xB3"
+	#define ENDL "\xC0"
+	#define STARTR "\xBF"
+	#define MIDR "\xB3"
+	#define ENDR "\xD9"
+#endif //OS Specific options
+#if defined(__debug__)
+    #define Destruct1 std::cout<<this->name<<"["<<i<<"] Destructing"<<std::endl
+    #define Destruct2 std::cout<<"Destructing "<<this->name<<std::endl;
+    #define Construct1 std::cout<<this->name<<"["<<i<<"] Constructing"<<std::endl
+    #define Construct2 std::cout<<"Constructing "<<this->name<<std::endl;
+#else
+    #define Destruct1 std::cout<<""
+    #define Destruct2 std::cout<<""
+    #define Construct1 std::cout<<""
+    #define Construct2 std::cout<<""
+#endif // Debug Options
+#ifdef _limit_
+    #define exception_handle_1 {std::cout<<"Buffer Overflow Error\n";throw col;}
+    #define exception_handle_2 {std::cout<<"Buffer Overflow Error\n";throw row;}
+    #if _limit_ > 2
+		#define exception_limit row<_limit_&&col<_limit_
+    #else
+		#define exception_limit row<5&&col<5
+	#endif
+    #define exception_handle_buuf {std::cout<<"Buffer Underflow Error\n";throw row+col;}
+#else
+	#define exception_limit true
+	#define exception_handle_buuf std::cout<<""
+    #define exception_handle_1 std::cout<<""
+    #define exception_handle_2 std::cout<<""
+    #warning "You should use limited element matrics to protect from buffer overflow by adding the next line"
+    #warning "#define _limit_ <max-index>+1"
+#endif // Buffer Limit
+#if defined(_noauto)
+#define indent if((this->ret[i][j]/100000)>0) std::cout<<" ";\
+else if((this->ret[i][j]/10000)>0) std::cout<<"  ";\
+else if((this->ret[i][j]/1000)>0) std::cout<<"   ";\
+else if((this->ret[i][j]/100)>0) std::cout<<"    ";\
+else if((this->ret[i][j]/10)>0) std::cout<<"     ";\
+else std::cout<<"      ";
+#else
+#define indent if(get_len(j)==int_len((long int)ret[i][j])) std::cout<<" ";\
+else{unsigned lenmax=get_len(j),len2=int_len((long int)ret[i][j]);for(unsigned z=0;z<lenmax-len2+1;z++) {std::cout<<" ";}}
+#endif //Auto Indent
+#include<cstdlib>
+#include<iostream>
+#include<cstring>
+#include<cstdio>
+#include<cstdlib>
+#include<iostream>
+#include<cstring>
+#include<stdexcept>
+template <class Matrics>
+class mat
+{
+    private:
+        unsigned int row,col;
+		Matrics** ret;
+		const char *name;bool dimerr;char *op1,*op2,*op;
+	public:
+		/*
+		* name: mat::mat()
+		* dsc: Constructor
+		* @param Row,Col,Name
+		* @return None
+		* @auto-call:yes
+		* @overloaded:no
+		*/		
+		mat(unsigned int row,unsigned int col,const char *name)
+		{
+			
+			this->dimerr=false;
+			this->row=row;this->col=col;this->name=name;
+			if(col<=0&&row<=0) exception_handle_buuf;
+			if(row>0&&exception_limit) this->ret=new Matrics* [row];
+			else exception_handle_2;
+			for(unsigned  int i=0;i<row;i++)
+			{
+				if(col>0)this->ret[i]=new Matrics[col];
+				Construct1;
+			}
+			Construct2;
+		}
+		/*
+		* name: mat::~mat()
+		* dsc: Destructor
+		* @param None
+		* @return None
+		* @auto-call:yes
+		* @overloaded:no
+		*/
+		~mat()
+		{
+			for(unsigned int i=0;i<this->row;i++)
+			{
+				delete [] this->ret[i];
+				Destruct1;
+        	}
+			Destruct2;
+		}
+		/*
+		* name: mat::in()
+		* dsc: Input function
+		* @param None
+		* @return None
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		void in(void)
+		{
+	        for(unsigned int i=0;i<this->row;i++)
+			{
+				std::cout<<"Row "<<i<<":";
+        	    for(unsigned int j=0;j<this->col;j++)
+        	    {
+                        	std::cin>>this->ret[i][j];
+				}
+        	}
+		}
+		/*
+		* name: mat::int_len()
+		* dsc: Int length calculator
+		* @param Int A
+		* @return unsigned length
+		* @auto-call: no
+		* @overloaded:no
+		*/
+		unsigned int_len(long int a)
+        {
+            long int n=abs(a);
+            if(n<10) return 1;
+            else
+            {
+                unsigned r=1;
+                while(n>=10)
+                {
+                    n/=10;
+                    r++;
+                }
+                return r;
+            }
+        }
+        /*
+		* name: mat::out()
+		* dsc: Output Function
+		* @param None
+		* @return None
+		* @auto-call: no
+		* @overloaded:no
+		*/
+		void out();
+		/*
+		* name: mat::add()
+		* dsc: Matrics Addition
+		* @param Mat A,Mat B
+		* @return None
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		void add(mat matA,mat matB)
+		{
+			if((!matA.dimerr)&&(!matB.dimerr))
+			{
+				if((matA.col==matB.col)&&(matA.row==matB.row)&&(this->col==matA.col)&&(this->row==matA.row))
+				{
+					for(unsigned int i=0;i<this->row;i++)
+					{
+						for(unsigned int j=0;j<this->col;j++)
+						{
+							this->ret[i][j]=matA.ret[i][j]+matB.ret[i][j];
+						}
+					}
+				}
+				else
+					this->dimerr=true;
+			}
+			else
+				this->dimerr=true;
+		}
+		/*
+		* name: mat::sub()
+		* dsc: Matrics Subtraction
+		* @param Mat A,Mat B
+		* @return None
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		void sub(mat matA,mat matB)
+		{
+			if((!matA.dimerr)&&(!matB.dimerr))
+			{
+				if((matA.col==matB.col)&&(matA.row==matB.row)&&(this->col==matA.col)&&(this->row==matA.row))
+				{
+					for(unsigned int i=0;i<this->row;i++)
+					{
+						for(unsigned int j=0;j<this->col;j++)
+						{
+							this->ret[i][j]=matA.ret[i][j]-matB.ret[i][j];
+						}
+					}
+				}
+				else
+					this->dimerr=true;
+			}
+			else
+				this->dimerr=true;
+		}
+		/*
+		* name: mat::mul()
+		* dsc: Matrics Multiplication(Scaler)
+		* @param Int Multiple
+		* @return None
+		* @auto-call:no
+		* @overloaded:yes
+		*/
+		void mul(int multiple)
+		{
+			if(!this->dimerr)
+			{
+				for(unsigned int i=0;i<this->row;i++)
+					for(unsigned int j=0;j<this->col;j++)
+						this->ret[i][j]=multiple*this->ret[i][j];
+			}
+		}
+		/*
+		* name:  mat::mul()
+		* dsc: Matrics Multiplication(Matrics)
+		* @param Mat A,Mat B
+		* @return None
+		* @auto-call:no
+		* @overloaded:yes
+		*/
+		void mul(mat matA,mat matB)
+		{
+			if((!matA.dimerr)&&(!matB.dimerr))
+			{
+				if((matA.col==matB.row)&&(this->row==matA.row)&&(this->col==matB.col))
+				{
+					for(unsigned int i=0;i<this->row;i++)
+						for(unsigned int j=0;j<this->col;j++)
+							for(unsigned int k=0;k<this->row;k++)
+								this->ret[i][j]+=matA.ret[i][k]*matB.ret[k][j];
+				}
+				else this->dimerr=true;
+			}
+			else this->dimerr=true;
+		}
+		/*
+		* name: mat::reinit()
+		* dsc: Reinitializer
+		* @param Row,Col,Name
+		* @return None
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		void reinit(unsigned int row,unsigned col,const char* name)
+		{
+			//Clear::TODO:exception handler
+			for(unsigned int i=0;i<this->row;i++)
+			{
+				delete [] this->ret[i];
+				Destruct1;
+        	}
+			Destruct2;
+			//Init::TODO:exception handler
+			this->dimerr=false;
+			this->row=row;this->col=col;this->name=name;
+			if(row>0&&exception_limit) this->ret=new Matrics* [row];
+			else exception_handle_2;
+			for(unsigned  int i=0;i<row;i++)
+			{
+				if(col>0) this->ret[i]=new Matrics[col];
+				else exception_handle_1;
+				Construct1;
+			}
+			Construct2;
+		}
+		/*
+		* name: mat::get_col()
+		* dsc: Column
+		* @param None
+		* @return Column
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		unsigned get_col(void)
+		{
+		    return this->col;
+		}
+		/*
+		* name: mat::row()
+		* dsc: Row
+		* @param none
+		* @return Row
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		int get_row(void)
+		{
+		    return this->row;
+		}
+		/*
+		* name: mat::get_raw()
+		* dsc: Raw Value
+		* @param None
+		* @return Matrics Array
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		Matrics** get_raw(void)
+		{
+		    return this->ret;
+		}
+		/*
+		* name: mat::get_len()
+		* dsc: Get length
+		* @param C
+		* @return Length
+		* @auto-call:no
+		* @overloaded:no
+		*/
+		unsigned get_len(unsigned c)
+        {
+            unsigned maxlen=0,len=0;
+            for(unsigned int i=0;i<row;i++)
+            {
+                len=int_len((long int)this->ret[i][c]);
+                maxlen=(len>=maxlen)?len:maxlen;
+            }
+            return maxlen;
+        }
+};
+template<class Matrics>
+void mat<Matrics>::out(void)
+{
+	if(!this->dimerr)
+	{
+		//if(strlen(op1)>0) std::cout<<op;else
+		std::cout<<this->name;
+		std::cout<<"\n";
+		for(unsigned int i=0;i<this->row;i++)
+		{
+			for(unsigned int j=0;j<this->col;j++)
+			{
+				if((this->row==1)&&(j==0)) std::cout<<"[";
+				if((this->row>=2)&&(j==0)&&(i==0)) std::cout<<STARTL;
+				if((this->row>1)&&(j==0)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDL;
+				if((this->row>=2)&&(j==0)&&(this->row-i==1)) std::cout<<ENDL;
+				if((j==0)&&(this->row>1)) std::cout<<" ";
+				std::cout<<this->ret[i][j];
+                indent
+				if((this->row==1)&&(j>=0)&&((this->col-j)>1)) std::cout<<" ";
+				if((this->row>=2)&&(j==this->col-1)&&(i==0)) std::cout<<STARTR;
+				if((this->row>1)&&(this->col-j==1)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDR;
+				if((this->row>=2)&&(j==this->col-1)&&(this->row-i==1)) std::cout<<ENDR;
+				if((j==this->col-1)&&(this->row>1)) std::cout<<"\n";
+				if((this->row==1)&&(j==this->col-1)) std::cout<<"]";
+			}
+		}
+	} else std::cout<<"Dimension error\n";
+}
+template<>
+void mat<double>::out()
+{
+	if(!this->dimerr)
+	{
+		//if(strlen(op1)>0) std::cout<<op;else
+		std::cout<<this->name;
+		std::cout<<"\n";
+		for(unsigned int i=0;i<this->row;i++)
+		{
+			for(unsigned int j=0;j<this->col;j++)
+			{
+				if((this->row==1)&&(j==0)) std::cout<<"[";
+				if((this->row>=2)&&(j==0)&&(i==0)) std::cout<<STARTL;
+				if((this->row>1)&&(j==0)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDL;
+				if((this->row>=2)&&(j==0)&&(this->row-i==1)) std::cout<<ENDL;
+				if((j==0)&&(this->row>1)) std::cout<<" ";
+				std::cout<<this->ret[i][j];
+                if(get_len(j)==int_len((long int)ret[i][j])) std::cout<<" ";
+                else
+                {
+                    unsigned lenmax=get_len(j),len2=int_len((long int)ret[i][j]);
+                    for(unsigned z=0;z<lenmax-len2+1;z++)
+                    std::cout<<" ";
+                }
+				if((this->row==1)&&(j>=0)&&((this->col-j)>1)) std::cout<<" ";
+				if((this->row>=2)&&(j==this->col-1)&&(i==0)) std::cout<<STARTR;
+				if((this->row>1)&&(this->col-j==1)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDR;
+				if((this->row>=2)&&(j==this->col-1)&&(this->row-i==1)) std::cout<<ENDR;
+				if((j==this->col-1)&&(this->row>1)) std::cout<<"\n";
+				if((this->row==1)&&(j==this->col-1)) std::cout<<"]";
+			}
+		}
+	} else std::cout<<"Dimension error\n";
+}
+template<>
+void mat<float>::out(void)
+{
+	if(!this->dimerr)
+	{
+		//if(strlen(op1)>0) std::cout<<op;else
+		std::cout<<this->name;
+		std::cout<<"\n";
+		for(unsigned int i=0;i<this->row;i++)
+		{
+			for(unsigned int j=0;j<this->col;j++)
+			{
+				if((this->row==1)&&(j==0)) std::cout<<"[";
+				if((this->row>=2)&&(j==0)&&(i==0)) std::cout<<STARTL;
+				if((this->row>1)&&(j==0)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDL;
+				if((this->row>=2)&&(j==0)&&(this->row-i==1)) std::cout<<ENDL;
+				if((j==0)&&(this->row>1)) std::cout<<" ";
+				else printf("%.2f",this->ret[i][j]);
+                if(get_len(j)==int_len((long int)ret[i][j])) std::cout<<" ";
+                else
+                {
+                    unsigned lenmax=get_len(j),len2=int_len((long int)ret[i][j]);
+                    for(unsigned z=0;z<lenmax-len2+1;z++)
+                    std::cout<<" ";
+                }
+				if((this->row==1)&&(j>=0)&&((this->col-j)>1)) std::cout<<" ";
+				if((this->row>=2)&&(j==this->col-1)&&(i==0)) std::cout<<STARTR;
+				if((this->row>1)&&(this->col-j==1)&&(i!=0)&&(this->row-i!=1)) std::cout<<MIDR;
+				if((this->row>=2)&&(j==this->col-1)&&(this->row-i==1)) std::cout<<ENDR;
+				if((j==this->col-1)&&(this->row>1)) std::cout<<"\n";
+				if((this->row==1)&&(j==this->col-1)) std::cout<<"]";
+			}
+		}
+	} else std::cout<<"Dimension error\n";
+}
+#endif // __MATRICS__