TRIQS / triqs_0.x

DEPRECATED -- This is the repository of the older versions of TRIQS
Other
11 stars 9 forks source link

SemiCircular inverse problem #71

Closed srbhp closed 12 years ago

srbhp commented 12 years ago

i was trying to use a script http://ipht.cea.fr/triqs/doc/user_manual/green/tutorial.html#a-slightly-more-complicated-example . Replace Wilson by SemiCircular and it is not working.
Giving a error --

Traceback (most recent call last): File "gr2.py", line 12, in g['s','s'] = inverse(SemiCircular(1.5) ) File "/home/saurabh/triqs/triqs_install/lib/python2.7/dist-packages/pytriqs/Base/GF_Local/_GFBloc_concept_impl.py", line 88, in setitem g <<= val File "/home/saurabh/triqs/triqs_install/lib/python2.7/dist-packages/pytriqs/Base/GF_Local/_GFBloc_base_data_tail.py", line 157, in ilshift self.copyFrom ( lazy_expressions.eval_lazy_expr(e_t, A2) ) File "/home/saurabh/triqs/triqs_install/lib/python2.7/dist-packages/pytriqs/Base/GF_Local/lazy_expressions.py", line 98, in eval_lazy_expr return f ( *map(lambda e :eval_lazy_expr(eval_term,e) , expr.childs[1:]) ) File "/home/saurabh/triqs/triqs_install/lib/python2.7/dist-packages/pytriqs/Base/GF_Local/inverse.py", line 30, in inverse cop.invert() File "/home/saurabh/triqs/triqs_install/lib/python2.7/dist-packages/pytriqs/Base/GF_Local/_GFBloc_base_data_tail.py", line 316, in invert d[:,:,om ] = numpy.linalg.inv(d[:,:,om]) File "/usr/lib/python2.7/dist-packages/numpy/linalg/linalg.py", line 445, in inv return wrap(solve(a, identity(a.shape[0], dtype=a.dtype))) File "/usr/lib/python2.7/dist-packages/numpy/linalg/linalg.py", line 328, in solve raise LinAlgError, 'Singular matrix' numpy.linalg.linalg.LinAlgError: Singular matrix

mferrero commented 12 years ago

Hi. Thanks for the input, we indeed didn't have a complete implementation of the SemiCircular on the real axis. I think I fixed it in 336d831. Let me know if it works!

srbhp commented 12 years ago

Thanks a lot. Yes, its working now.