<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">On 21/06/17 23:12, Matthew Knepley
wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4GmYPx5+qoEQhM9N2axoZ9yWrw+j2fjx9Gj1knE-wRe9ZQ@mail.gmail.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote"><br>
<div>From the above description, its not clear to me that
you want this topology. For example, in the case that each
cell gets an internal cell,</div>
<div>it seems like you could just handle this in the
function space for each cell. Even for multiple internal
cells, it seems like function space</div>
<div>parameterization is a better option. Topology is
supposed to indicate the support of basis functions, but
here that job is done. I would</div>
<div>just treat it as some sort of augmented approximation
space.</div>
</div>
</div>
</div>
</blockquote>
<br>
If I understand what you mean, I considered doing something like
that- basically just defining extra degrees of freedom in the cells
where dual porosity is to be applied.<br>
<br>
It seemed to me that if I then went ahead and created the Jacobian
matrix using DMCreateMatrix(), it would give me extra nonzero
entries that shouldn't be there - interactions between the dual
porosity variables in neighbouring cells. Is there any way to avoid
that?<br>
<br>
- Adrian<br>
<pre class="moz-signature" cols="72">--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: <a class="moz-txt-link-abbreviated" href="mailto:a.croucher@auckland.ac.nz">a.croucher@auckland.ac.nz</a>
tel: +64 (0)9 923 4611
</pre>
</body>
</html>