Matrix manipulation code performing unknown operations -


i having trouble understanding of following code. think middle section of code point out performing multiplication, not sure, can please explain me middle part of code trying accomplish?

i understand input , output parts of array elements, not actual manipulations.

i understand this: 

 #include "stdafx.h"     #include<iostream>     using namespace std;     int main()     {        int i,j,k,n;        float a[100][200],t;        cout<<"enter order of matrix-";        cin>>n;        cout<<"enter elements of matrix"<<endl;        for(i=0;i<n;i++)           for(j=0;j<n;j++)              cin>>a[i][j];

but following lines confusing:

      for(i=0;i<n;i++)        {           for(j=n;j<2*n;j++)           {               if(i==j-n)                  a[i][j]=1;              else                  a[i][j]=0;            }        }        for(i=0;i<n;i++)        {           t=a[i][i];           for(j=i;j<2*n;j++)               a[i][j]=a[i][j]/t;           for(j=0;j<n;j++)           {              if(i!=j)              {                 t=a[j][i];                 for(k=0;k<2*n;k++)                     a[j][k]=a[j][k]-t*a[i][k];               }           }        }

again, following lines understandable:

cout<<"\n\ninverse matrix\n\n";        for(i=0;i<n;i++)        {           for(j=n;j<2*n;j++)              cout<<"\t"<<a[i][j];           cout<<"\n";         }     return 0;     }

your code attempting invert matrix. here code doing.

assuming have nxn matrix a, can create new nx2n matrix

b = ai

where nxn identity matrix. creation of b these lines trying do:

for(i=0;i<n;i++)        {           for(j=n;j<2*n;j++) // makes sure accessing columns in right hand part of matrix           {               if(i==j-n)                  a[i][j]=1;              else                  a[i][j]=0;            }        }

you perform series of row operations on 2n wide matrix identity matrix on left hand side resulting in third matrix

c=id

where again nxn identity matrix. creation of c through row operations part of code trying do: (though won't speak correctness of algorithm being used has been while since have done , typically use svd algorithm inverses)

for(i=0;i<n;i++) //transform b matrix c matrix        {           t=a[i][i];           for(j=i;j<2*n;j++)               a[i][j]=a[i][j]/t;           for(j=0;j<n;j++)           {              if(i!=j)              {                 t=a[j][i];                 for(k=0;k<2*n;k++)                     a[j][k]=a[j][k]-t*a[i][k];               }           }        }

now when complete row operations , identity matrix on left hand side of c matrix, turns out nxn d matrix makes right hand side of c matrix inverse of a. have

ad = i

which why last bit of code grabs right hand side of called c matrix.

that's code trying do. if want understand more theory behind it, take @ "numerical recipes in c" book. it's gold standard kind of thing , freely available online. http://www.nrbook.com/a/bookcpdf.html


-

Comments

Post a Comment

Popular posts from this blog

PySide and Qt Properties: Connecting signals from Python to QML -

i trying run simple pyside appliccation design made in qml. call function in qml python in separated thread. i have following qml view: import qt 4.7 rectangle { visible: true width: 100; height: 100 color: "lightgrey" text { id: datestr objectname: "datestr" x: 10 y: 10 text: "init text" font.pixelsize: 12 function updatedata(text) { datestr.text = qstr(text) console.log("you said: " + text) } } } here application: some imports.... #!/usr/bin/python import sys import time pyside import qtdeclarative, qtcore pyside.qtgui import qapplication pyside.qtdeclarative import qdeclarativeview the qthread worker.... #worker thread class worker(qtcore.qthread): updateprogress = qtcore.signal(int) def __init__(self, child): qtcore.qthread.__init__(self) self.child = child def run(self):
Read more

c# - DevExpress.Wpf.Grid.InfiniteGridSizeException was unhandled -

Image
when using devexpress, see error: devexpress.wpf.grid.infinitegridsizeexception unhandled message="by default, infinite grid height not allowed since grid rows rendered , hence grid work slowly. fix issue, should place grid container give finite height grid, or should manually specify grid's height or maxheight. note can avoid exception setting gridcontrol.allowinfinitegridsize static property true, in case grid run slowly." the problem dxgrid has infinite height. to fix, set height non-infinite. snoop absolutely invaluable this: if "height" xaml element infinite (i.e. 0 or nan ), can set using one of following: option 1: height="{binding path=actualheight, relativesource={relativesource mode=findancestor, ancestortype=uielement}}" option 2: verticalalignment="stretch" option 3: height="auto" hint: use verticalalignment="stretch" if child of grid <rowdefinition height=
Read more

scala - 'wrong top statement declaration' when using slick in IntelliJ -

Image
i'm porting project scala 2.9 2.10, therefore have use slick instead of scalaquery. i'm using slick 2.1.0 since supports ms access. according this tutorial , upgrade guide changed robs object class , added val robs : why error message wrong top statement declaration , how rid of it? edit: i'm new scala... according this question seems can't put val outside of methods or classes, right? code above directly in packge. right approach slick then? should move code in class or trait? according this changed val robs = tablequery[robs] to object robs extends tablequery(new robs(_)) {} . no warnings or errors anymore. :)
Read more