# simpleMatrix code

 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 21, 2013, 05:23 simpleMatrix code #1 Senior Member   Dongyue Li Join Date: Jun 2012 Location: Torino, Italy Posts: 675 Rep Power: 8 Hi, In OpenFOAM's test---simpleMatrix.I have this code: Code: ```int main(int argc, char *argv[]) { simpleMatrix hmm(3); hmm[0][0] = -3.0; hmm[0][1] = 10.0; hmm[0][2] = -4.0; hmm[1][0] = 2.0; hmm[1][1] = 3.0; hmm[1][2] = 10.0; hmm[2][0] = 2.0; hmm[2][1] = 6.0; hmm[2][2] = 1.0; hmm.source()[0] = vector(2.0, 1.0, 3.0); hmm.source()[1] = vector(1.0, 4.0, 3.0); hmm.source()[2] = vector(0.0, 5.0, 2.0); Info<< hmm << endl; //Info<< hmm.solve() << endl; //Info<< hmm << endl; //Info<< hmm.LUsolve() << endl; //Info<< hmm << endl; //Info<< "End\n" << endl; return 0; }``` What does the codes in red mean?Thanks in advance. This is the result in terminal: Code: ```3 3((-3 10 -4)(2 3 10)(2 6 1)) 3((2 1 3) (1 4 3) (0 5 2))```

 May 21, 2013, 10:18 #2 New Member   Pierre HORGUE Join Date: May 2009 Posts: 17 Rep Power: 8 Your system "hmm" is a linear system of vectors So you have A x = b with : - A the matrix - "x" the solution and "b" the source term (both are vector of vectors). The code in red fill the source vector "b" (composed by vector because your system is a simpleMatrix of vectors) Code: ``` hmm.source()[0] = vector(2.0, 1.0, 3.0); hmm.source()[1] = vector(1.0, 4.0, 3.0); hmm.source()[2] = vector(0.0, 5.0, 2.0);``` Pierre sharonyue likes this.

August 21, 2014, 23:46
#3
Senior Member

Wen Xu
Join Date: May 2014
Posts: 134
Rep Power: 3
Quote:
 Originally Posted by Pedro24 Your system "hmm" is a linear system of vectors So you have A x = b with : - A the matrix - "x" the solution and "b" the source term (both are vector of vectors). The code in red fill the source vector "b" (composed by vector because your system is a simpleMatrix of vectors) Code: ``` hmm.source()[0] = vector(2.0, 1.0, 3.0); hmm.source()[1] = vector(1.0, 4.0, 3.0); hmm.source()[2] = vector(0.0, 5.0, 2.0);``` Pierre
hi，Pedro24。but A=3 3((-3 10 -4)(2 3 10)(2 6 1)) multiply with x1=(-0.452599 1.12232 -0.149847) (the result got by the test code) do not equal to b1=(2 1 3). why? i have misunderstand this? please help me. Thank you!!!

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post holger_marschall OpenFOAM 242 March 7, 2013 13:30 cwang5 OpenFOAM Programming & Development 1 May 30, 2011 04:47 MechE OpenFOAM 28 May 16, 2011 11:02 Zdravko Stojanovic Main CFD Forum 2 July 19, 2010 10:11 Heinz Wilkening Main CFD Forum 38 March 5, 1999 12:44

All times are GMT -4. The time now is 17:49.