Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Getting LinAlg error that should not happen #3

Open
daob opened this issue May 5, 2009 · 0 comments
Open

Getting LinAlg error that should not happen #3

daob opened this issue May 5, 2009 · 0 comments

Comments

@daob
Copy link
Owner

daob commented May 5, 2009

Searching /home/daob/work/automtmm/tests/ess-round3/SK/IMSMETN...
Found r3imsmetnSK.LS8.
Running LISREL for r3imsmetnSK.LS8,Command is "wine Lisrel85.exe /home/daob/work/automtmm/tests/ess-round3/SK/IMSMETN/r3imsmetnSK.LS8 OUT".
LISREL terminated normally.
Calculating standardized matrices for r3imsmetnSK.LS8...
Reading matrix BE...BE is a FUll matrix.
Reading matrix AL...AL is a VEctor.
Reading matrix GA...GA is a FUll matrix.
Reading matrix PH...PH is a SYmmetric matrix.
Reading matrix PS...PS is a SYmmetric matrix.
Reading matrix KA...KA is a VEctor.
Reading matrix TX...TX is a VEctor.
Reading matrix TY...TY is a VEctor.
Reading matrix TD...TD is a DIagonal matrix.
Reading matrix TE...TE is a SYmmetric matrix.
Reading matrix LX...LX is a FUll matrix.
Reading matrix LY...LY is a FUll matrix.

Getting matrices for group 1 (igrp=0)

LinAlgError Traceback (most recent call last)

/home/daob/work/automtmm/ in ()

/home/daob/work/automtmm/walk_and_run.py in walk_and_run(top_dir, tempdir, action)
158 break
159 if solution_obtained(os.path.join(tempdir, 'OUT')):
--> 160 smats = lisfile.standardize_matrices()
161 for igrp in range(len(smats)):
162 for matname, stanmat in smats[igrp].iteritems():

/home/daob/work/automtmm/parse_lisrel.pyc in standardize_matrices(self, path)
262 C = np.matrix(np.sqrt(np.diag(np.diag(Eee))))
263
--> 264 ga_s = C.I * ga * np.sqrt(np.diag(np.diag(ph)))
265 ly_s = np.sqrt(np.diag(1/np.diag(Eyy))) * ly * C
266

/usr/lib/python2.6/dist-packages/numpy/core/defmatrix.pyc in getI(self)
513 else:
514 from numpy.dual import pinv as func
--> 515 return asmatrix(func(self))
516
517 def getA(self):

/usr/lib/python2.6/dist-packages/numpy/linalg/linalg.pyc in inv(a)
348 """
349 a, wrap = _makearray(a)
--> 350 return wrap(solve(a, identity(a.shape[0], dtype=a.dtype)))
351
352

/usr/lib/python2.6/dist-packages/numpy/linalg/linalg.pyc in solve(a, b)
247 results = lapack_routine(n_eq, n_rhs, a, n_eq, pivots, b, n_eq, 0)
248 if results['info'] > 0:
--> 249 raise LinAlgError, 'Singular matrix'
250 if one_eq:
251 return wrap(b.ravel().astype(result_t))

LinAlgError: Singular matrix

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant