You can subscribe to this list here.
2003 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}
(2) 
_{Oct}
(2) 
_{Nov}
(27) 
_{Dec}
(31) 

2004 
_{Jan}
(6) 
_{Feb}
(15) 
_{Mar}
(33) 
_{Apr}
(10) 
_{May}
(46) 
_{Jun}
(11) 
_{Jul}
(21) 
_{Aug}
(15) 
_{Sep}
(13) 
_{Oct}
(23) 
_{Nov}
(1) 
_{Dec}
(8) 
2005 
_{Jan}
(27) 
_{Feb}
(57) 
_{Mar}
(86) 
_{Apr}
(23) 
_{May}
(37) 
_{Jun}
(34) 
_{Jul}
(24) 
_{Aug}
(17) 
_{Sep}
(50) 
_{Oct}
(24) 
_{Nov}
(10) 
_{Dec}
(60) 
2006 
_{Jan}
(47) 
_{Feb}
(46) 
_{Mar}
(127) 
_{Apr}
(19) 
_{May}
(26) 
_{Jun}
(62) 
_{Jul}
(47) 
_{Aug}
(51) 
_{Sep}
(61) 
_{Oct}
(42) 
_{Nov}
(50) 
_{Dec}
(33) 
2007 
_{Jan}
(60) 
_{Feb}
(55) 
_{Mar}
(77) 
_{Apr}
(102) 
_{May}
(82) 
_{Jun}
(102) 
_{Jul}
(169) 
_{Aug}
(117) 
_{Sep}
(80) 
_{Oct}
(37) 
_{Nov}
(51) 
_{Dec}
(43) 
2008 
_{Jan}
(71) 
_{Feb}
(94) 
_{Mar}
(98) 
_{Apr}
(125) 
_{May}
(54) 
_{Jun}
(119) 
_{Jul}
(60) 
_{Aug}
(111) 
_{Sep}
(118) 
_{Oct}
(125) 
_{Nov}
(119) 
_{Dec}
(94) 
2009 
_{Jan}
(109) 
_{Feb}
(38) 
_{Mar}
(93) 
_{Apr}
(88) 
_{May}
(29) 
_{Jun}
(57) 
_{Jul}
(53) 
_{Aug}
(48) 
_{Sep}
(68) 
_{Oct}
(151) 
_{Nov}
(23) 
_{Dec}
(35) 
2010 
_{Jan}
(84) 
_{Feb}
(60) 
_{Mar}
(184) 
_{Apr}
(112) 
_{May}
(60) 
_{Jun}
(90) 
_{Jul}
(23) 
_{Aug}
(70) 
_{Sep}
(119) 
_{Oct}
(27) 
_{Nov}
(47) 
_{Dec}
(54) 
2011 
_{Jan}
(22) 
_{Feb}
(19) 
_{Mar}
(92) 
_{Apr}
(93) 
_{May}
(35) 
_{Jun}
(91) 
_{Jul}
(32) 
_{Aug}
(61) 
_{Sep}
(7) 
_{Oct}
(69) 
_{Nov}
(81) 
_{Dec}
(23) 
2012 
_{Jan}
(64) 
_{Feb}
(95) 
_{Mar}
(35) 
_{Apr}
(36) 
_{May}
(63) 
_{Jun}
(98) 
_{Jul}
(70) 
_{Aug}
(171) 
_{Sep}
(149) 
_{Oct}
(64) 
_{Nov}
(67) 
_{Dec}
(126) 
2013 
_{Jan}
(108) 
_{Feb}
(104) 
_{Mar}
(171) 
_{Apr}
(133) 
_{May}
(108) 
_{Jun}
(100) 
_{Jul}
(93) 
_{Aug}
(126) 
_{Sep}
(74) 
_{Oct}
(59) 
_{Nov}
(145) 
_{Dec}
(93) 
2014 
_{Jan}
(38) 
_{Feb}
(45) 
_{Mar}
(26) 
_{Apr}
(41) 
_{May}
(125) 
_{Jun}
(70) 
_{Jul}
(28) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 






1
(1) 
2

3
(2) 
4

5
(3) 
6
(3) 
7
(13) 
8
(4) 
9

10

11

12
(3) 
13
(1) 
14

15
(7) 
16

17

18

19
(5) 
20
(2) 
21
(1) 
22
(1) 
23
(1) 
24

25
(5) 
26

27

28
(5) 
29
(4) 
30

From: HaeWon Choi <haewon@uc...>  20060929 14:53:48

Hi, actually John clarify what I exactly want. Thank you John, For PetSc, MatZeroRows set 1 for all diagonal entities of all Dirichlet boundary rows (other entities will be zero) as I =20 mentioned last time. In fact you can find examples using this approach from PetSc since I =20 have learned this from PetSc. What you have to do is modify RHS as John's example shows. I have used this method for my other codes using PetSc (but not for =20 libMesh yet). This approach gives same results as reduced matrix method. HaeWon On Sep 29, 2006, at 5:45 AM, John Peterson wrote: > Just as an addendum to Tim's note, you can maintain any > symmetry originally present in the problem by "subtracting" > the column entries multiplied by the Dirichlet value from > the right hand side vector. > > If Au=3Db, where > > [a_11 a_12 a_13] [u1] =3D [b1] > [a_21 a_22 a_23] [u2] =3D [b2] > [a_31 a_32 a_33] [u3] =3D [b3] > > and u1 =3D g1, a nonhomogeneous BC val, we can modify Au=3Db as: > > [ 1 0 0 ] [u1] =3D [g1] > [ 0 a_22 a_23] [u2] =3D [b2  g1*a21] > [ 0 a_32 a_33] [u3] =3D [b3  g1*a31] > > > This imposes u1=3Dg1, and maintains any original symmetry of A. > > > John > > > > Tim Kr=F6ger writes: >> Dear all, >> >> On Fri, 29 Sep 2006, li pan wrote: >> >>> I'm also thinking about this. Are you sure that you >>> only need to zero all the enties of rows? I read >>> somewhere that columes should also be zeroed. Could >>> somebody confirm this? >> >> If you zero only the row entries, the matrix will no longer be >> symmetric. This is often considered as a drawback (provided that the >> matrix was symmetric before) because certain solvers (e.g. CG) cannot >> be used then any more. >> >> On the other hand, if the column entries are also zeroed, the =20 >> solution >> of the system will be wrong  except for the case of *homogeneous* >> Dirichlet conditions. For this reason, some people transform their >> problems to homogeneous boundary conditions, i.e. they do in practice >> the same thing that is usually done theoretically anyway, i.e. they >> subtract a function that fulfills all Dirichlet boundary conditions >> but not the PDE. Note that in the discrete case, such a function is >> trivial to find. >> >> Because I find this all quite unsatisfactory, I was glad to see that >> libMesh uses the penalty method (which I did not know before) because >> it is easy to implement, works in all cases, and does not destroy >> symmetry or positivity of the matrix. >> >> Best Regards, >> >> Tim > > =20= >  > Take Surveys. Earn Cash. Influence the Future of IT > Join SourceForge.net's Techsay panel and you'll get the chance to =20 > share your > opinions on IT & business topics through brief surveys  and earn =20 > cash > http://www.techsay.com/default.php?=20 > page=3Djoin.php&p=3Dsourceforge&CID=3DDEVDEV____________________________= ____=20 > _______________ > Libmeshusers mailing list > Libmeshusers@... > https://lists.sourceforge.net/lists/listinfo/libmeshusers 
From: John Peterson <peterson@cf...>  20060929 11:45:34

Just as an addendum to Tim's note, you can maintain any symmetry originally present in the problem by "subtracting" the column entries multiplied by the Dirichlet value from the right hand side vector. If Au=3Db, where [a_11 a_12 a_13] [u1] =3D [b1] [a_21 a_22 a_23] [u2] =3D [b2] [a_31 a_32 a_33] [u3] =3D [b3] and u1 =3D g1, a nonhomogeneous BC val, we can modify Au=3Db as: [ 1 0 0 ] [u1] =3D [g1] [ 0 a_22 a_23] [u2] =3D [b2  g1*a21] [ 0 a_32 a_33] [u3] =3D [b3  g1*a31] This imposes u1=3Dg1, and maintains any original symmetry of A. John Tim Kr=F6ger writes: > Dear all, >=20 > On Fri, 29 Sep 2006, li pan wrote: >=20 > > I'm also thinking about this. Are you sure that you > > only need to zero all the enties of rows? I read > > somewhere that columes should also be zeroed. Could > > somebody confirm this? >=20 > If you zero only the row entries, the matrix will no longer be=20 > symmetric. This is often considered as a drawback (provided that th= e=20 > matrix was symmetric before) because certain solvers (e.g. CG) canno= t=20 > be used then any more. >=20 > On the other hand, if the column entries are also zeroed, the soluti= on=20 > of the system will be wrong  except for the case of *homogeneous*=20= > Dirichlet conditions. For this reason, some people transform their=20= > problems to homogeneous boundary conditions, i.e. they do in practic= e=20 > the same thing that is usually done theoretically anyway, i.e. they=20= > subtract a function that fulfills all Dirichlet boundary conditions=20= > but not the PDE. Note that in the discrete case, such a function is= =20 > trivial to find. >=20 > Because I find this all quite unsatisfactory, I was glad to see that= =20 > libMesh uses the penalty method (which I did not know before) becaus= e=20 > it is easy to implement, works in all cases, and does not destroy=20= > symmetry or positivity of the matrix. >=20 > Best Regards, >=20 > Tim 
From: <tim@ce...>  20060929 08:08:17

Dear all, On Fri, 29 Sep 2006, li pan wrote: > I'm also thinking about this. Are you sure that you > only need to zero all the enties of rows? I read > somewhere that columes should also be zeroed. Could > somebody confirm this? If you zero only the row entries, the matrix will no longer be symmetric. This is often considered as a drawback (provided that the matrix was symmetric before) because certain solvers (e.g. CG) cannot be used then any more. On the other hand, if the column entries are also zeroed, the solution of the system will be wrong  except for the case of *homogeneous* Dirichlet conditions. For this reason, some people transform their problems to homogeneous boundary conditions, i.e. they do in practice the same thing that is usually done theoretically anyway, i.e. they subtract a function that fulfills all Dirichlet boundary conditions but not the PDE. Note that in the discrete case, such a function is trivial to find. Because I find this all quite unsatisfactory, I was glad to see that libMesh uses the penalty method (which I did not know before) because it is easy to implement, works in all cases, and does not destroy symmetry or positivity of the matrix. Best Regards, Tim 
From: li pan <li76pan@ya...>  20060929 07:45:27

hi, I'm also thinking about this. Are you sure that you only need to zero all the enties of rows? I read somewhere that columes should also be zeroed. Could somebody confirm this? pan  HaeWon Choi <haewon@...> wrote: > Hi, all: > I'm wondering why libmesh does not implement > MatZeroRows option > given in PetSc. > In fact this B.C. treatment is more appropriate > since you only need > to set > diagonal entities as 1.0 for all Dirichlet boundary > rows. > Of course, in this case you need to add RHS > with additional boundary contributions from LHS. > So implementation is exactly same when you deal with > reduced matrix > boundary implementation but instead of physically > remove boundary > nodes you need to add MatZeroRows option after > MatAssembly. > I would like to have this option in libmesh. > > Thanks. > > >  > Take Surveys. Earn Cash. Influence the Future of IT > Join SourceForge.net's Techsay panel and you'll get > the chance to share your > opinions on IT & business topics through brief > surveys  and earn cash > http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV > _______________________________________________ > Libmeshusers mailing list > Libmeshusers@... > https://lists.sourceforge.net/lists/listinfo/libmeshusers > __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com 