Tuesday, 5 September 2017

C-Program to Check whether given Matrix is Orthogonal or not



#include<stdio.h>

int main()
{
 int i,j,k,size,flag=1;
 printf("Enter size of square matrix i.e 3 for 3x3 and 4 for 4x4 etc\n");
 scanf("%d",&size);
 int a[size][size],temp[size][size],identityCheck[size][size],sum=0;
 printf("Enter Elements in matrix\n");
 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
   scanf("%d",&a[i][j]);
  }
 }
 //Transpose of given matrix
 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
   temp[i][j]=a[j][i];
  }
 }

 //Multiplying given A with A^T
 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
   for(k=0;k<size;k++)
   {
    sum+=(a[i][k]*temp[k][j]);
   }
   identityCheck[i][j]=sum;
   sum=0;
  }
 }

 //Checking whether A*A^T resulted in Identity Matrix or not

 for(i=0;i<size;i++)
 {
  for(j=0;j<size;j++)
  {
   if(i==j && identityCheck[i][j]!=1)
   flag=0;
   if(i!=j && identityCheck[i][j]!=0)
   flag=0;
  }
 }

 //Checking for Orthogonal Matrix or not
 if(flag==1)
 {
  printf("Given Matrix is Orthogonal Matrix\n");
 }
 else
 {
  printf("Given Matrix is not an Orthogonal Matrix\n");
 }
return 0;
}

No comments:

Post a Comment