Doxygen Source Code Documentation
Main Page Alphabetical List Data Structures File List Data Fields Globals Search
eis_figi2.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006 #include "f2c.h"
00007
00008 int figi2_(integer *nm, integer *n, doublereal *t,
00009 doublereal *d__, doublereal *e, doublereal *z__, integer *ierr)
00010 {
00011
00012 integer t_dim1, t_offset, z_dim1, z_offset, i__1, i__2;
00013
00014
00015 double sqrt(doublereal);
00016
00017
00018 static doublereal h__;
00019 static integer i__, j;
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
00068
00069
00070
00071 t_dim1 = *nm;
00072 t_offset = t_dim1 + 1;
00073 t -= t_offset;
00074 z_dim1 = *nm;
00075 z_offset = z_dim1 + 1;
00076 z__ -= z_offset;
00077 --e;
00078 --d__;
00079
00080
00081 *ierr = 0;
00082
00083 i__1 = *n;
00084 for (i__ = 1; i__ <= i__1; ++i__) {
00085
00086 i__2 = *n;
00087 for (j = 1; j <= i__2; ++j) {
00088
00089 z__[i__ + j * z_dim1] = 0.;
00090 }
00091
00092 if (i__ == 1) {
00093 goto L70;
00094 }
00095 h__ = t[i__ + t_dim1] * t[i__ - 1 + t_dim1 * 3];
00096 if (h__ < 0.) {
00097 goto L900;
00098 } else if (h__ == 0) {
00099 goto L60;
00100 } else {
00101 goto L80;
00102 }
00103 L60:
00104 if (t[i__ + t_dim1] != 0. || t[i__ - 1 + t_dim1 * 3] != 0.) {
00105 goto L1000;
00106 }
00107 e[i__] = 0.;
00108 L70:
00109 z__[i__ + i__ * z_dim1] = 1.;
00110 goto L90;
00111 L80:
00112 e[i__] = sqrt(h__);
00113 z__[i__ + i__ * z_dim1] = z__[i__ - 1 + (i__ - 1) * z_dim1] * e[i__] /
00114 t[i__ - 1 + t_dim1 * 3];
00115 L90:
00116 d__[i__] = t[i__ + (t_dim1 << 1)];
00117
00118 }
00119
00120 goto L1001;
00121
00122
00123 L900:
00124 *ierr = *n + i__;
00125 goto L1001;
00126
00127
00128 L1000:
00129 *ierr = (*n << 1) + i__;
00130 L1001:
00131 return 0;
00132 }
00133