[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

down the equivalence here: <https://lh3.googleusercontent.com/-BfOP_i7Vx_s/WOYKfP1BodI/UE0/0_3yExW3mGAvbRKhOElaF5rTccfwkHcwgCLcB/s1600/NonlinearNS_TensorContraction.PNG> On Thursday, April 6, 2017 at 10:41:07 AM UTC+2, Alex Zimmerman wrote: > > Correction! When I tested you

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

, but you saved me from using that ugly loop. Thanks! On Thursday, April 6, 2017 at 10:33:46 AM UTC+2, Alex Zimmerman wrote: > > Thanks for following up. > > I don't think that > double sum = v*gradz*w > > is equivalent to > double sum = 0.; > for (unsign

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

Thanks for following up. I don't think that double sum = v*gradz*w is equivalent to double sum = 0.; for (unsigned int i = 0; i < dim ; ++i) { for (unsigned int j = 0; j < dim; ++j) { sum += _w[j]*_gradz[i][j]*_v[i]; } } Most importantly, making the change definitely

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

solved fields if no extractor is used, so you'd have to > reconstruct "S" manually, in a similar way to that which is done in step-18. > > I hope that this helps clarify things a little bit. > > Best regards, > Jean-Paul > > > On Wednesday, April 5, 2017 at 1:10:0

[deal.II] Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

I am debugging a code where part of my weak form requires indexing into the Tensor<2,dim> that results from getting the gradient of a Tensor<1,dim> (in this case, a velocity vector). I found the convention for the index notation here:

Re: [deal.II] How can I resize a std::vector<Functions::ParsedFunction>?

and > assigning them. > > Think of this example: > std::vector<int*> v; > v.resize(4); > now v[0] is a pointer to an int, but it is NULL unless you do something > like > v[0] = new int; > > > > On Wed, Mar 22, 2017 at 3:09 PM, Alex Zimmerman > &l

[deal.II] How can I resize a std::vector<Functions::ParsedFunction>?

Fundamentally I am trying to allow for a variable number of ParsedFunction objects to be specified in a parameter input file. Maybe there is a better approach which circumvents my issue below. I can continue my work for some time with this being a constant; but as soon as I want to extend to 3D

Re: [deal.II] Why is this asserted to be impossible?

Bangerth wrote: > > On 09/17/2016 05:20 AM, Alex Zimmerman wrote: > > > > Not understanding why this is impossible, I began digging through the > code a > > bit and thought I would write my own 1D implementation. I noticed that > the > > code in the existing te

[deal.II] Why is this asserted to be impossible?

I have a code that is an extension of step-26 with non-homogeneous Neumann boundary conditions and zero forcing function, among other changes. Until recently I had only actually run the code in 2D, though I have been continuing the "dimension independent" programming. Yesterday I tried running

[deal.II] Re: Untraceable "not implemented" error in GDB

Yep that does it! Thanks for your patience :) On Friday, September 16, 2016 at 7:28:41 PM UTC+2, Bruno Turcksin wrote: > > Alex, > > before the run command, you need to use: catch throw > > Best, > > Bruno > > On Friday, September 16, 2016 at 1:21:41

Re: [deal.II] Re: Failing to set up simple test configuration

ments so that they can be considered in the > discussion > > about getting deal.II to properly on windows. > > > > Cheers, > > J-P > > > > > > On Wednesday, September 7, 2016 at 2:56:10 PM UTC+2, Alex Zimmerman > wrote: >

Re: [deal.II] Re: Failing to set up simple test configuration

on the Windows filesystem (from the bash prompt). On Wednesday, September 7, 2016 at 10:36:21 PM UTC+2, Wolfgang Bangerth wrote: > > On 09/07/2016 06:56 AM, Alex Zimmerman wrote: > > > > After some brief Googling, I think I might even be able to fix this > issue if I > >

[deal.II] Re: Failing to set up simple test configuration

ii/dealii/issues/3060> related to running deal.II > on windows. Its interesting to hear that you've got it running on Windows > 10. Out of curiosity, how have you configured deal.II? Do you link it > against any external libraries? > > On Wednesday, September 7, 2016 at 10:07:06 AM UTC

[deal.II] Re: Failing to set up simple test configuration

2016 at 10:07:06 AM UTC+2, Alex Zimmerman wrote: > > This might be a problem with my file system. I've been successfully using > the Ubuntu subsystem on Windows 10 for the past few weeks, and that's where > the odd "/mnt/c/Users..." path comes from in my question. > &

[deal.II] Re: Failing to set up simple test configuration

and in path? > > Regards, > J-P > > On Wednesday, September 7, 2016 at 9:57:33 AM UTC+2, Alex Zimmerman wrote: >> >> Maybe this is a bug, but I think it's more likely that I'm doing >> something wrong here. >> >> I had tried to set up a test for my use

[deal.II] Re: Failing to set up simple test configuration

ts/my_test.debug .. Passed5.51 sec > 100% tests passed, 0 tests failed out of 1 > Total Test time (real) = 5.52 sec One of these days I'll ask a question that I don't answer myself a few minutes later :) On Wednesday, September 7, 2016 at 9:57:33 AM UTC+2, Alex Zimmerman wrote: > >

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

, Alex Zimmerman wrote: > > 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 f

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

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