[deal.II] Extracting Shape Laplacian

2017-08-09 Thread kylew
Hello all, I am working on implementing some stabilization methods for the Navier-Stokes equations, these techniques often require the use of the laplacian of the shape functions. I know I get the hessian of the k'th shape function at quadratue point q via the following Tensor<3, dim> hessian

[deal.II] Re: Extracting Shape Laplacian

2017-08-09 Thread kylew
I found that my idea above was not what I had typed in my code, I am no longer receiving the run time error. However I would still like some verification that the indices are as I think they are. > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, se

Re: [deal.II] Extracting Shape Laplacian

2017-08-09 Thread kylew
Thanks, I did try looking at the documentation for FEValuesViews::Vector < dim, spacedim >::hessian(const unsigned int shape_function, const unsigned int q_point) before coming here, and all it provides is "Return

[deal.II] Interpolated Boundary Conditions from Distributed Solution

2018-11-29 Thread kylew
Hello All, I am trying to interpolate a finite element solution to use as dirichlet boundary conditions for an auxiliary problem. I have the following working for serial runs: Functions::FEFieldFunction, TrilinosWrappers::MPI:: BlockVector> boundaryconditions(dof_handler, curl_vel); VectorT

Re: [deal.II] Interpolated Boundary Conditions from Distributed Solution

2018-12-13 Thread kylew
Thank you very much! I will report back with a minimal test case. I am still curious if this would generally be the recommended approach for what I'm trying to accomplish? On Saturday, December 1, 2018 at 5:05:27 AM UTC-6, Jean-Paul Pelteret wrote: > > Dear Kyle, > > I’m not too familiar with F

Re: [deal.II] Interpolated Boundary Conditions from Distributed Solution

2018-12-13 Thread kylew
Thank you very much! I will report back with a minimal test case. I am still curious though, would this generally be the correct way to implement what I am trying to achieve? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.

[deal.II] Trouble Installing: "undefined symbol: cblas_ctrmv"

2019-01-24 Thread kylew
from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout. ### # # deal.II configuration: #CMAKE_BUILD_TYPE: DebugRelease #BUILD_SHARED_LIBS: ON #CMAKE_INST

[deal.II] Re: Trouble Installing: "undefined symbol: cblas_ctrmv"

2019-01-24 Thread kylew
Thank you very much, this will work fine. -Kyle williams -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User G

[deal.II] Trouble Running Exercises

2019-01-24 Thread kylew
with deal.II-9.1.0-pre found at # /home/kylew/local/deal.II I'm then able to compile just fine, but when I run the resulting executable i get the following error ./step-1: error while loading shared libraries: libdeal_II.so.9.1.0-pre: cannot open shared object file: No such file or dire

[deal.II] Re: Trouble Installing: "undefined symbol: cblas_ctrmv"

2019-01-25 Thread kylew
Hi Denis, I did see both of the issues raised on github before I posted, but it wasn't really clear what the resolution was other than its a gsl dependency problem. nm -g /usr/lib64/libgslcblas.so nm -g /usr/lib64/libgsl.so both tell me there are no symbols, even without grepping for the one i

[deal.II] Re: Trouble Installing: "undefined symbol: cblas_ctrmv"

2019-01-25 Thread kylew
I see now that the symbols are actually there, but they are dynamic symbols which don't seem to get picked up (Though I don't really know what the difference is). -Kyle -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/

[deal.II] Mistake in Step 16?

2019-02-05 Thread kylew
Hi All, I was working through example step 16 and I think I found a mistake. In LaplaceProblem::setup_system() We have the following std::set dirichlet_boundary_ids; Functions::ZeroFunction