00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016
00017
00018 static integer c__7 = 7;
00019 static real c_b5 = 0.f;
00020 static real c_b6 = 1.f;
00021
00022 doublereal stzt02_(integer *m, integer *n, real *af, integer *lda, real *tau,
00023 real *work, integer *lwork)
00024 {
00025
00026 integer af_dim1, af_offset, i__1, i__2;
00027 real ret_val;
00028
00029
00030 integer i__;
00031 real rwork[1];
00032 extern doublereal slamch_(char *), slange_(char *, integer *,
00033 integer *, real *, integer *, real *);
00034 extern int xerbla_(char *, integer *), slaset_(
00035 char *, integer *, integer *, real *, real *, real *, integer *), slatzm_(char *, integer *, integer *, real *, integer *,
00036 real *, real *, real *, integer *, real *);
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097 af_dim1 = *lda;
00098 af_offset = 1 + af_dim1;
00099 af -= af_offset;
00100 --tau;
00101 --work;
00102
00103
00104 ret_val = 0.f;
00105
00106 if (*lwork < *n * *n + *n) {
00107 xerbla_("STZT02", &c__7);
00108 return ret_val;
00109 }
00110
00111
00112
00113 if (*m <= 0 || *n <= 0) {
00114 return ret_val;
00115 }
00116
00117
00118
00119 slaset_("Full", n, n, &c_b5, &c_b6, &work[1], n);
00120
00121
00122
00123 for (i__ = *m; i__ >= 1; --i__) {
00124 i__1 = *n - *m + 1;
00125 slatzm_("Left", &i__1, n, &af[i__ + (*m + 1) * af_dim1], lda, &tau[
00126 i__], &work[i__], &work[*m + 1], n, &work[*n * *n + 1]);
00127
00128 }
00129
00130
00131
00132 i__1 = *m;
00133 for (i__ = 1; i__ <= i__1; ++i__) {
00134 i__2 = *n - *m + 1;
00135 slatzm_("Left", &i__2, n, &af[i__ + (*m + 1) * af_dim1], lda, &tau[
00136 i__], &work[i__], &work[*m + 1], n, &work[*n * *n + 1]);
00137
00138 }
00139
00140
00141
00142 i__1 = *n;
00143 for (i__ = 1; i__ <= i__1; ++i__) {
00144 work[(i__ - 1) * *n + i__] += -1.f;
00145
00146 }
00147
00148 ret_val = slange_("One-norm", n, n, &work[1], n, rwork) / (
00149 slamch_("Epsilon") * (real) max(*m,*n));
00150 return ret_val;
00151
00152
00153
00154 }