{"action":"created","issue":{"url":"https://api.github.com/repos/libMesh/libmesh/issues/830","repository_url":"https://api.github.com/repos/libMesh/libmesh","labels_url":"https://api.github.com/repos/libMesh/libmesh/issues/830/labels{/name}","comments_url":"https://api.github.com/repos/libMesh/libmesh/issues/830/comments","events_url":"https://api.github.com/repos/libMesh/libmesh/issues/830/events","html_url":"https://github.com/libMesh/libmesh/issues/830","id":134598109,"number":830,"title":"Infinite elments revitalisation","user":{"login":"BalticPinguin","id":2708704,"avatar_url":"https://avatars.githubusercontent.com/u/2708704?v=3","gravatar_id":"","url":"https://api.github.com/users/BalticPinguin","html_url":"https://github.com/BalticPinguin","followers_url":"https://api.github.com/users/BalticPinguin/followers","following_url":"https://api.github.com/users/BalticPinguin/following{/other_user}","gists_url":"https://api.github.com/users/BalticPinguin/gists{/gist_id}","starred_url":"https://api.github.com/users/BalticPinguin/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/BalticPinguin/subscriptions","organizations_url":"https://api.github.com/users/BalticPinguin/orgs","repos_url":"https://api.github.com/users/BalticPinguin/repos","events_url":"https://api.github.com/users/BalticPinguin/events{/privacy}","received_events_url":"https://api.github.com/users/BalticPinguin/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":1,"created_at":"2016-02-18T14:42:01Z","updated_at":"2016-02-18T14:50:57Z","closed_at":null,"body":"The first part I already mentioned some time ago via email: If you run the miscellaneous ex1 in debug-mode, you get the following message:\r\n``` \r\n Building Infinite Elements:\r\n  updating element neighbor tables...\r\n  collecting boundary sides...\r\n  found 96 inner and 0 outer boundary faces\r\n  added 0 infinite elements and 0 nodes to the mesh \r\n```\r\nSo in fact, there are no infinite elements. To fix this, I would suggest the following patch:\r\n[infinite.patch.txt](https://github.com/libMesh/libmesh/files/136317/infinite.patch.txt)\r\n\r\nWith this, two problems occur: First, exodus refuses to work here; Roy Stogner explained me that \" we use one block per subdomain id, and Exodus doesn't allow more than one element type\r\nper block.\" The easiest way probably would be to remove the statement or to use an other format at this point; maybe in a longer term one could change the format such that exodus is possible here as well!?\r\n\r\nHowever, even after this the program fails, with the error code\r\n```/usr/include/c++/4.8/debug/vector:353:error: attempt to subscript container \r\n    with out-of-bounds index 0, but container only holds 0 elements.\r\n\r\nObjects involved in the operation:\r\nsequence \"this\" @ 0x0x1475d30 {\r\n  type = NSt7__debug6vectorINS0_IdSaIdEEESaIS2_EEE;\r\n}\r\nAborted (core dumped)```\r\nwhich is because the elements Ss and St in combine_base_radial() (src/fe/inf_fe.C) are empty:\r\n`(gdb) p Ss.size()\r\n$1 = 0\r\n(gdb) p St.size()\r\n$2 = 0\r\n(gdb) p mode.size()\r\n$3 = 4`\r\nAs a fix, I would suggest the following:\r\n[infinite2.patch.txt](https://github.com/libMesh/libmesh/files/136339/infinite2.patch.txt)\r\n\r\nwhich made the example working for me finally.\r\nBest,\r\nHubert\r\n\r\n"},"comment":{"url":"https://api.github.com/repos/libMesh/libmesh/issues/comments/185753124","html_url":"https://github.com/libMesh/libmesh/issues/830#issuecomment-185753124","issue_url":"https://api.github.com/repos/libMesh/libmesh/issues/830","id":185753124,"user":{"login":"permcody","id":3036960,"avatar_url":"https://avatars.githubusercontent.com/u/3036960?v=3","gravatar_id":"","url":"https://api.github.com/users/permcody","html_url":"https://github.com/permcody","followers_url":"https://api.github.com/users/permcody/followers","following_url":"https://api.github.com/users/permcody/following{/other_user}","gists_url":"https://api.github.com/users/permcody/gists{/gist_id}","starred_url":"https://api.github.com/users/permcody/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/permcody/subscriptions","organizations_url":"https://api.github.com/users/permcody/orgs","repos_url":"https://api.github.com/users/permcody/repos","events_url":"https://api.github.com/users/permcody/events{/privacy}","received_events_url":"https://api.github.com/users/permcody/received_events","type":"User","site_admin":false},"created_at":"2016-02-18T14:50:56Z","updated_at":"2016-02-18T14:50:56Z","body":"On Thu, Feb 18, 2016 at 7:42 AM Hubert Weißmann <notifications@github.com>\nwrote:\n\n> The first part I already mentioned some time ago via email: If you run the\n> miscellaneous ex1 in debug-mode, you get the following message:\n>\n>  Building Infinite Elements:\n>   updating element neighbor tables...\n>   collecting boundary sides...\n>   found 96 inner and 0 outer boundary faces\n>   added 0 infinite elements and 0 nodes to the mesh\n>\n> So in fact, there are no infinite elements. To fix this, I would suggest\n> the following patch:\n> infinite.patch.txt\n> <https://github.com/libMesh/libmesh/files/136317/infinite.patch.txt>\n>\n> With this, two problems occur: First, exodus refuses to work here; Roy\n> Stogner explained me that \" we use one block per subdomain id, and Exodus\n> doesn't allow more than one element type\n> per block.\" The easiest way probably would be to remove the statement or\n> to use an other format at this point; maybe in a longer term one could\n> change the format such that exodus is possible here as well!?\n>\nHow about simply creating another subdomain block if you have another\nelement type? That's the normal usage case for handling exodus format. It's\nnot an arbitrary restriction. It's so that exodus can reduce the\ninformation it packs for each element. If you know the type for the whole\nblock, then you can simply list the connectivity id list for each element\nand no other information is required.\n\n\n> However, even after this the program fails, with the error code\n> ```/usr/include/c++/4.8/debug/vector:353:error: attempt to subscript\n> container\n> with out-of-bounds index 0, but container only holds 0 elements.\n>\n> Objects involved in the operation:\n> sequence \"this\" @ 0x0x1475d30 {\n> type = NSt7__debug6vectorINS0_IdSaIdEEESaIS2_EEE;\n> }\n> Aborted (core dumped)``\n> which is because the elements Ss and St in combine_base_radial()\n> (src/fe/inf_fe.C) are empty:\n> (gdb) p Ss.size()\n> $1 = 0\n> (gdb) p St.size()\n> $2 = 0\n> (gdb) p mode.size()\n> $3 = 4`\n> As a fix, I would suggest the following:\n> infinite2.patch.txt\n> <https://github.com/libMesh/libmesh/files/136339/infinite2.patch.txt>\n>\n> which made the example working for me finally.\n> Best,\n> Hubert\n>\n> —\n> Reply to this email directly or view it on GitHub\n> <https://github.com/libMesh/libmesh/issues/830>.\n>\n"}}