On Fri, Jul 8, 2016 at 10:18 AM, David Knezevic <david.kneze...@akselos.com>
wrote:
>
> 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.
>
Sounds good.
I just double-checked and PETSc provides a LAPACKgesvd_ in the complex
numbers case, so you should still be able to go through their API.
--
John
------------------------------------------------------------------------------
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
