{"action":"created","issue":{"url":"https://api.github.com/repos/chan-y-park/loom/issues/2","labels_url":"https://api.github.com/repos/chan-y-park/loom/issues/2/labels{/name}","comments_url":"https://api.github.com/repos/chan-y-park/loom/issues/2/comments","events_url":"https://api.github.com/repos/chan-y-park/loom/issues/2/events","html_url":"https://github.com/chan-y-park/loom/issues/2","id":100135941,"number":2,"title":"get_ramification_points() misses some ramification points","user":{"login":"neitzke","id":8607462,"avatar_url":"https://avatars.githubusercontent.com/u/8607462?v=3","gravatar_id":"","url":"https://api.github.com/users/neitzke","html_url":"https://github.com/neitzke","followers_url":"https://api.github.com/users/neitzke/followers","following_url":"https://api.github.com/users/neitzke/following{/other_user}","gists_url":"https://api.github.com/users/neitzke/gists{/gist_id}","starred_url":"https://api.github.com/users/neitzke/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/neitzke/subscriptions","organizations_url":"https://api.github.com/users/neitzke/orgs","repos_url":"https://api.github.com/users/neitzke/repos","events_url":"https://api.github.com/users/neitzke/events{/privacy}","received_events_url":"https://api.github.com/users/neitzke/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":21,"created_at":"2015-08-10T18:26:49Z","updated_at":"2015-09-02T17:40:50Z","closed_at":null,"body":"If I use the Seiberg-Witten curve\r\n```\r\nx**4 - 10*z*x**2 + 4*x + 9*z**2 = 0\r\n```\r\nand all other parameters as in default.ini, I expect to get 6 branch points (as in Figure 16 of 1209.0866) but I only get 4. It seems that some of them are being discarded by the check \"if mx > 1\" in get_ramification_points()."},"comment":{"url":"https://api.github.com/repos/chan-y-park/loom/issues/comments/137183386","html_url":"https://github.com/chan-y-park/loom/issues/2#issuecomment-137183386","issue_url":"https://api.github.com/repos/chan-y-park/loom/issues/2","id":137183386,"user":{"login":"neitzke","id":8607462,"avatar_url":"https://avatars.githubusercontent.com/u/8607462?v=3","gravatar_id":"","url":"https://api.github.com/users/neitzke","html_url":"https://github.com/neitzke","followers_url":"https://api.github.com/users/neitzke/followers","following_url":"https://api.github.com/users/neitzke/following{/other_user}","gists_url":"https://api.github.com/users/neitzke/gists{/gist_id}","starred_url":"https://api.github.com/users/neitzke/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/neitzke/subscriptions","organizations_url":"https://api.github.com/users/neitzke/orgs","repos_url":"https://api.github.com/users/neitzke/repos","events_url":"https://api.github.com/users/neitzke/events{/privacy}","received_events_url":"https://api.github.com/users/neitzke/received_events","type":"User","site_admin":false},"created_at":"2015-09-02T17:40:50Z","updated_at":"2015-09-02T17:40:50Z","body":"For polynomials without repeated zeroes, the Python command ```mpmath.polyroots``` seems to do the job, e.g.\r\n```\r\nmpmath.mp.dps = 100\r\nmpmath.polyroots([589824,0,0,-143360,0,0,-6912])\r\n```\r\nfinds zeroes of the discriminant in the example above.\r\n\r\nMaybe this method is OK for us actually: at least in many examples, the discriminant will be a rational function, and we could just use the above to find roots of the numerator. It saves the trouble of shelling out to sage/maxima/whatever and seems to work to arbitrary precision."}}