Doxygen Source Code Documentation
Main Page Alphabetical List Data Structures File List Data Fields Globals Search
eis_elmbak.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006 #include "f2c.h"
00007
00008 int elmbak_(integer *nm, integer *low, integer *igh,
00009 doublereal *a, integer *int__, integer *m, doublereal *z__)
00010 {
00011
00012 integer a_dim1, a_offset, z_dim1, z_offset, i__1, i__2, i__3;
00013
00014
00015 static integer i__, j;
00016 static doublereal x;
00017 static integer la, mm, mp, kp1, mp1;
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
00066
00067 --int__;
00068 a_dim1 = *nm;
00069 a_offset = a_dim1 + 1;
00070 a -= a_offset;
00071 z_dim1 = *nm;
00072 z_offset = z_dim1 + 1;
00073 z__ -= z_offset;
00074
00075
00076 if (*m == 0) {
00077 goto L200;
00078 }
00079 la = *igh - 1;
00080 kp1 = *low + 1;
00081 if (la < kp1) {
00082 goto L200;
00083 }
00084
00085 i__1 = la;
00086 for (mm = kp1; mm <= i__1; ++mm) {
00087 mp = *low + *igh - mm;
00088 mp1 = mp + 1;
00089
00090 i__2 = *igh;
00091 for (i__ = mp1; i__ <= i__2; ++i__) {
00092 x = a[i__ + (mp - 1) * a_dim1];
00093 if (x == 0.) {
00094 goto L110;
00095 }
00096
00097 i__3 = *m;
00098 for (j = 1; j <= i__3; ++j) {
00099
00100 z__[i__ + j * z_dim1] += x * z__[mp + j * z_dim1];
00101 }
00102
00103 L110:
00104 ;
00105 }
00106
00107 i__ = int__[mp];
00108 if (i__ == mp) {
00109 goto L140;
00110 }
00111
00112 i__2 = *m;
00113 for (j = 1; j <= i__2; ++j) {
00114 x = z__[i__ + j * z_dim1];
00115 z__[i__ + j * z_dim1] = z__[mp + j * z_dim1];
00116 z__[mp + j * z_dim1] = x;
00117
00118 }
00119
00120 L140:
00121 ;
00122 }
00123
00124 L200:
00125 return 0;
00126 }
00127