René's URL Explorer Experiment

{"@context":"https://schema.org","@type":"DiscussionForumPosting","headline":"Installing Slycot (openmp) breaks numpy","articleBody":"After installing slycot in a conda environment, a kernel crash occurs when trying to solve a linear system with a singular matrix.\r\n\r\nI narrowed the issue down to the following:\r\n- slycot from conda-forge depends on libflang \u003e=5.0.0,\u003c6.0.0.a0\r\n- libflang 5.0.0 depends on openmp 5.0.0\r\n- openmp 5.0.0 overwrites \"libiomp5md.dll\" that was installed by numpy. The two dlls are very different in size (2mb for the numpy one, 600kb for the openmp one)\r\n\r\nOpenmp places its \"libiomp5md.dll\" in Library\\bin.\r\nThe version from numpy is either in Library\\bin (channel main) or directly in bin (conda-forge).\r\n\r\n### Actual Behavior\r\n\r\nUsing numpy.linalg.solve on singular matrix with dimension over 149 crashes the Kernel.\r\n\r\n### Expected Behavior\r\n\r\nExpected a \"LinAlgError: Singular matrix\" as it is shown with matrices up to the dimension of 149. \r\n\r\n### Steps to Reproduce\r\n\r\n- Create a new conda environment with python 3.11, numpy (and ipython for convenience)\r\n  - `mamba create -n slycot_crash python=3.11 numpy ipython -c conda-forge`\r\n- activate the environment\r\n  - `mamba activate slycot_crash`\r\n- enter ipython\r\n  - `ipython`\r\n- paste and run the code below to confirm it works as expected (raises LinAlgError)\r\n- exit ipython\r\n  -  `exit`\r\n- install openmp=5\r\n  - `mamba install openmp=5 -c conda-forge`\r\n- enter ipython and run the code below again to confirm kernel crash\r\n- (optional) change n=149 and it only fails if the size is 150 or larger\r\n\r\n```\r\nimport numpy as np\r\n\r\nn = 150\r\nA = np.random.rand(n,n)\r\nA[-1,:] = np.zeros((1,n)) # make matrix singular\r\n\r\nb = np.random.rand(n)\r\n\r\nx = np.linalg.solve(A, b)\r\n```","author":{"url":"https://github.com/chnce","@type":"Person","name":"chnce"},"datePublished":"2024-01-26T09:32:06.000Z","interactionStatistic":{"@type":"InteractionCounter","interactionType":"https://schema.org/CommentAction","userInteractionCount":12},"url":"https://github.com/230/Slycot/issues/230"}