Re: [deal.II] dealii/spack on apple m1

2021-04-16 Thread Praveen C
Thanks for pointing this out. I think these recent fixes should be the solution to compiling on M1. I will give it a try. best praveen > On 17-Apr-2021, at 12:26 AM, Jean-Paul Pelteret wrote: > > Dear Praveen, > > Even though you’re trying to install deal.II, it would seem that the issue >

Re: [deal.II] LineMinimization and p-Laplace

2021-04-16 Thread Jean-Paul Pelteret
Hi Julie, So as you know from the error message, the assertion on this line is being triggered. The fit function is called from this section of code

Re: [deal.II] Reading a list from parameter file

2021-04-16 Thread Wolfgang Bangerth
Paras, I am trying to read a list of integers from the parameter file using the dealii::Patterns::List class. However, I get a runtime error which can be reproduced using the attached MWE. I also tried the approach described (for tensor) in

Re: [deal.II] interpolate FE_Nelelec

2021-04-16 Thread Jean-Paul Pelteret
To add to Wolfgang’s comment, you can have a look at the compute_edge_projection_l2() function that forms a part of the function that computes

Re: [deal.II] Re: interpolate FE_Nelelec

2021-04-16 Thread Wolfgang Bangerth
On 4/16/21 11:36 AM, John Smith wrote: However, the FE_Nedelec is different. It implements edge elements. They are vector-based. That is, functions are represented by a superposition of vector-valued shape functions: \vec{A} = \sum u_i vec{N}_i . Therefore, the output of the "value" method

Re: [deal.II] interpolate FE_Nelelec

2021-04-16 Thread John Smith
Dear Jean-Paul, Thank you for your reply. I am solving a very standard curl(1/mu(curl(A)) = J problem as described in the paper of Oszkar Biro. https://ieeexplore.ieee.org/document/497322 No gauge. I need interpolation to project the vector current potential on the shape functions. I just

[deal.II] Re: interpolate FE_Nelelec

2021-04-16 Thread John Smith
Ah! Now I got it. The function you suggested indeed vector-valued. It works just fine now. Thanks a lot! I was trying to output the function values via “vector_value” method like this: class CustomFunction : public Function<3> { public: virtual void vector_value( const Point<3> & p,

Re: [deal.II] dealii/spack on apple m1

2021-04-16 Thread Jean-Paul Pelteret
Dear Praveen, Even though you’re trying to install deal.II, it would seem that the issue that you documented sysctl: unknown oid 'machdep.cpu.leaf7_features' sysctl: unknown oid 'machdep.cpu.vendor' sysctl: unknown oid 'machdep.cpu.model' sysctl: unknown oid 'machdep.cpu.leaf7_features' sysctl:

Re: [deal.II] interpolate FE_Nelelec

2021-04-16 Thread Jean-Paul Pelteret
Dear John, I’m not sure that there is an interpolation function that would work in that way for Nedelec elements (there is VectorTools::project_boundary_values_curl_conforming_l2() for

Re: [deal.II] Re: interpolate FE_Nelelec

2021-04-16 Thread John Smith
Dear Simon, Thank you for your reply. Your suggestion definitely works. I used functions that output one component when I was interpolating scalar potentials. It works well. However, the FE_Nedelec is different. It implements edge elements. They are vector-based. That is, functions are

[deal.II] Re: interpolate FE_Nelelec

2021-04-16 Thread simon...@gmail.com
Hi, I don't have any experience with FE_Nedelec, but the error message states that the function you are trying to interpolate has 1 component, while the element has dim components. If you implemented your own Function, you need to make sure it has dim-components by calling the constructor of

[deal.II] interpolate FE_Nelelec

2021-04-16 Thread John Smith
Hello, It seems I am unable to find a function similar to VectorTools::interpolate for FE_Nedelec finite elements. The existing implementation of this functions gives the following run-time error: *An error occurred in line <556> of file in function* * void