RE: [deal.II] Questions on Step-18

2021-07-07 Thread Michael Li
Hi Jean-Paul, Thank you for your patience and detailed explanation! It did help me understand the formulations. The only thing after reading the Wikipedia (https://en.wikipedia.org/wiki/Rotation_matrix#Axis_and_angle ) I still have doubts is that the angle should be calculated by ½ mag(curl) and

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Jean-Paul Pelteret
Hi Matthew, Great, I’m glad to hear that! > So my last question is what is the difference between a flag being ON vs set > up with external dependencies ( which is what trilinos is set at now) When you set one of these flags to “ON”, you are demanding from CMake that this option be enabled

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Matthew Rich
Thanks that worked! I needed to install the full trilinos-all-dev to get all the packages. I also needed to install libptscotch-dev. So my last question is what is the difference between a flag being ON vs set up with external dependencies ( which is what trilinos is set at now) Will I be

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Jean-Paul Pelteret
Hi again Matthew, I’ve also found this build log for the Debian deal.II apt package https://buildd.debian.org/status/fetch.php?pkg=deal.ii=all=9.2.0-3=1606990185=0 and you can see that in amongst all of the

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Jean-Paul Pelteret
Hi Matthew, If you look at either log file, then you’ll see that Trilinos was in fact not detected for some reason. > DEAL_II_WITH_TRILINOS = OFF I’m no longer experienced in building deal.II using apt packages, but it might be that you have to install the entire trilinos-dev to get all of

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Matthew Rich
So I tried running cmake again to get trilinos picked up. It does not seem to be detected. I even try supplying the directory with cmake -/usr/local ./ -DTRILINOS_DIR=/usr/include/trilinos But it still does not get detected. Should I build trilinos from source as well? I attached another

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Matthew Rich
Hi Thanks for the prompt reply. Yes I built deal.II from source but libtrilinos-sacado-dev was installed via the debian package and apt (log files attached). I have all the headers in my /usr/include/trilinos directory. Just skimming the summary file it appears that deal.II did not detect

Re: [deal.II] New location for sacado.hpp

2021-07-07 Thread Jean-Paul Pelteret
Hi Matthew, I don’t quite understand exactly how you installed deal.II and Trilinos. Perhaps you can provide some more explicit details so that we can help you out. * Is Trilinos installed using the apt package, and deal.II installed from source? If you built deal.II from source, may you

[deal.II] New location for sacado.hpp

2021-07-07 Thread Matthew Rich
Hi all, I am a bit confused on step-33 which uses the differentiation of sacado. I installed the latest version. 9.3 which has trilinios in it, and therefore sacado. Browsing that directory I see other header files for the package but not sacado.hpp The tutorial has the line #include Was this

Re: [deal.II] triangulation.signals.post_refinement :connection of a member function

2021-07-07 Thread Simon Wiesheier
I will look that up :-) Just one last question: In the constructed a class I store a DoFHandler and Triangulation object as member variables like this: SmartPointer> dof_handler; SmartPointer> triangulation; I created the instance at the same time when constructing the Solid instance (pM =

[deal.II] Re: Solving an interface problem using matrix-free

2021-07-07 Thread 'David Schneider' via deal.II User Group
I can answer it myself after a bit of work, maybe there is someone wondering about similar things in the future: > (1) will the SolverControl emit convergence although I operate only on some of the DoFs (do I need to reset the remaining ones) Yes it will. Just reset the irrelevant DoFs and

[deal.II] Unexpected data output with cell data vector

2021-07-07 Thread vachanpo...@gmail.com
Dear all, I am having some confusion regarding the index used in DataOut for a cell vector. I have created a distributed vector for storing cell data using p::d::Triangulation::global_active_cell_index_partitioner()