Doxygen Source Code Documentation
Main Page Alphabetical List Data Structures File List Data Fields Globals Search
eis_cbabk2.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006 #include "f2c.h"
00007
00008 int cbabk2_(integer *nm, integer *n, integer *low, integer *
00009 igh, doublereal *scale, integer *m, doublereal *zr, doublereal *zi)
00010 {
00011
00012 integer zr_dim1, zr_offset, zi_dim1, zi_offset, i__1, i__2;
00013
00014
00015 static integer i__, j, k;
00016 static doublereal s;
00017 static integer ii;
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
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 --scale;
00066 zi_dim1 = *nm;
00067 zi_offset = zi_dim1 + 1;
00068 zi -= zi_offset;
00069 zr_dim1 = *nm;
00070 zr_offset = zr_dim1 + 1;
00071 zr -= zr_offset;
00072
00073
00074 if (*m == 0) {
00075 goto L200;
00076 }
00077 if (*igh == *low) {
00078 goto L120;
00079 }
00080
00081 i__1 = *igh;
00082 for (i__ = *low; i__ <= i__1; ++i__) {
00083 s = scale[i__];
00084
00085
00086
00087 i__2 = *m;
00088 for (j = 1; j <= i__2; ++j) {
00089 zr[i__ + j * zr_dim1] *= s;
00090 zi[i__ + j * zi_dim1] *= s;
00091
00092 }
00093
00094
00095 }
00096
00097
00098 L120:
00099 i__1 = *n;
00100 for (ii = 1; ii <= i__1; ++ii) {
00101 i__ = ii;
00102 if (i__ >= *low && i__ <= *igh) {
00103 goto L140;
00104 }
00105 if (i__ < *low) {
00106 i__ = *low - ii;
00107 }
00108 k = (integer) scale[i__];
00109 if (k == i__) {
00110 goto L140;
00111 }
00112
00113 i__2 = *m;
00114 for (j = 1; j <= i__2; ++j) {
00115 s = zr[i__ + j * zr_dim1];
00116 zr[i__ + j * zr_dim1] = zr[k + j * zr_dim1];
00117 zr[k + j * zr_dim1] = s;
00118 s = zi[i__ + j * zi_dim1];
00119 zi[i__ + j * zi_dim1] = zi[k + j * zi_dim1];
00120 zi[k + j * zi_dim1] = s;
00121
00122 }
00123
00124 L140:
00125 ;
00126 }
00127
00128 L200:
00129 return 0;
00130 }
00131