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

2

3
(2) 
4
(2) 
5
(6) 
6

7

8

9

10
(8) 
11
(1) 
12
(14) 
13
(10) 
14
(9) 
15

16
(2) 
17
(1) 
18
(7) 
19
(2) 
20
(6) 
21

22

23
(11) 
24

25

26
(7) 
27

28
(1) 
29
(2) 
30
(2) 
31




From: Roy Stogner <roystgnr@ic...>  20130723 21:14:08

On Tue, 23 Jul 2013, Roy Stogner wrote: > On Tue, 23 Jul 2013, Kathrin Smetana wrote: > >> I am using DenseVector<Number> and I would like to multiply it with >> a DenseMatrix<Real> to avoid unnecessary computations. > > This should actually be pretty quick to enable. Let me give it a try. Let me know if https://github.com/libMesh/libmesh/pull/114 works for you?  Roy 
From: Roy Stogner <roystgnr@ic...>  20130723 18:51:32

On Tue, 23 Jul 2013, Kathrin Smetana wrote: > I am using DenseVector<Number> and I would like to multiply it with > a DenseMatrix<Real> to avoid unnecessary computations. This should actually be pretty quick to enable. Let me give it a try.  Roy 
From: Kathrin Smetana <ksmetana@MIT.EDU>  20130723 18:39:52

Hi John, I am using DenseVector<Number> and I would like to multiply it with a DenseMatrix<Real> to avoid unnecessary computations. Best, Kathrin Am 23.07.2013 um 12:15 schrieb John Peterson <jwpeterson@...>: > > > > On Tue, Jul 23, 2013 at 10:10 AM, Kathrin Smetana <ksmetana@...> wrote: > Hi all, > > I am working on a elastodynamicsproblem with a complexvalued displacement field. Hence I want to multiply a dense matrix of type real with a vector of type number. As the DenseMatrix class only allows multiplication of a matrix and a vector of the same type, I was wondering if there might be a more elegant to do this than always changing the type of the matrix before applying the function vector_mult. > > More elegant to me would just be to use DenseVector<Number> instead of whatever you are currently using. > > Multiplication is defined between two members of R or two members of C, but not one of each :) > >  > John 
From: Kathrin Smetana <ksmetana@MIT.EDU>  20130723 17:48:49

Hi Manav, Many thanks for the files. I was also thinking of modifying the DenseMatrix class, but I was just wondering if there is a way to avoid that. Best, Kathrin On 07/23/2013 12:31 PM, Manav Bhatia wrote: > I have been doing this for a while using some modifications to the > dense_matrix class. > Attached are three files dense_matrix header and source, and the > dense_matrix_blas_lapack to facilitate the multiplication of real > matrix with complex vector. > Hope this helps. > Manav > > On Tue, Jul 23, 2013 at 12:15 PM, John Peterson <jwpeterson@... > <mailto:jwpeterson@...>> wrote: > > On Tue, Jul 23, 2013 at 10:10 AM, Kathrin Smetana > <ksmetana@... <mailto:ksmetana@...>> wrote: > > > Hi all, > > > > I am working on a elastodynamicsproblem with a complexvalued > > displacement field. Hence I want to multiply a dense matrix of > type real > > with a vector of type number. As the DenseMatrix class only allows > > multiplication of a matrix and a vector of the same type, I was > wondering > > if there might be a more elegant to do this than always changing > the type > > of the matrix before applying the function vector_mult. > > > > More elegant to me would just be to use DenseVector<Number> instead of > whatever you are currently using. > > Multiplication is defined between two members of R or two members > of C, but > not one of each :) > >  > John >  > See everything from the browser to the database with AppDynamics > Get endtoend visibility with application monitoring from AppDynamics > Isolate bottlenecks and diagnose root cause in seconds. > Start your free trial of AppDynamics Pro today! > http://pubads.g.doubleclick.net/gampad/clk?id=48808831&iu=/4140/ostg.clktrk > _______________________________________________ > Libmeshusers mailing list > Libmeshusers@... > <mailto:Libmeshusers@...> > https://lists.sourceforge.net/lists/listinfo/libmeshusers > > 
From: Manav Bhatia <bhatiamanav@gm...>  20130723 17:18:55

On Tue, Jul 23, 2013 at 12:57 PM, John Peterson <jwpeterson@...>wrote: > >>> Manav, >>> >>> There is some good stuff in here, particularly >>> the _nonhermitian_eig_lapack stuff. >>> >>> I'd like to get it into libmesh, do you have any test codes that verify >>> it is working? >>> >>>  >>> John >>> >> >> I am not sure if I have verified this code. I had started to add this a >> few months ago, but then it became unnecessary and I never got around to >> veryfying the code. >> I can test it at my end, and send you an updated patch if needed? >> > > That would be good, but it would be even better if we had a very short > libmesh code that also demonstrated it working (not a full application > code). > > It looks like it's designed to work when Number=Real, but I wonder if > you've tested that it compiles with Number=std::complex? > > Also, what do you envision the userfacing interface for _nonhermitian_eig_lapack > should be? For example, lu_solve, cholesky_solve, and svd call some > underlying Lapack routines. DenseMatrix::eig() maybe? > > >  > John > Good observation about the real/complex design. After looking at the code, I now realize that I never really instantiate the method (shows that I never got around to using it). I think that I was porting this implementation from some of my other codes, but stopped because I was able to get analytical expressions for the eigenvalues/eigenvectors of the matrix of interest. So, to answer your question, I have not verified that it works, and certainly not verified that it compiles with complex numbers. Both of these should not be difficult to do. A specialization to interface with complex lapack routines could also be added, but I have no experience with interfacing std::complex with the complex number representation for lapack. Any thoughts here? I like DenseMatrix::eig(). Perhps with arguments to return the eigenvalues and eigenvectors. I will try to put together this code and a small example and send it your way. Manav 
From: John Peterson <jwpeterson@gm...>  20130723 16:58:12

On Tue, Jul 23, 2013 at 10:47 AM, Manav Bhatia <bhatiamanav@...>wrote: > > > On Tue, Jul 23, 2013 at 12:42 PM, John Peterson <jwpeterson@...>wrote: > >> >> Manav, >> >> There is some good stuff in here, particularly >> the _nonhermitian_eig_lapack stuff. >> >> I'd like to get it into libmesh, do you have any test codes that verify >> it is working? >> >>  >> John >> > > I am not sure if I have verified this code. I had started to add this a > few months ago, but then it became unnecessary and I never got around to > veryfying the code. > I can test it at my end, and send you an updated patch if needed? > That would be good, but it would be even better if we had a very short libmesh code that also demonstrated it working (not a full application code). It looks like it's designed to work when Number=Real, but I wonder if you've tested that it compiles with Number=std::complex? Also, what do you envision the userfacing interface for _nonhermitian_eig_lapack should be? For example, lu_solve, cholesky_solve, and svd call some underlying Lapack routines. DenseMatrix::eig() maybe?  John 
From: Manav Bhatia <bhatiamanav@gm...>  20130723 16:47:34

On Tue, Jul 23, 2013 at 12:42 PM, John Peterson <jwpeterson@...>wrote: > > Manav, > > There is some good stuff in here, particularly > the _nonhermitian_eig_lapack stuff. > > I'd like to get it into libmesh, do you have any test codes that verify it > is working? > >  > John > I am not sure if I have verified this code. I had started to add this a few months ago, but then it became unnecessary and I never got around to veryfying the code. I can test it at my end, and send you an updated patch if needed? Manav 
From: John Peterson <jwpeterson@gm...>  20130723 16:42:41

On Tue, Jul 23, 2013 at 10:31 AM, Manav Bhatia <bhatiamanav@...>wrote: > I have been doing this for a while using some modifications to the > dense_matrix class. > > Attached are three files dense_matrix header and source, and the > dense_matrix_blas_lapack to facilitate the multiplication of real matrix > with complex vector. > > Hope this helps. > Manav, There is some good stuff in here, particularly the _nonhermitian_eig_lapack stuff. I'd like to get it into libmesh, do you have any test codes that verify it is working?  John 
From: Manav Bhatia <bhatiamanav@gm...>  20130723 16:31:12

I have been doing this for a while using some modifications to the dense_matrix class. Attached are three files dense_matrix header and source, and the dense_matrix_blas_lapack to facilitate the multiplication of real matrix with complex vector. Hope this helps. Manav On Tue, Jul 23, 2013 at 12:15 PM, John Peterson <jwpeterson@...>wrote: > On Tue, Jul 23, 2013 at 10:10 AM, Kathrin Smetana <ksmetana@...> > wrote: > > > Hi all, > > > > I am working on a elastodynamicsproblem with a complexvalued > > displacement field. Hence I want to multiply a dense matrix of type real > > with a vector of type number. As the DenseMatrix class only allows > > multiplication of a matrix and a vector of the same type, I was wondering > > if there might be a more elegant to do this than always changing the type > > of the matrix before applying the function vector_mult. > > > > More elegant to me would just be to use DenseVector<Number> instead of > whatever you are currently using. > > Multiplication is defined between two members of R or two members of C, but > not one of each :) > >  > John > >  > See everything from the browser to the database with AppDynamics > Get endtoend visibility with application monitoring from AppDynamics > Isolate bottlenecks and diagnose root cause in seconds. > Start your free trial of AppDynamics Pro today! > http://pubads.g.doubleclick.net/gampad/clk?id=48808831&iu=/4140/ostg.clktrk > _______________________________________________ > Libmeshusers mailing list > Libmeshusers@... > https://lists.sourceforge.net/lists/listinfo/libmeshusers > 
From: John Peterson <jwpeterson@gm...>  20130723 16:16:06

On Tue, Jul 23, 2013 at 10:10 AM, Kathrin Smetana <ksmetana@...> wrote: > Hi all, > > I am working on a elastodynamicsproblem with a complexvalued > displacement field. Hence I want to multiply a dense matrix of type real > with a vector of type number. As the DenseMatrix class only allows > multiplication of a matrix and a vector of the same type, I was wondering > if there might be a more elegant to do this than always changing the type > of the matrix before applying the function vector_mult. > More elegant to me would just be to use DenseVector<Number> instead of whatever you are currently using. Multiplication is defined between two members of R or two members of C, but not one of each :)  John 
From: Kathrin Smetana <ksmetana@MIT.EDU>  20130723 16:10:16

Hi all, I am working on a elastodynamicsproblem with a complexvalued displacement field. Hence I want to multiply a dense matrix of type real with a vector of type number. As the DenseMatrix class only allows multiplication of a matrix and a vector of the same type, I was wondering if there might be a more elegant to do this than always changing the type of the matrix before applying the function vector_mult. Thanks, Kathrin 