I am not entirely sure what that version is, but the assembly of point
integrals was not supported in 1.4.0 (see the ChangeLog). This means
that any assembly of point integrals would return zero quietly. Consider
upgrading to 1.5.
--
Marie
On 02/23/2015 01:43 PM, Wilhelm Braun wrote:
Hi Marie,
I am using the following DOLFIN version:
1.4.0+dfsg-2~ppa3~precise1
--Wilhelm
On 21/02/15 19:05, Marie E. Rognes wrote:
Hi Wilhelm,
Which FEniCS version are you using?
--
Marie
On 20 Feb 2015, at 18:22, Wilhelm Braun <[email protected]
<mailto:[email protected]>> wrote:
I tried to implement the goal function using the point integral as
outlined in the point integral demo in DOLFIN.
However, the point integral evaluated to zero, and hence, the
adaptive algorithm stopped after one iteration, but did not throw an
error.
Please find the code for the definition of the goal function M below.
--Wilhelm
----------------------------------------------------------------------------------------------------
code for goal function M
----------------------------------------------------------------------------------------------------
def center_func(x):
return (vbar-0.05 <= x[0] and x[0] <= vbar + 0.05 and
near(x[1], 0.0)) or -0.05 <= x[1] and x[1] <= 0.05 and near(x[0], vbar)
# Define domain for point integral
center_domain = VertexFunction("size_t", mesh, 0)
center = AutoSubDomain(center_func)
center.mark(center_domain, 1)
dPP = dP[center_domain]
M = u* dPP(1)
----------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------
On 20/02/15 11:23, Anders Logg wrote:
I suggest looking at the point integral demo in DOLFIN:
demo/undocumented/point-integral/python/demo_point-integral.py
Then use the point integral to define you goal functional. Warning:
I haven't tested this myself and don't know if it will break the
adaptive algorithm somehow. It would be interesting to get feedback
on whether it works (without any promise that I will fix it if it
does not... ;-)
--
Anders
Fri Feb 20 2015 at 11:47:10 AM skrev Wilhelm Braun
<[email protected] <mailto:[email protected]>>:
Dear all,
I have posted a question in the FEniCS QA forum, and it was
recommended that I bring up the issue on this mailing list.
Please have a look here for my question.
http://fenicsproject.org/qa/6555/automated-oriented-adaptivity-point-computational-domain
Thank you very much!
Wilhelm Braun
This message and any attachment are intended solely for the addressee
and may contain confidential information. If you have received this
message in error, please send it back to me, and immediately delete it.
Please do not use, copy or disclose the information contained in this
message or in any attachment. Any views or opinions expressed by the
author of this email do not necessarily reflect the views of the
University of Nottingham.
This message has been checked for viruses but the contents of an
attachment may still contain software viruses which could damage your
computer system, you are advised to perform your own checks. Email
communications with the University of Nottingham may be monitored as
permitted by UK legislation.
This message and any attachment are intended solely for the addressee
and may contain confidential information. If you have received this
message in error, please send it back to me, and immediately delete it.
Please do not use, copy or disclose the information contained in this
message or in any attachment. Any views or opinions expressed by the
author of this email do not necessarily reflect the views of the
University of Nottingham.
This message has been checked for viruses but the contents of an
attachment may still contain software viruses which could damage your
computer system, you are advised to perform your own checks. Email
communications with the University of Nottingham may be monitored as
permitted by UK legislation.
_______________________________________________
fenics mailing list
[email protected] <mailto:[email protected]>
http://fenicsproject.org/mailman/listinfo/fenics
This message and any attachment are intended solely for the addressee
and may contain confidential information. If you have received this
message in error, please send it back to me, and immediately delete it.
Please do not use, copy or disclose the information contained in this
message or in any attachment. Any views or opinions expressed by the
author of this email do not necessarily reflect the views of the
University of Nottingham.
This message has been checked for viruses but the contents of an
attachment may still contain software viruses which could damage your
computer system, you are advised to perform your own checks. Email
communications with the University of Nottingham may be monitored as
permitted by UK legislation.
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
