{"action":"closed","issue":{"url":"https://api.github.com/repos/pmgbergen/porepy/issues/75","repository_url":"https://api.github.com/repos/pmgbergen/porepy","labels_url":"https://api.github.com/repos/pmgbergen/porepy/issues/75/labels{/name}","comments_url":"https://api.github.com/repos/pmgbergen/porepy/issues/75/comments","events_url":"https://api.github.com/repos/pmgbergen/porepy/issues/75/events","html_url":"https://github.com/pmgbergen/porepy/issues/75","id":274812616,"number":75,"title":"error in extrusion + non_conforming","user":{"login":"alessiofumagalli","id":3007391,"avatar_url":"https://avatars2.githubusercontent.com/u/3007391?v=4","gravatar_id":"","url":"https://api.github.com/users/alessiofumagalli","html_url":"https://github.com/alessiofumagalli","followers_url":"https://api.github.com/users/alessiofumagalli/followers","following_url":"https://api.github.com/users/alessiofumagalli/following{/other_user}","gists_url":"https://api.github.com/users/alessiofumagalli/gists{/gist_id}","starred_url":"https://api.github.com/users/alessiofumagalli/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/alessiofumagalli/subscriptions","organizations_url":"https://api.github.com/users/alessiofumagalli/orgs","repos_url":"https://api.github.com/users/alessiofumagalli/repos","events_url":"https://api.github.com/users/alessiofumagalli/events{/privacy}","received_events_url":"https://api.github.com/users/alessiofumagalli/received_events","type":"User","site_admin":false},"labels":[],"state":"closed","locked":false,"assignee":{"login":"keileg","id":1837598,"avatar_url":"https://avatars1.githubusercontent.com/u/1837598?v=4","gravatar_id":"","url":"https://api.github.com/users/keileg","html_url":"https://github.com/keileg","followers_url":"https://api.github.com/users/keileg/followers","following_url":"https://api.github.com/users/keileg/following{/other_user}","gists_url":"https://api.github.com/users/keileg/gists{/gist_id}","starred_url":"https://api.github.com/users/keileg/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/keileg/subscriptions","organizations_url":"https://api.github.com/users/keileg/orgs","repos_url":"https://api.github.com/users/keileg/repos","events_url":"https://api.github.com/users/keileg/events{/privacy}","received_events_url":"https://api.github.com/users/keileg/received_events","type":"User","site_admin":false},"assignees":[{"login":"keileg","id":1837598,"avatar_url":"https://avatars1.githubusercontent.com/u/1837598?v=4","gravatar_id":"","url":"https://api.github.com/users/keileg","html_url":"https://github.com/keileg","followers_url":"https://api.github.com/users/keileg/followers","following_url":"https://api.github.com/users/keileg/following{/other_user}","gists_url":"https://api.github.com/users/keileg/gists{/gist_id}","starred_url":"https://api.github.com/users/keileg/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/keileg/subscriptions","organizations_url":"https://api.github.com/users/keileg/orgs","repos_url":"https://api.github.com/users/keileg/repos","events_url":"https://api.github.com/users/keileg/events{/privacy}","received_events_url":"https://api.github.com/users/keileg/received_events","type":"User","site_admin":false}],"milestone":null,"comments":6,"created_at":"2017-11-17T10:05:32Z","updated_at":"2017-12-13T10:37:29Z","closed_at":"2017-12-13T10:37:29Z","author_association":"COLLABORATOR","body":"Hi, I have a problem when I couple the extrusion and the non-conforming algorithm. However, I think the main problem is due to the latter. Note: if considering the extrusion couple with a conforming algorithm it works.\r\nMinimal example:\r\n\r\n```python\r\nimport numpy as np                                                               \r\nimport scipy.sparse as sps                                                       \r\n                                                                                 \r\nfrom porepy.viz.exporter import Exporter                                         \r\nfrom porepy.fracs import importer                                                \r\n                                                                                 \r\nfrom porepy.fracs import extrusion                                               \r\nfrom porepy.fracs.fractures import Fracture, FractureNetwork                     \r\nfrom porepy.fracs import meshing                                                 \r\n                                                                                                                                                               \r\nfolder_export = 'example_4_vem/'                                                 \r\nfile_export = 'vem'                                                              \r\ntol = 1e-8                                                                       \r\n                                                                                 \r\n# The fractures are not aligned with the mesh                                    \r\npts = np.array([[0, 0], [2, 2], [1, 1], [1, 3]]).T                               \r\nedges = np.array([[0, 1], [2, 3]]).T                                             \r\n                                                                                 \r\nextrustion_kwargs = {}                                                           \r\nextrustion_kwargs['tol'] = tol                                                   \r\nextrustion_kwargs['exposure_angle'] = 2*np.pi/3.*np.ones(edges.shape[1])         \r\n                                                                                 \r\nfractures = extrusion.fractures_from_outcrop(pts, edges, **extrustion_kwargs)    \r\nnetwork = FractureNetwork(fractures, tol=tol)                                    \r\n                                                                                 \r\n# If uncomment the code work with the new network                                \r\n# f_1 = Fracture(np.array([[-1, -1, 0], [1, -1, 0], [1, 1, 0], [-1, 1, 0]]).T)   \r\n# f_2 = Fracture(np.array([[-1, 0, -1], [1, 0, -1], [1, 0, 1], [-1, 0, 1]]).T)   \r\n# f_3 = Fracture(np.array([[-0.75, 0.25, -0.75], [0.75, 0.25, -0.75], [0.75,     \r\n#                            0.25, 0.75], [-0.75, 0.25, 0.75]]).T)               \r\n# network = FractureNetwork([f_1, f_2, f_3])                                     \r\n                                                                                 \r\nmesh_kwargs = {}                                                                 \r\nmesh_kwargs['mesh_size'] = {'mode': 'weighted',                                  \r\n                            'value': 0.5,                                        \r\n                            'bound_value': 0.5,                                  \r\n                            'tol': tol}                                          \r\n                                                                                 \r\ngb = meshing.dfn(network, conforming=False, **mesh_kwargs)                       \r\n                                                                                 \r\ngb.remove_nodes(lambda g: g.dim == 0)                                            \r\ngb.compute_geometry()                                                            \r\ngb.assign_node_ordering()                                                        \r\n                                                                                 \r\nsave = Exporter(gb, file_export, folder_export)                                  \r\nsave.write_vtk()\r\n```"}}
{"actor":{"display_login":"keileg"}}