add missing files
authorGuillaume Horel <guillaume.horel@gmail.com>
Tue, 3 Sep 2019 18:45:43 +0000 (14:45 -0400)
committerGuillaume Horel <guillaume.horel@gmail.com>
Sun, 8 Sep 2019 15:14:49 +0000 (11:14 -0400)
lapack/trtrs/Makefile
lapack/trtrs/trtrs_parallel.c [new file with mode: 0644]
lapack/trtrs/trtrs_single.c
lapack/trtrs/ztrtrs_parallel.c [new file with mode: 0644]
lapack/trtrs/ztrtrs_single.c [new file with mode: 0644]

index 400b8b6..f9faaf9 100644 (file)
@@ -305,4 +305,148 @@ xtrtrs_LCU_single.$(SUFFIX) : xtrtrs_single.c
 xtrtrs_LCN_single.$(SUFFIX) : xtrtrs_single.c
        $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=4 -DDIAG $< -o $(@F)
 
+ctrtrs_UNU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+ctrtrs_UNN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+ctrtrs_UTU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+ctrtrs_UTN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+ctrtrs_URU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+ctrtrs_URN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+ctrtrs_UCU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+ctrtrs_UCN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -UUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
+ctrtrs_LNU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+ctrtrs_LNN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+ctrtrs_LTU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+ctrtrs_LTN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+ctrtrs_LRU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+ctrtrs_LRN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+ctrtrs_LCU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+ctrtrs_LCN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -UDOUBLE -DUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
+ztrtrs_UNU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+ztrtrs_UNN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+ztrtrs_UTU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+ztrtrs_UTN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+ztrtrs_URU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+ztrtrs_URN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+ztrtrs_UCU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+ztrtrs_UCN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -UUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
+ztrtrs_LNU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+ztrtrs_LNN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+ztrtrs_LTU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+ztrtrs_LTN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+ztrtrs_LRU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+ztrtrs_LRN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+ztrtrs_LCU_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+ztrtrs_LCN_parallel.$(SUFFIX) : ztrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DDOUBLE -DUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
+xtrtrs_UNU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+xtrtrs_UNN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+xtrtrs_UTU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+xtrtrs_UTN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+xtrtrs_URU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+xtrtrs_URN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+xtrtrs_UCU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+xtrtrs_UCN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -UUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
+xtrtrs_LNU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=1 -UDIAG $< -o $(@F)
+
+xtrtrs_LNN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=1 -DDIAG $< -o $(@F)
+
+xtrtrs_LTU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=2 -UDIAG $< -o $(@F)
+
+xtrtrs_LTN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=2 -DDIAG $< -o $(@F)
+
+xtrtrs_LRU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=3 -UDIAG $< -o $(@F)
+
+xtrtrs_LRN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=3 -DDIAG $< -o $(@F)
+
+xtrtrs_LCU_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=4 -UDIAG $< -o $(@F)
+
+xtrtrs_LCN_parallel.$(SUFFIX) : xtrtrs_parallel.c
+       $(CC) -c $(CFLAGS) -DCOMPLEX -DXDOUBLE -DUPLO -DTRANS=4 -DDIAG $< -o $(@F)
+
 include ../../Makefile.tail
diff --git a/lapack/trtrs/trtrs_parallel.c b/lapack/trtrs/trtrs_parallel.c
new file mode 100644 (file)
index 0000000..52f42f6
--- /dev/null
@@ -0,0 +1,111 @@
+/*********************************************************************/
+/* Copyright 2009, 2010 The University of Texas at Austin.           */
+/* All rights reserved.                                              */
+/*                                                                   */
+/* Redistribution and use in source and binary forms, with or        */
+/* without modification, are permitted provided that the following   */
+/* conditions are met:                                               */
+/*                                                                   */
+/*   1. Redistributions of source code must retain the above         */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer.                                                  */
+/*                                                                   */
+/*   2. Redistributions in binary form must reproduce the above      */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer in the documentation and/or other materials       */
+/*      provided with the distribution.                              */
+/*                                                                   */
+/*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */
+/*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */
+/*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */
+/*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */
+/*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */
+/*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */
+/*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */
+/*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */
+/*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */
+/*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */
+/*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */
+/*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */
+/*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */
+/*    POSSIBILITY OF SUCH DAMAGE.                                    */
+/*                                                                   */
+/* The views and conclusions contained in the software and           */
+/* documentation are those of the authors and should not be          */
+/* interpreted as representing official policies, either expressed   */
+/* or implied, of The University of Texas at Austin.                 */
+/*********************************************************************/
+
+#include <stdio.h>
+#include "common.h"
+
+#if   !defined(TRANS) && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNUU
+#define TRSV TRSV_NUU
+#elif !defined(TRANS) && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNUN
+#define TRSV TRSV_NUN
+#elif !defined(TRANS) && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNLU
+#define TRSV TRSV_NLU
+#elif !defined(TRANS) && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNLN
+#define TRSV TRSV_NLN
+#elif defined(TRANS) && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTUU
+#define TRSV TRSV_TUU
+#elif defined(TRANS) && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTUN
+#define TRSV TRSV_TUN
+#elif defined(TRANS) && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTLU
+#define TRSV TRSV_TLU
+#elif defined(TRANS) && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTLN
+#define TRSV TRSV_TLN
+#endif
+
+static int inner_thread(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n,
+                        FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
+
+  TRSM (args, range_m, range_n, sa, sb, 0);
+
+  return 0;
+}
+
+blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
+
+  int mode;
+
+#ifndef TRANS
+    if (args -> n  == 1){
+      TRSV (args -> m, args -> a, args -> lda, args -> b, 1, sb);
+    } else {
+#ifdef XDOUBLE
+      mode  =  BLAS_XDOUBLE | BLAS_REAL;
+#elif defined(DOUBLE)
+      mode  =  BLAS_DOUBLE  | BLAS_REAL;
+#else
+      mode  =  BLAS_SINGLE  | BLAS_REAL;
+#endif
+
+      gemm_thread_n(mode, args, NULL, NULL, inner_thread, sa, sb,  args -> nthreads);
+    }
+#else
+    if (args -> n == 1){
+      TRSV (args -> m, args -> a, args -> lda, args -> b, 1, sb);
+    } else {
+#ifdef XDOUBLE
+      mode  =  BLAS_XDOUBLE | BLAS_REAL | (1 << BLAS_TRANSA_SHIFT);
+#elif defined(DOUBLE)
+      mode  =  BLAS_DOUBLE  | BLAS_REAL | (1 << BLAS_TRANSA_SHIFT);
+#else
+      mode  =  BLAS_SINGLE  | BLAS_REAL | (1 << BLAS_TRANSA_SHIFT);
+#endif
+
+      gemm_thread_n(mode, args, NULL, NULL, inner_thread, sa, sb,  args -> nthreads);
+    }
+#endif
+
+    return 0;
+  }
index a690d4a..c82b813 100644 (file)
 #include <stdio.h>
 #include "common.h"
 
+#if   !defined(TRANS) && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNUU
+#define TRSV TRSV_NUU
+#elif !defined(TRANS) && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNUN
+#define TRSV TRSV_NUN
+#elif !defined(TRANS) && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNLU
+#define TRSV TRSV_NLU
+#elif !defined(TRANS) && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNLN
+#define TRSV TRSV_NLN
+#elif defined(TRANS) && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTUU
+#define TRSV TRSV_TUU
+#elif defined(TRANS) && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTUN
+#define TRSV TRSV_TUN
+#elif defined(TRANS) && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTLU
+#define TRSV TRSV_TLU
+#elif defined(TRANS) && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTLN
+#define TRSV TRSV_TLN
+#endif
+
 blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
 
-#ifndef UPLO
-#ifndef DIAG
-#ifndef DIAG
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#else
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#endif
-#else
-#ifndef DIAG
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#else
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#endif
-#else
-#ifndef DIAG
-#ifndef DIAG
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#else
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#endif
-#else
-#ifndef DIAG
     if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
+        TRSV (args -> m, args -> a, args -> lda, args -> b, 1, sb);
     } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
+        TRSM (args, range_m, range_n, sa, sb, 0);
     }
-#else
-    if (args -> n == 1){
-        TRSV_NLU (args -> m, args -> a, args -> lda, args -> b, 1, sb);
-    } else {
-        TRSM_LNLU (args, range_m, range_n, sa, sb, 0);
-    }
-#endif
-#endif
   return 0;  }
diff --git a/lapack/trtrs/ztrtrs_parallel.c b/lapack/trtrs/ztrtrs_parallel.c
new file mode 100644 (file)
index 0000000..d5248f2
--- /dev/null
@@ -0,0 +1,118 @@
+/*********************************************************************/
+/* Copyright 2009, 2010 The University of Texas at Austin.           */
+/* All rights reserved.                                              */
+/*                                                                   */
+/* Redistribution and use in source and binary forms, with or        */
+/* without modification, are permitted provided that the following   */
+/* conditions are met:                                               */
+/*                                                                   */
+/*   1. Redistributions of source code must retain the above         */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer.                                                  */
+/*                                                                   */
+/*   2. Redistributions in binary form must reproduce the above      */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer in the documentation and/or other materials       */
+/*      provided with the distribution.                              */
+/*                                                                   */
+/*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */
+/*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */
+/*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */
+/*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */
+/*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */
+/*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */
+/*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */
+/*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */
+/*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */
+/*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */
+/*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */
+/*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */
+/*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */
+/*    POSSIBILITY OF SUCH DAMAGE.                                    */
+/*                                                                   */
+/* The views and conclusions contained in the software and           */
+/* documentation are those of the authors and should not be          */
+/* interpreted as representing official policies, either expressed   */
+/* or implied, of The University of Texas at Austin.                 */
+/*********************************************************************/
+
+#include <stdio.h>
+#include "common.h"
+
+#if   TRANS == 1 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNUU
+#define ZTRSV ZTRSV_NUU
+#elif TRANS == 1 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNUN
+#define ZTRSV ZTRSV_NUN
+#elif TRANS == 1 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNLU
+#define ZTRSV ZTRSV_NLU
+#elif TRANS == 1 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNLN
+#define ZTRSV ZTRSV_NLN
+#elif TRANS == 2 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTUU
+#define ZTRSV ZTRSV_TUU
+#elif TRANS == 2 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTUN
+#define ZTRSV ZTRSV_TUN
+#elif TRANS == 2 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTLU
+#define ZTRSV ZTRSV_TLU
+#elif TRANS == 2 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTLN
+#define ZTRSV ZTRSV_TLN
+#elif TRANS == 3 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LRUU
+#define ZTRSV ZTRSV_RUU
+#elif TRANS == 3 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LRUN
+#define ZTRSV ZTRSV_RUN
+#elif TRANS == 3 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LRLU
+#define ZTRSV ZTRSV_RLU
+#elif TRANS == 3 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LRLN
+#define ZTRSV ZTRSV_RLN
+#elif TRANS == 4 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LCUU
+#define ZTRSV ZTRSV_CUU
+#elif TRANS == 4 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LCUN
+#define ZTRSV ZTRSV_CUN
+#elif TRANS == 4 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LCLU
+#define ZTRSV ZTRSV_CLU
+#elif TRANS == 4 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LCLN
+#define ZTRSV ZTRSV_CLN
+#endif
+
+static int inner_thread(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n,
+                        FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
+
+  TRSM (args, range_m, range_n, sa, sb, 0);
+  return 0;
+}
+
+blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
+
+  int mode;
+
+    if (args -> n == 1){
+      ZTRSV (args -> m, args -> a, args -> lda, args -> b, 1, sb);
+    } else {
+#ifdef XDOUBLE
+      mode  =  BLAS_XDOUBLE | BLAS_COMPLEX;
+#elif defined(DOUBLE)
+      mode  =  BLAS_DOUBLE  | BLAS_COMPLEX;
+#else
+      mode  =  BLAS_SINGLE  | BLAS_COMPLEX;
+#endif
+
+      gemm_thread_n(mode, args, NULL, NULL, inner_thread, sa, sb, args -> nthreads);
+    }
+
+   return 0;
+  }
diff --git a/lapack/trtrs/ztrtrs_single.c b/lapack/trtrs/ztrtrs_single.c
new file mode 100644 (file)
index 0000000..f39d729
--- /dev/null
@@ -0,0 +1,98 @@
+/*********************************************************************/
+/* Copyright 2009, 2010 The University of Texas at Austin.           */
+/* All rights reserved.                                              */
+/*                                                                   */
+/* Redistribution and use in source and binary forms, with or        */
+/* without modification, are permitted provided that the following   */
+/* conditions are met:                                               */
+/*                                                                   */
+/*   1. Redistributions of source code must retain the above         */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer.                                                  */
+/*                                                                   */
+/*   2. Redistributions in binary form must reproduce the above      */
+/*      copyright notice, this list of conditions and the following  */
+/*      disclaimer in the documentation and/or other materials       */
+/*      provided with the distribution.                              */
+/*                                                                   */
+/*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */
+/*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */
+/*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */
+/*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */
+/*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */
+/*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */
+/*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */
+/*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */
+/*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */
+/*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */
+/*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */
+/*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */
+/*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */
+/*    POSSIBILITY OF SUCH DAMAGE.                                    */
+/*                                                                   */
+/* The views and conclusions contained in the software and           */
+/* documentation are those of the authors and should not be          */
+/* interpreted as representing official policies, either expressed   */
+/* or implied, of The University of Texas at Austin.                 */
+/*********************************************************************/
+
+#include <stdio.h>
+#include "common.h"
+
+#if   TRANS == 1 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNUU
+#define ZTRSV ZTRSV_NUU
+#elif TRANS == 1 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNUN
+#define ZTRSV ZTRSV_NUN
+#elif TRANS == 1 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LNLU
+#define ZTRSV ZTRSV_NLU
+#elif TRANS == 1 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LNLN
+#define ZTRSV ZTRSV_NLN
+#elif TRANS == 2 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTUU
+#define ZTRSV ZTRSV_TUU
+#elif TRANS == 2 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTUN
+#define ZTRSV ZTRSV_TUN
+#elif TRANS == 2 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LTLU
+#define ZTRSV ZTRSV_TLU
+#elif TRANS == 2 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LTLN
+#define ZTRSV ZTRSV_TLN
+#elif TRANS == 3 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LRUU
+#define ZTRSV ZTRSV_RUU
+#elif TRANS == 3 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LRUN
+#define ZTRSV ZTRSV_RUN
+#elif TRANS == 3 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LRLU
+#define ZTRSV ZTRSV_RLU
+#elif TRANS == 3 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LRLN
+#define ZTRSV ZTRSV_RLN
+#elif TRANS == 4 && !defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LCUU
+#define ZTRSV ZTRSV_CUU
+#elif TRANS == 4 && !defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LCUN
+#define ZTRSV ZTRSV_CUN
+#elif TRANS == 4 && defined(UPLO) && !defined(DIAG)
+#define TRSM TRSM_LCLU
+#define ZTRSV ZTRSV_CLU
+#elif TRANS == 4 && defined(UPLO) && defined(DIAG)
+#define TRSM TRSM_LCLN
+#define ZTRSV ZTRSV_CLN
+#endif
+
+blasint CNAME(blas_arg_t *args, BLASLONG *range_m, BLASLONG *range_n, FLOAT *sa, FLOAT *sb, BLASLONG mypos) {
+    if (args -> n == 1){
+        ZTRSV (args -> m, args -> a, args -> lda, args -> b, 1, sb);
+    } else {
+        TRSM (args, range_m, range_n, sa, sb, 0);
+    }
+  return 0;  }