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

Reply via email to