Re: Who uses clapack?

From: James Van Buskirk (not_valid_at_comcast.net)
Date: 12/13/04

  • Next message: kia: "Re: Who uses clapack?" 
    Date: Mon, 13 Dec 2004 04:36:37 GMT

    "Jentje Goslinga" <goslinga@telus.net> wrote in message
    news:41BD020D.8060503@telus.net...

    > Yes, I had understood that, but don't you have to write
    > separate implementations of those template functions for
    > complex numbers for example. OK, I understand now, the
    > complex float and complex double are the same code and the
    > real float and double are another code, understood.

    Now I don't think I get it. The Fortran version needs two
    separate codes to handle the 6 possibilities: a real code
    for single, double, and quad precision real and a complex
    code for single, double, and quad precision complex. But
    the previous poster seems to be saying that in the C++
    version, a single code suffices for the four cases of
    (single vs. double precision) and (real vs. complex).
    There are at least some level 1 blas subprograms where
    there are two complex versions and only one real version,
    e.g. cdotc/zdotc/ydotc & cdotu/zdotu/ydotu for complex
    data but only sdot/ddot/qdot for real data. It seems to
    me that you aren't quite going to get a 4X code reduction
    factor for several of the blas subprograms. 

    -- 
write(*,*) transfer((/17.392111325966148d0,6.5794487871554595D-85, &
6.0134700243160014d-154/),(/'x'/)); end
  • Next message: kia: "Re: Who uses clapack?"

    Relevant Pages

    • Re: Who uses clapack?
      ... > complex float and complex double are the same code and the ... code for single, double, and quad precision complex. ... data but only sdot/ddot/qdot for real data. ...
      (sci.math.num-analysis)
    • Re: Who uses clapack?
      ... > complex float and complex double are the same code and the ... code for single, double, and quad precision complex. ... data but only sdot/ddot/qdot for real data. ...
      (comp.lang.c)