Re: [deal.II] Re: Trouble getting UMFPACK to work

[陈敏]

You can install lapack/blas of netlib version or openblas version or intel
mkl version

best
chen

On Wed, Jul 6, 2022 at 1:05 PM ME20D503 NEWTON  wrote:

> Dear Daniel,
> i am getting this while running cmake-DDEAL_II_WITH_UMFPACK=ON, can you
> please see  the highlighted output?
>
> 
>
>  No source or binary directory provided.  Both will be assumed to be the
>   same as the current working directory, but note that this warning will
>   become a fatal error in future CMake releases.
>
>
> -- This is CMake 3.16.3
> --
> -- Include /home/newton/fem/dealii-9.2.0/cmake/setup_external_macros.cmake
> -- Include /home/newton/fem/dealii-9.2.0/cmake/macros/macro_add_flags.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_compiler_setup.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_cxx_compiler_bug.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_mpi_interface.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_clear_cmake_required.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_configure_feature.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_definitions.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_library.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_test.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_file.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_library.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_path.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_program.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_initialize_cached_variables.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_insource_setup_target.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_invoke_autopilot.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_package_handle.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_pickup_tests.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_query_git_information.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_setup_target.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_decorate_with_stars.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_if_links.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_if_supported.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_language_optional.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_evaluate_expression.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_expand_instantiations.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_filter_system_libraries.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_find_package.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_find_system_library.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_item_matches.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_purge_feature.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_register_feature.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_remove_duplicates.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_replace_flag.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_reset_cmake_required.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_set_if_empty.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_setup_source_list.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_setup_unity_target.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_strip_flag.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_switch_library_preference.cmake
> -- Include /home/newton/fem/dealii-9.2.0/cmake/macros/macro_to_string.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_to_string_and_add_prefix.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_unset_if_changed.cmake
> -- Include
> /home/newton/fem/dealii-9.2.0/cmake/macros/macro_verbose_include.cmake
> --
> -- Include /home/newton/fem/dealii-9.2.0/cmake/setup_cached_variables.cmake
> -- Prepending ${CUDA_NVCC_FLAGS} to ${DEAL_II_CUDA_FLAGS}
> -- The CXX compiler identification is GNU 9.4.0
> -- The C compiler identification is GNU 9.4.0
> -- Check for working CXX compiler: /usr/bin/c++
> -- Check for working CXX compiler: /usr/bin/c++ -- works
> -- Detecting CXX compiler ABI info
> -- Detecting CXX compiler 

Re: [deal.II] Re: Trouble getting UMFPACK to work

[ME20D503 NEWTON]

Dear Daniel,
i am getting this while running cmake-DDEAL_II_WITH_UMFPACK=ON, can you 
please see  the highlighted output?



 No source or binary directory provided.  Both will be assumed to be the
  same as the current working directory, but note that this warning will
  become a fatal error in future CMake releases.


-- This is CMake 3.16.3
-- 
-- Include /home/newton/fem/dealii-9.2.0/cmake/setup_external_macros.cmake
-- Include /home/newton/fem/dealii-9.2.0/cmake/macros/macro_add_flags.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_compiler_setup.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_cxx_compiler_bug.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_check_mpi_interface.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_clear_cmake_required.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_configure_feature.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_definitions.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_library.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_add_test.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_file.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_library.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_path.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_find_program.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_initialize_cached_variables.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_insource_setup_target.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_invoke_autopilot.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_package_handle.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_pickup_tests.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_query_git_information.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_deal_ii_setup_target.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_decorate_with_stars.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_if_links.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_if_supported.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_enable_language_optional.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_evaluate_expression.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_expand_instantiations.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_filter_system_libraries.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_find_package.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_find_system_library.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_item_matches.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_purge_feature.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_register_feature.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_remove_duplicates.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_replace_flag.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_reset_cmake_required.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_set_if_empty.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_setup_source_list.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_setup_unity_target.cmake
-- Include /home/newton/fem/dealii-9.2.0/cmake/macros/macro_strip_flag.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_switch_library_preference.cmake
-- Include /home/newton/fem/dealii-9.2.0/cmake/macros/macro_to_string.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_to_string_and_add_prefix.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_unset_if_changed.cmake
-- Include 
/home/newton/fem/dealii-9.2.0/cmake/macros/macro_verbose_include.cmake
-- 
-- Include /home/newton/fem/dealii-9.2.0/cmake/setup_cached_variables.cmake
-- Prepending ${CUDA_NVCC_FLAGS} to ${DEAL_II_CUDA_FLAGS}
-- The CXX compiler identification is GNU 9.4.0
-- The C compiler identification is GNU 9.4.0
-- Check for working CXX compiler: /usr/bin/c++
-- Check for working CXX compiler: /usr/bin/c++ -- works
-- Detecting CXX compiler ABI info
-- Detecting CXX compiler ABI info - done
-- Detecting CXX compile features
-- Detecting CXX compile features - done
-- Check for working C compiler: /usr/bin/cc
-- Check for working C compiler: /usr/bin/cc -- works
-- Detecting C compiler ABI info
-- Detecting C compiler ABI info - done
-- Detecting C compile features
-- Detecting C compile 

Re: [deal.II] UMFPACK ISSUE.

[陈敏]

After run cmake, there is a log to tell you whether the UMFPACK have
configure with the dealii, you should give the log to us. The cmake log
file can be found in build/cmakefile/**.log as well.

Best
Chen

ME20D503 NEWTON  于2022年7月6日周三 12:38写道：

> Dear Luca,
>
> can help me out, what this error is trying to tell me?
>
> To call this function you need UMFPACK, but you configured deal.II
> without passing the necessary switch to 'cmake'. Please consult the
> installation instructions in doc/readme.html.
>
>
> On Tuesday, July 5, 2022 at 12:41:03 PM UTC+5:30 luca@gmail.com wrote:
>
>> This part of the message contains all the information you need to fix the
>> issue.
>>
>> Luca
>>
>> > Il giorno 5 lug 2022, alle ore 06:46, ME20D503 NEWTON <
>> me20...@iittp.ac.in> ha scritto:
>> >
>> > Additional information:
>> > To call this function you need UMFPACK, but you configured deal.II
>> > without passing the necessary switch to 'cmake'. Please consult the
>> > installation instructions in doc/readme.html.
>>
> --
> 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 Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/3f7a8855-f572-4145-8454-0a8dc97bb532n%40googlegroups.com
> 
> .
>

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CADr3OdKVk-QjkCSr5wVz_3%3DiF2jg9cH6Ww7J7E7VvTXkV7--dg%40mail.gmail.com.

Re: [deal.II] UMFPACK ISSUE.

[ME20D503 NEWTON]

Dear Luca,

can help me out, what this error is trying to tell me?

To call this function you need UMFPACK, but you configured deal.II
without passing the necessary switch to 'cmake'. Please consult the
installation instructions in doc/readme.html.


On Tuesday, July 5, 2022 at 12:41:03 PM UTC+5:30 luca@gmail.com wrote:

> This part of the message contains all the information you need to fix the 
> issue. 
>
> Luca
>
> > Il giorno 5 lug 2022, alle ore 06:46, ME20D503 NEWTON <
> me20...@iittp.ac.in> ha scritto:
> > 
> > Additional information: 
> > To call this function you need UMFPACK, but you configured deal.II
> > without passing the necessary switch to 'cmake'. Please consult the
> > installation instructions in doc/readme.html.
>

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/3f7a8855-f572-4145-8454-0a8dc97bb532n%40googlegroups.com.

Re: [deal.II] discretization of a analytical function on a second triangulation

[Daniel Arndt]

Simon,

defining a class derived from dealii::Function sounds sensible.
Whether calling VectorTools::interpolate makes sense, depends on what you
want to do with the interpolated function. In may cases, it's enough to
just evaluate it locally during assembly using functionality like
FEValues::get_function_values.

Best,
Daniel

On Mon, Jul 4, 2022 at 10:26 AM Simon  wrote:

> Dear all:
>
>
> I am solving a nonlinear PDE on a Triangulation T_1.
>
> I know the analytical representation of a scalar function of two variables
> (no space co-ordinates, but two invariants of a quantity)
> and my goal is to find a discretized version of the analytical function
> (on a second, two-dimensional triangulation T_2).
> The nodal values of T_2 should be the exact values corresponding to the
> analytical solution.
> As for the values between the nodes, I want to start with a bilinear
> interpolation.
>
> The coupling to my original PDE works as follows:
> For each quadrature point defined on T_1, I compute the two input
> variables of the above function based on my solution vector, find the
> corresponding cell of T_2 and want to retrieve the interpolated nodal
> values.
>
> My idea is to use the function
> VectorTools::interpolate ,
> but I do not know a suitable Function object to hand over.
>
> Is this an appropriate approach at all?
>
>
> Thank you,
> Simon
>
> --
> 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 Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/14b1d8f6-f902-4414-a5c6-348ae3f5003en%40googlegroups.com
> 
> .
>

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAOYDWbKNTV%3D8Bipyq_gzbMoLvxdvv-%2BmHLipn27DbOcK-eWL_A%40mail.gmail.com.

Re: [deal.II] Re: Trouble getting UMFPACK to work

[Daniel Arndt]

What does the log for running
cmake -DDEAL_II_WITH_UMFPACK=ON .
look like? The summary printed at the end should tell you if `UMFPACK`
really was found. N=Any problems in finding a compatible library should be
contained in the output printed when running the CMake command.

Best,
Daniel

On Mon, Jul 4, 2022 at 12:25 PM ME20D503 NEWTON 
wrote:

> Hello everyone,
> which command i have to use to check whether UMFPACK is enabled or not??
> if it is OFF how to turn it ON??
>
> Thank you in advance.
>
> On Thursday, August 31, 2017 at 9:28:52 PM UTC+5:30 John wrote:
>
>> Hello everyone,
>>
>> I am currently using Deal II version 8.0.0. When I try to run a project I
>> am getting the following:
>>
>> 
>> Exception on processing:
>>
>> 
>> An error occurred in line <379> of file 
>> 
>> in function
>> void dealii::SparseDirectUMFPACK::factorize(const Matrix&) [with
>> Matrix = dealii::SparseMatrix]
>> The violated condition was:
>> false
>> The name and call sequence of the exception was:
>> ExcMessage("To call this function you need UMFPACK, but configured
>> deal.II without passing the necessary switch to 'cmake'. Please consult the
>> installation instructions in doc/readme.html.")
>> Additional Information:
>> To call this function you need UMFPACK, but configured deal.II without
>> passing the necessary switch to 'cmake'. Please consult the installation
>> instructions in doc/readme.html.
>>
>> Stacktrace:
>> ---
>> #0  /Downloads/deal.II/lib/libdeal_II.so.8.0.0: void
>> dealii::SparseDirectUMFPACK::factorize
>> >(dealii::SparseMatrix const&)
>> #1  step-3: TmsProblem::solve()
>> #2  step-3: TmsProblem::run()
>> #3  step-3: main
>> 
>>
>> Aborting!
>> 
>> CMake Error at CMakeFiles/run_target.cmake:6 (MESSAGE):
>>
>>
>>   Program terminated with exit code: 1
>>
>>
>> CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run'
>> failed
>> make[3]: *** [CMakeFiles/run] Error 1
>> CMakeFiles/Makefile2:264: recipe for target 'CMakeFiles/run.dir/all'
>> failed
>> make[2]: *** [CMakeFiles/run.dir/all] Error 2
>> CMakeFiles/Makefile2:271: recipe for target 'CMakeFiles/run.dir/rule'
>> failed
>> make[1]: *** [CMakeFiles/run.dir/rule] Error 2
>> Makefile:196: recipe for target 'run' failed
>> make: *** [run] Error 2
>>
>>
>>
>> I have tried "*cmake -DDEAL_II_WITH_UMFPACK=ON*" and "*cmake
>> -DDEAL_II_WITH_UMFPACK=ON -DUMFPACK_DIR=/Downloads/deal.II/bundled/umfpack
>> .*" and also  "*cmake -DDEAL_II_WITH_UMFPACK=ON
>> -DUMFPACK_DIR=/Downloads/deal.II/bundled/umfpack/UMFPACK .*" and then
>> following each of these with "*make*" and "*make run*". However, all of
>> these still end up giving me the same error as above. Can anyone help me
>> with getting umfpack to work? Thank you in advance.
>>
> --
> 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 Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/f8629e26-974f-4dd9-aa0e-99a198de0f96n%40googlegroups.com
> 
> .
>

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAOYDWbLFsJycNqXUO44z52WHtGCWC42YCgjicXNdsui3A7wWVg%40mail.gmail.com.

Re: [deal.II] shape function gradient at arbitrary points in the element

[Simon Wiesheier]

Dear Martin,

thanks for pointing out to the FEPointEvaluation class.
However, my intent is to compute the gradients of shape functions in real
coordinates.

" so in that case the recommendation would be to use FEValues, despite that
being terribly inefficient."

What is the approach to achieve this?
My only idea is to create a new FEValues object for each quadrature point
and pass an appropriate Quadrature object. But seems to be very inefficient.
Is there an alternative way in combination with FEValues?

Best,
Simon


Am Di., 5. Juli 2022 um 08:16 Uhr schrieb Martin Kronbichler <
kronbichler.mar...@gmail.com>:

> Dear Simon,
>
> We have the class FEPointEvalation,
> https://dealii.org/developer/doxygen/deal.II/classFEPointEvaluation.html
> , which implements the operation of evaluating a solution at arbitrary
> points (as handed in as an array to points in unit coordinates). To do this
> for positions in real coordinates, you typically need to invert the
> mapping, Mapping::transform_real_to_unit_cell or
> Mapping::transform_points_real_to_unit_cell, see
> https://dealii.org/developer/doxygen/deal.II/classMapping.html , always
> assuming that you have located the correct points. Most tasks can be
> reduced to this setup.
>
> If your intent is not just the solution interpolation but the value of
> shape functions, e.g. because you want to assemble a local matrix,
> FEPointEvaluation is not directly providing that info, so in that case the
> recommendation would be to use FEValues, despite that being terribly
> inefficient. We could possibly integrate such functionality into
> FEPointEvaluation if you have interest in that direction, so let us know
> what your needs are.
>
> Best,
> Martin
> On 05.07.22 15:56, Simon wrote:
>
> Dear all:
>
> I have to compute the gradient with respect to real space co-ordinates of
> the shape functions belonging to a FE_Q element at arbitrary points in the
> element, that is, not only at quadrature points.
>
> What is the way to do this in dealii?
> FEValues provides shape function values, gradients,... only at quadrature
> points.
>
> Thank you,
> Simon
> --
> 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 Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/16afc8fa-bcfe-4a81-9ade-040c18f85596n%40googlegroups.com
> 
> .
>
> --
> 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 Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/34843cd0-22a6-d4b5-7cb8-f530715a6108%40gmail.com
> 
> .
>

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAM50jEteHa8vhryPshNi4SGP3jhm%3DsEj_A8K6wq1TSX-aamKmQ%40mail.gmail.com.

Re: [deal.II] shape function gradient at arbitrary points in the element

[Martin Kronbichler]

Dear Simon,

We have the class FEPointEvalation, 
https://dealii.org/developer/doxygen/deal.II/classFEPointEvaluation.html 
, which implements the operation of evaluating a solution at arbitrary 
points (as handed in as an array to points in unit coordinates). To do 
this for positions in real coordinates, you typically need to invert the 
mapping, Mapping::transform_real_to_unit_cell or 
Mapping::transform_points_real_to_unit_cell, see 
https://dealii.org/developer/doxygen/deal.II/classMapping.html , always 
assuming that you have located the correct points. Most tasks can be 
reduced to this setup.


If your intent is not just the solution interpolation but the value of 
shape functions, e.g. because you want to assemble a local matrix, 
FEPointEvaluation is not directly providing that info, so in that case 
the recommendation would be to use FEValues, despite that being terribly 
inefficient. We could possibly integrate such functionality into 
FEPointEvaluation if you have interest in that direction, so let us know 
what your needs are.


Best,
Martin

On 05.07.22 15:56, Simon wrote:

Dear all:

I have to compute the gradient with respect to real space co-ordinates 
of the shape functions belonging to a FE_Q element at arbitrary points 
in the element, that is, not only at quadrature points.


What is the way to do this in dealii?
FEValues provides shape function values, gradients,... only at 
quadrature points.


Thank you,
Simon
--
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send 
an email to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/16afc8fa-bcfe-4a81-9ade-040c18f85596n%40googlegroups.com 
.


--
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 Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/34843cd0-22a6-d4b5-7cb8-f530715a6108%40gmail.com.

[deal.II] shape function gradient at arbitrary points in the element

[Simon]

Dear all:

I have to compute the gradient with respect to real space co-ordinates of 
the shape functions belonging to a FE_Q element at arbitrary points in the 
element, that is, not only at quadrature points. 

What is the way to do this in dealii?
FEValues provides shape function values, gradients,... only at quadrature 
points. 

Thank you,
Simon

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/16afc8fa-bcfe-4a81-9ade-040c18f85596n%40googlegroups.com.

[deal.II] PhD project opportunity at Heidelberg University

[Michał Wichrowski]

We (Guido Kanschat and I) offer a PhD project devising a code for 
fluid-structure interaction for cell movement, which should allow for 
experiments with various fluid and solid mechanics models in cooperation 
with experimentalists.

We are looking for somebody with experience with deal.II and preferably 
fluid or solid mechanics.

A master's degree is required when the person starts. Heidelberg University 
is an equal opportunity employer where women are particularly invited to 
apply.

If you are interested, don't hesitate to get in touch with me at 
wichrow...@uni-heidelberg.de

Best,
Michał Wichrowski

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/2db1906d-a0a0-4899-95db-faef144468f8n%40googlegroups.com.

[deal.II] Re: DoFTools::make_periodicity_constraints error when working with FEColleciton

[jose.a...@gmail.com]

Hi Peter,

Thanks! I will install your dealii fork and test it out.

Cheers,
Jose

On Monday, July 4, 2022 at 7:25:51 PM UTC+2 peterr...@gmail.com wrote:

> Hi Rose,
>
> the assert was definitely not positioned at the right place and checked on 
> every face - even on non-active ones. I have moved it in the PR 
> https://github.com/dealii/dealii/pull/14091. Now this assert is not 
> triggered.
>
> Generally, I am wondering if hp and PBC is working. At least, I cannot 
> find a test in the hp test folder. Such a test should have triggered the 
> assert.
>
> Hope the patch works!
>
> Peter
>
> On Monday, 4 July 2022 at 18:09:41 UTC+2 jose.a...@gmail.com wrote:
>
>> Dear dealii community,
>>
>> I am still trying to decipher what causes the error. It occurs in the 
>> first assert of make_periodicity_constraints 
>> 
>> :
>>
>> AssertDimension(
>> face_1->get_fe(face_1->nth_active_fe_index(0)).n_unique_faces(), 1);
>>
>> So I created vectors to store the active_fe_index and the 
>> nth_active_fe_index(0), replicated the fe_index_is_active assert and the 
>> get_fe() call:
>>
>> active_fe_index.reinit(triangulation.n_active_cells());
>> nth_active_fe_index.reinit(triangulation.n_active_cells());
>>
>> active_fe_index = -1.0;
>> nth_active_fe_index = -1.0;
>>
>> for (const auto  :
>> dof_handler.active_cell_iterators())
>> if (cell->is_locally_owned())
>> {
>> active_fe_index(cell->active_cell_index()) = cell->active_fe_index();
>>
>> if (cell->at_boundary())
>> for (const auto _id : cell->face_indices())
>> if (cell->face(face_id)->at_boundary())
>> {
>> AssertThrow(cell->face(face_id)->get_active_fe_indices().size() == 1, 
>> dealii::ExcMessage("Error!"));
>>
>> nth_active_fe_index(cell->active_cell_index()) = cell->face(face_id)->
>> nth_active_fe_index(0);
>>
>> AssertThrow(cell->face(face_id)->fe_index_is_active(cell->face(face_id)->
>> nth_active_fe_index(0)), dealii::ExcMessage("Error!"));
>>
>> cell->face(face_id)->get_fe(cell->face(face_id)->nth_active_fe_index(0));
>> }
>> }
>>
>> I don't get any error from this loop and the from the visual inspection 
>> of .vtu file containing both vectors one can see that they indeed 
>> coincide. Maybe the problem is coming from GridTools::
>> collect_periodic_faces as the face iterators are created by it?
>>
>> Furthermore I also noticed that if no global refinement is performed the 
>> executable runs without error in serial but still gets the error in 
>> parallel.
>>
>> Attached is the updated MWE. A true/false argument can be passed to the 
>> executable to perform one global refinement or not.
>>
>> Cheers,
>> Jose
>>
>> On Wednesday, June 29, 2022 at 4:45:04 PM UTC+2 jose.a...@gmail.com 
>> wrote:
>>
>>> Dear dealii community,
>>>
>>> I am working on a gradient enhanced crystal plasticity code where the 
>>> displacement field is continuous but the slip fields can be discontinuous. 
>>> To this end, I am using a hp::FECollection. For example, in the 
>>> case of two crystals with one slip system I have a hp::FECollection<2> 
>>> with
>>>
>>> fe_collection[0] = [FE_Q<2> FE_Q<2> FE_Q<2>  FE_Nothing<2>]
>>> fe_collection[1] = [FE_Q<2> FE_Q<2> FE_Nothing<2> FE_Q<2> ]
>>>
>>> where the first two FE_Q<2> correspond to the displacement field in 2D 
>>> and the latter to the slip field.
>>>
>>> My reference problem consist of the unit square under simple shear with 
>>> periodic boundary conditions on the x-direction. When I call the 
>>> DoFTools::make_periodicity_constraints>> dim> method, I am getting the following error when I have subdomains 
>>> with different values of fe_index
>>>
>>> An error occurred in line <1690> of file 
>>> 
>>>  
>>> in function
>>> const dealii::FiniteElement& 
>>> dealii::DoFAccessor::get_fe(unsigned int) 
>>> const [with int structdim = 1; int dim = 2; int spacedim = 2; bool 
>>> level_dof_access = false]
>>> The violated condition was: 
>>> fe_index_is_active(fe_index) == true
>>> Additional information: 
>>> This function can only be called for active FE indices
>>>
>>> I wrote a minimum working example based on my code to showcase the 
>>> error, which I have attached to this post. The executable takes a true or 
>>> false argument to assign different material identifiers to different 
>>> subdomains or not. If there are no subdomains and therefore, and more 
>>> importantly, the size of the hp::FECollection is 1, the code runs 
>>> with no problem.
>>>
>>> My question is if my algorithm is wrong or does the method 
>>> DoFTools::make_periodicity_constraints>> dim> does not currently support hp::FECollection with sizes bigger 
>>> than one.
>>>
>>> Cheers,
>>> Jose
>>>
>>

-- 
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 Group" group.
To 

Re: [deal.II] UMFPACK ISSUE.

[Luca Heltai]

This part of the message contains all the information you need to fix the 
issue. 

Luca

> Il giorno 5 lug 2022, alle ore 06:46, ME20D503 NEWTON  
> ha scritto:
> 
> Additional information: 
> To call this function you need UMFPACK, but you configured deal.II
> without passing the necessary switch to 'cmake'. Please consult the
> installation instructions in doc/readme.html.

-- 
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 Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/568CF12F-1424-4A00-9F63-8FE1ACBD2D4A%40gmail.com.