On Fri, Jul 8, 2016 at 12:10 PM, John Peterson <jwpeter...@gmail.com> wrote:
>
>
> On Fri, Jul 8, 2016 at 9:28 AM, David Knezevic <david.kneze...@akselos.com
> > wrote:
>
>> _svd_helper (and a few other functions) in dense_matrix_blas_lapack.C are
>> only available with LIBMESH_USE_REAL_NUMBERS. I'd like to use the SVD in
>> the complex-valued case.
>>
>> Any idea if there is a good reason for restricting to the real case, or
>> was it just that no one has needed the complex-valued case so far?
>>
>
> You added the _svd_lapack code with the real number restrictions baked in
> back in 2011:
>
> 39db53db94f54ab416a40da1db86f2f4f6b2fe4b
>
> I guess only 2011 David can answer this question ;-)
>
haha, oops :) Though the same question applies to all the _lapack stuff in
that file, not just the svd case.
In any case, I've looked at it a bit more closely. The lapack API has a
different number of arguments in the real vs. complex cases, which I think
is why we never bothered to support complex (it requires various extra
#ifdefs). I'll see if I can at least get it working for the svd case.
David
------------------------------------------------------------------------------
Attend Shape: An AT&T Tech Expo July 15-16. Meet us at AT&T Park in San
Francisco, CA to explore cutting-edge tech and listen to tech luminaries
present their vision of the future. This family event has something for
everyone, including kids. Get more information and register today.
http://sdm.link/attshape
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
