Re: [deal.II] More than one instance of _axpy

Hi Victor, I don't know exactly what is causing this but I would bet that it's related to the fix we applied in https://github.com/dealii/dealii/pull/10287 namely, some other package (I'm working from memory here, but I think it's Trilinos) defines its own LAPACK wrappers in a slightly differe

[deal.II] More than one instance of _axpy

``` [ 43%] Building CXX object CMakeFiles/aspect.dir/source/particle/interpolator/bilinear_least_squares.cc.o In file included from /opt/apps/intel18/cray_mpich_7_7/dealii/9.2.0/include/deal.II/lac/full_matrix.templates.h(33), from /work/00434/eijkhout/lonestar/aspect/aspect-2.2

Re: [deal.II] Unable to visualise high order polynomial variation in ParaView

One thing that might not be obvious is that it is not enough to have a manifold id for boundary faces. You need to assign a manifold to the cells as well. On Wednesday, January 20, 2021 at 11:49:59 AM UTC-5 luca@gmail.com wrote: > Did you call > > data_out.build_patches(mapping, > subdivi

Re: [deal.II] Unable to visualise high order polynomial variation in ParaView

Did you call data_out.build_patches(mapping, subdivisions,DataOut::curved_inner_cells); with subdivisions > 1? L. > On 20 Jan 2021, at 12:06, vachanpo...@gmail.com > wrote: > > Dear all, > > I have a working code that outputs pvtu/vtu files for visualisation. > > I have read the wiki page

[deal.II] Unable to visualise high order polynomial variation in ParaView

Dear all, I have a working code that outputs pvtu/vtu files for visualisation. I have read the wiki page for writing and visualising high order output. Currently, I only have straight-edge quadrilaterals so I the o