Re: [deal.II] Catching exception for find_active_cell_around_point?

2016-08-17 Thread Alex Zimmerman
I'm new to deal.II and I'm having a similar problem. I'm trying to extend the FEFieldFunction class with an extrapolated_value method that handles this exception. Here's the relevant method: template double ExtrapolatedField:: extrapolated_value(const Point &point,

Re: [deal.II] Catching exception for find_active_cell_around_point?

2016-08-17 Thread Alex Zimmerman
I left out something important: I'm on a 8.5.0-pre version. (I don't see a way to edit my post.) On Wednesday, August 17, 2016 at 12:18:30 PM UTC+2, Alex Zimmerman wrote: > > I'm new to deal.II and I'm having a similar problem. I'm trying to extend > the FEFieldFunction class with an extrapolat

Re: [deal.II] Catching exception for find_active_cell_around_point?

2016-08-17 Thread Alex Zimmerman
I just stepped through my program with qtcreator and I found a mistake that undermines my question. I thought I was using this extrapolated_value method, but I was still using the parent class's value method. I imagine that this will work once I fix that mistake. I'll update this thread if I ma

Re: [deal.II] Catching exception for find_active_cell_around_point?

2016-08-17 Thread Alex Zimmerman
That was indeed my problem. Silly mistake. This is why we use IDE's before posting to the mailing list. Also an FYI for anyone who stumbles upon this conversation with myself: I extended Function instead, adding a FEFieldFunction as a member. On Wednesday, August 17, 2016 at 2:08:57 PM UTC+2,

Re: [deal.II] weakly enforcing an identity - vector mean curvature as the surface laplacian of the identity on a codim-1 manifold

2016-08-17 Thread Wolfgang Bangerth
On 08/17/2016 10:22 AM, thomas stephens wrote: *On to the question: *It is a fact that the vector mean curvature k_bar is equal to the surface laplacian of the identity function id_X on the manifold X, k_bar = laplace_beltrami id_X. Equation (4.4) in [1] is the weak form of this identity: (k_b

[deal.II] Re: [dealii-developers] dealii.org website down?

2016-08-17 Thread Wolfgang Bangerth
On 08/17/2016 11:30 AM, thomas stephens wrote: Currently unable to load https://www.dealii.org/ along with tutorials (problem began at approximately 12:30pm ETD) Yes, there is a power outage on campus in Heidelberg, Germany, where the website is hosted. Guido says that people are working on it

[deal.II] deal.II now also in Colorado

2016-08-17 Thread Wolfgang Bangerth
All, it is my pleasure to announce that since yesterday, my share of the deal.II development will continue from Colorado State University. Future user and developer workshops will therefore no longer have to happen in the sweltering heat of Texas, but may instead be in the beautiful state of Col

[deal.II] Re: Data exchanges in IO

2016-08-17 Thread bjpalmer58
On Tuesday, August 16, 2016 at 2:47:18 AM UTC-7, Daniel Arndt wrote: > > Bruce, > > [...] 3. Like in 2., but for adding values to individual elements. Use >> VectorOperation::add. >> >> I guess 3) is obscure. Does this mean that you are adding a new element >> where none existed before or that

[deal.II] Re: Data exchanges in IO

2016-08-17 Thread Daniel Arndt
Bruce, You were just missing p_rho.compress(VectorOperation::insert); after assemble_system(). Best, Daniel Am Mittwoch, 17. August 2016 21:01:22 UTC+2 schrieb bjpal...@gmail.com: > > > > On Tuesday, August 16, 2016 at 2:47:18 AM UTC-7, Daniel Arndt wrote: >> >> Bruce, >> >> [...] 3. Like in 2.