map <char,CMatrix> map2;
int  type (string r, char A[],int y )
{    int t=0,h=0,d;
     int x = r.length()-1;
	for (int i=0;i<y;i++)
    {
	if (r[0]==A[i])
         t=1;
	}
	for (int i=0;i<y;i++)
   {
		if (r[x]==A[i])
			 h=1;}
	if(t==1 && h==1)
		d=1;
	else if (t==1 || h==1)
	d=2;
	else
		d=3;
	return d;}


void get_matrix(string &s1,char A[],int y,char CHAR)
{   string r;

	//int flag=0,FirstOperation=0;               //to count the number of the operations between two brackets ()
	char elements [100];                    //contains the operands of the input string * , - , + , / , ^
	int places[100];                        //contains the index of each operand
    string values[100];                    //contains all the numbers in the input string in order in string form
	//double Values[100]={};					   //contains all the numbers in the input string in order in double form
	//double Result;                             //contains the final result of the input string as a double form
	string result;                             //contains the final result of the input string as a string form
	int j=0;
	for (int i = 0; i < s1.length(); i++)
	{
		if (s1[i]=='+'||s1[i]=='-'||s1[i]=='*'||s1[i]=='/'||s1[i]=='^')
		{
			//flag++;
			elements[j]=s1[i];
			places[j]=i;
			if (j==0)
			{
				values[j]=s1.substr(0,i);
				//Values[j]=atof(values[j].c_str());
			}
			else
			{
				values[j]=s1.substr(places[j-1]+1,places[j]-places[j-1]-1);
				//Values[j]=atof(values[j].c_str());
			}
			j++;
		}
	}
	values[j]=s1.substr(places[j-1]+1,s1.size()-places[j-1]-1);
	/*if (flag==0)
	{
		break;
	}*/

		for (int i = 0; i < j; i++)
	{
			if (elements[i]=='^')
            {

		         r=values[i]+string(1, elements[i])+values[i+1];
			     int typee;
                 typee= type ( r,A, y );

			if(typee==2)

		   {
		       if(values[1]=='.')

			{
			 map<char,CMatrix>::iterator name;
             name=map2.[Values[i]];
			if(name!=map2.end())
			{

			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i]);
		    double   b=Values[j+1]=atof(values[j+1].c_str());
            (it1->second).elepow(b);
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
			}
			}
			else
			{

			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i+1]);
		    double b =Values[j]=atof(values[j+1].c_str());
			map2[CHAR]= (it1->second).elepow(b);
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
			}
		   }
		   }

			else if(typee==3)

			{
			Values[i]=atof(values[j].c_str());
			Values[i+1]=atof(values[j].c_str());
			values[i+1]= pow(Values[i],Values[i+1]);

			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			CMatrix M(values[i+1]);
			map2[CHAR] = M;
			break;
			/*if(i!=0)
			{arrayModification(i,j,Values);
			arrayModification(i,j,elements);
			arrayModification(i,j,places);}*/
			}}




			if (elements[i]=='*')

		{
			r=Values[i]+string(1, elements[i])+Values[i+1];
			int typee;
            typee= type ( r,  A, y );
			if(typee==1)
			{
			CMatrix M();
			map<char,CMatrix>::iterator it1;           // iterator that iterates over the map to find the desired value(CMatrix) of a certain key(char)
			map<char,CMatrix>::iterator it2;
			it1 =map2.find(Values[i]);      //find the CMatrix of the first operand
			it2 =map2.find(Values[i+1]);
            map2[CHAR]= (it1->second) * (it2->second) ;
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;

			/*if(i!=0)
			{arrayModification(i,j,Values);
			arrayModification(i,j,elements);
			arrayModification(i,j,places);
            }*/


        // mymap.insert( pair < char , CMatrix > (CHAR,CMatrix(M)));


			}

			else if(typee==2)

		   {  map<char,CMatrix>::iterator name;
              name=map2.[Values[i]];
			if(name !=mymap.end())
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i]);
		    double   b=Values[j+1]=atof(values[j+1].c_str());
			//if(b==1)
			map2[CHAR]= (it1->second).elemul(b);
			//else mymap[CHAR]= (it1->second).elemul(b);
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
           /* if(i!=0)
			{arrayModification(i,j,Values);
			arrayModification(i,j,elements);
			arrayModification(i,j,places);}*/

			}

			else if(typee==3)
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i+1]);
		    double   b =Values[j]=atof(values[j+1].c_str());

			map2[CHAR]= (it1->second).elemul(b);
			//else mymap[CHAR]= (it1->second).elemul(b);

			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
            /*if(i!=0)
			{arrayModification(i,j,Values);
			arrayModification(i,j,elements);
			arrayModification(i,j,places);}*/
			}

		   }

			else if(typee==3)

			{
			Values[i]=atof(values[j].c_str());
			Values[i+1]=atof(values[j].c_str());
			values[i+1]= Values[i]*Values[i+1];
			s1.replace(s1.find(r),r.length(),string(1, CHAR));

			CMatrix M(values[i+1]);
			map2[CHAR] = M;
			break;
           /* if (i!=0)
			arrayModification(i+1,j,Values);
			arrayModification(i+1,j,elements);
			arrayModification(i+1,j,places);}*/
			}}



		if (elements[i]=='/')
		{   r=Values[i]+string(1, elements[i])+Values[i+1];
			int typee;
            typee= type ( r,  A , y );
			if(typee==1)
			{
			CMatrix M();
			map<char,CMatrix>::iterator it1;           // iterator that iterates over the map to find the desired value(CMatrix) of a certain key(char)
			map<char,CMatrix>::iterator it2;
			it1 =map2.find(Values[i]);      //find the CMatrix of the first operand
			it2 =map2.find(Values[i+1]);
            map2[CHAR]= (it1->second) / (it2->second) ;
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
            /*if(i!=0)
			{arrayModification(i,j,Values);
			arrayModification(i,j,elements);
			arrayModification(i,j,places);
            }*/



        // mymap.insert( pair < char , CMatrix > (CHAR,CMatrix(M)));


			}

			else if(typee==2)

		   {  map<char,CMatrix>::iterator name;
              name=map2.[Values[i]];
			if(name !=map2.end())
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i]);
		    double   b=Values[j+1]=atof(values[j+1].c_str());
			if(b==1)
			{map2[CHAR]= (it1->second).elediv();

			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;}
			else
				{map2[CHAR]= (it1->second).elediv(b);
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
			/*arrayModification(i+1,j,Values);
			arrayModification(i+1,j,elements);
			arrayModification(i+1,j,places);*/}
			}

			else
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i+1]);
		    double   b =Values[j]=atof(values[j+1].c_str());
			if(b==1)
			{map2[CHAR]= (it1->second).elediv();
			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;
			/*arrayModification(i+1,j,Values);
			arrayModification(i+1,j,elements);
			arrayModification(i+1,j,places);*/}


			else{
				map2[CHAR]= (it1->second).elediv(b);
				s1.replace(s1.find(r),r.length(),string(1, CHAR));
				break;
			/*arrayModification(i+1,j,Values);
			arrayModification(i+1,j,elements);
			arrayModification(i+1,j,places);*/}

			}


			}



			else if(typee ==3)

			{
			Values[i]=atof(values[j].c_str());
			Values[i+1]=atof(values[j].c_str());
			values[i+1]= Values[i]/Values[i+1];
			s1.replace(s1.find(r),r.length(),string(1, CHAR));

			CMatrix M(values[i+1]);
			map2[CHAR] = M;
			break;

			}}


		if (elements[i]=='+')
		{   r=Values[i]+string(1, elements[i])+Values[i+1];
			int typee;
            typee= type ( r,  A, y );
			if(typee==1)
			{
			CMatrix M();
			map<char,CMatrix>::iterator it1;           // iterator that iterates over the map to find the desired value(CMatrix) of a certain key(char)
			map<char,CMatrix>::iterator it2;
			it1 =map2.find(Values[i]);      //find the CMatrix of the first operand
			it2 =map2.find(Values[i+1]);
            map2[CHAR]= (it1->second) + (it2->second) ;

			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;


        // mymap.insert( pair < char , CMatrix > (CHAR,CMatrix(M)));


			}

			else if(typee==2)

		   {  map<char,CMatrix>::iterator name;
              name=map2.[Values[i]];
			if(name !=map2.end())
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i]);
		    double   b=Values[j+1]=atof(values[j+1].c_str());


			 map2[CHAR]= (it1->second).eleadd(b);
			 s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;

			/*else
			{
			map<char,CMatrix>::iterator it1;
			it1 =mymap.find(Values[i+1]);
		    double   b =Values[j]=atof(values[j+1].c_str());
			mymap[CHAR]= (it1->second).eleadd(b);
			//else mymap[CHAR]= (it1->second).eleadd(b);


			}*/

		   }

			else if(typee==3)

			{
			Values[i]=atof(values[j].c_str());
			Values[i+1]=atof(values[j].c_str());
			values[i+1]= Values[i]+  Values[i+1];
			s1.replace(s1.find(r),r.length(),string(1, CHAR));

			CMatrix M(values[i+1]);
			map2[CHAR] = M;
			break;
			}}


		if (elements[i]=='-')
		{   r=Values[i]+string(1, elements[i])+Values[i+1];
			int typee;
            typee= type ( r,  A, y );
			if(typee==1)
			{
			CMatrix M();
			map<char,CMatrix>::iterator it1;           // iterator that iterates over the map to find the desired value(CMatrix) of a certain key(char)
			map<char,CMatrix>::iterator it2;
			it1 =map2.find(Values[i]);      //find the CMatrix of the first operand
			it2 =map2.find(Values[i+1]);
            map2[CHAR]= (it1->second) - (it2->second) ;


			s1.replace(s1.find(r),r.length(),string(1, CHAR));
			break;



        // mymap.insert( pair < char , CMatrix > (CHAR,CMatrix(M)));


			}

			else if(typee==2)

		   {  map<char,CMatrix>::iterator name;
              name=map2.[Values[i]];
			if(name !=map2.end())
			{
			map<char,CMatrix>::iterator it1;
			it1 =map2.find(Values[i]);
		    double   b=Values[j+1]=atof(values[j+1].c_str());
			//if(b==1)
			//mymap[CHAR]= (it1->second).eleadd();
			 map2[CHAR]= (it1->second).elesub(b);
			  map2[CHAR]= (it1->second).eleadd(b);
			 s1.replace(s1.find(r),r.length(),string(1, CHAR));

			break;

			/*else
			{
			map<char,CMatrix>::iterator it1;
			it1 =mymap.find(Values[i+1]);
		    double   b =Values[j]=atof(values[j+1].c_str());
			//mymap[CHAR]= (it1->second).eleadd();
			 mymap[CHAR]= (it1->second).elesub(b);


			}*/

		   }

			else if(typee==3)

			{
			Values[i]=atof(values[j].c_str());
			Values[i+1]=atof(values[j].c_str());
			values[i+1]= Values[i]-Values[i+1];
			s1.replace(s1.find(r),r.length(),string(1, CHAR));

			CMatrix M(values[i+1]);
			map2[CHAR] = M;
			break;
			}}
}