Re: [deal.II] How to do DoFRenumbering::component_wise for multilevel grid

2022-11-28 Thread 'yy.wayne' via deal.II User Group
I was debugging in debug mode. Just changed make to release mode and the 
error gone! Why running in debug mode
causes this error?

As for the side note, that's very helpful! I worried about the write 
mesh price but don't know what techniques can I use.
That's exactly what I want. Thank you.

Best,
Wayne

在2022年11月29日星期二 UTC+8 12:30:15 写道:

>
> > The error is from compute_component_wise on level = 0, the result it 
> returns 
> > don't equal to dof_handler.n_dofs(level).
> > I added a MPI_Barrier between first and second loop but error is not 
> eliminated.
> > Snipaste_2022-11-28_12-53-02.png
>
> I do not see this error. Are you running in debug mode? And with how many 
> MPI 
> processes?
>
> As a side note: The detour by writing the data to a file and then reading 
> it 
> again is very inefficient, in particular if you have many MPI processes. 
> If 
> you want to use GridOut/GridIn, let the former write the data into an 
> object 
> of type std::ostringstream, then use Utilities::MPI::broadcast() to send 
> the 
> string you get out of std::ostringstream to the other processes, and give 
> it 
> to a std::istringstream object. Then use GridIn on the std::istringstream 
> object. This way, no data transfer to or from disk is necessary -- it all 
> happens in memory.
>
> Best
> Wolfgang
>
> -- 
> 
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
>

-- 
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/9be828ee-3e6d-4cee-ae54-6fcf390752c6n%40googlegroups.com.


[deal.II] How can I integrate along a surface/edge within the domain? Can I tag the surface/edge with an ID?

2022-11-28 Thread Abbas
Hello, 

As per the screenshot below, I created and imported a mesh using gmsh, with 
two material IDs. I want to integrate on the edge that separate the two 
materials for my weak form. I tried tagging the internal edge with an ID in 
gmsh to loop on it, but of course dealii throws an error saying that I 
cannot tag an internal edge with an ID. 

How can I integrate on an arbitrary edge like this one? 

-- 
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/2f2f050b-ca85-43b4-a88b-9829070a846bn%40googlegroups.com.


Re: [deal.II] How to do DoFRenumbering::component_wise for multilevel grid

2022-11-28 Thread Wolfgang Bangerth



The error is from compute_component_wise on level = 0, the result it returns 
don't equal to dof_handler.n_dofs(level).

I added a MPI_Barrier between first and second loop but error is not eliminated.
Snipaste_2022-11-28_12-53-02.png


I do not see this error. Are you running in debug mode? And with how many MPI 
processes?


As a side note: The detour by writing the data to a file and then reading it 
again is very inefficient, in particular if you have many MPI processes. If 
you want to use GridOut/GridIn, let the former write the data into an object 
of type std::ostringstream, then use Utilities::MPI::broadcast() to send the 
string you get out of std::ostringstream to the other processes, and give it 
to a std::istringstream object. Then use GridIn on the std::istringstream 
object. This way, no data transfer to or from disk is necessary -- it all 
happens in memory.


Best
 Wolfgang

--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/


--
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/35e96ca4-0b33-1870-cde9-2e42c94eddf4%40colostate.edu.


Re: [deal.II] Re: How to solve 2d elasticity problem for two different materials.

2022-11-28 Thread Deepika Kushwah
Thank you for your response.


Thanks & Regards,
Deepika



On Mon, Nov 28, 2022 at 7:39 PM Bruno Turcksin 
wrote:

> Deepika,
>
> You can take a look at step-28. We use material_id to represent different
> materials.
>
> Best,
>
> Bruno
>
> On Monday, November 28, 2022 at 1:56:22 AM UTC-5 deepika...@iitgoa.ac.in
> wrote:
>
>> Hello Everyone,
>>
>> I have solved the 2d elasticity problem for DBCs and NBCs for a single
>> material. Now I wanted to solve it for different materials (Ex. steel & Al)
>> in the same domain. Can you please suggest which tutorial will help to do
>> the same?
>>
>>
>> Thanks & Regard,
>> Deepika
>>
> --
> 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/791842ec-2591-4db5-8b3c-9801e843ea91n%40googlegroups.com
> 
> .
>

-- 
**

This e-mail is for the sole use of the intended recipient(s) and may

contain confidential and privileged information. If you are not the

intended recipient, please contact the sender by reply e-mail and destroy

all copies and the original message. Any unauthorized review, use,

disclosure, dissemination, forwarding, printing or copying of this email


is strictly prohibited and appropriate legal action will be taken. 




-- 
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/CAKTdrprqe%3DZOjCC5JmBfFMcc4rOQSThhqyiP%2BcjP5CzKFWbErw%40mail.gmail.com.


[deal.II] deal.II Newsletter #234

2022-11-28 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone!

This is deal.II newsletter #234.
It automatically reports recently merged features and discussions about the 
deal.II finite element library.


## Below you find a list of recently proposed or merged features:

#14486: Add more reinit() to LA::distributed::BlockVector. (proposed by 
marcfehling; merged) https://github.com/dealii/dealii/pull/14486

#14485: Utilities::MPI::scatter: improve asserts (proposed by mschreter; 
merged) https://github.com/dealii/dealii/pull/14485

#14484: Some AffineConstraints instantiations. (proposed by marcfehling; 
merged) https://github.com/dealii/dealii/pull/14484

#14483: CMake: Update forgotten find_package() call (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14483

#14482: CMake: Rework process_feature() macro (proposed by tamiko) 
https://github.com/dealii/dealii/pull/14482

#14481: A few improvements to the PETSc classes (proposed by stefanozampini) 
https://github.com/dealii/dealii/pull/14481

#14480: CMake: Update CMake cache config documentation (proposed by tamiko) 
https://github.com/dealii/dealii/pull/14480

#14479: CMake: ensure some more, modern style rules (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14479

#14478: Remove obsolete doxygen tags. (proposed by marcfehling; merged) 
https://github.com/dealii/dealii/pull/14478

#14477: Inline most of the remaining functions in VtuStream into 
write_vtu_main(). (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/14477

#14476: Cleanup reinit() of PETSc/Trilinos BlockVector. (proposed by 
marcfehling) https://github.com/dealii/dealii/pull/14476

#14475: Introduce DEAL_II_HOST_DEV_ALWAYS_INLINE (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/14475

#14474: Update documentation for deal.II+Cuda (proposed by masterleinad; 
merged) https://github.com/dealii/dealii/pull/14474

#14473: add a changelog for SparsityPatternBase. (proposed by drwells; merged) 
https://github.com/dealii/dealii/pull/14473

#14472: CMake: rename all find modules to DEAL_II_* (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14472

#14471: CMake: Remove _USER variables (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/14471

#14470: [9.4] Take over #14462 and #14465 (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/14470

#14469: add a news/changes entry (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/14469

#14468: CMake: Bugfix: only configure arborx with matching MPI configuration 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14468

#14467: CMake: Introduce a dealii:: namespace for targets and harmonize naming 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14467

#14466: CMake: search quietly for kokkos and avoid an unnecessary warning 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14466

#14465: CMake: always export compile_commands.json in deal.II and user projects 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14465

#14464: Add Utilities::MPI::scatter (proposed by mschreter; merged) 
https://github.com/dealii/dealii/pull/14464

#14463: CMake: lowercase all function and macro names (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14463

#14462: gitignore: ignore clangd files and directories (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14462

#14461: Test Kokkos on GPU (proposed by zjiaqi2018) 
https://github.com/dealii/dealii/pull/14461

#14460: Fix more warnings about unused variables (proposed by kronbichler; 
merged) https://github.com/dealii/dealii/pull/14460

#14459: clang-tidy: use line breaks in config, remove nonexistent check 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14459

#14457: Add Kokkos link libraries to Trilinos link libraries (proposed by 
masterleinad) https://github.com/dealii/dealii/pull/14457

#14456: Fix PETScWrappers::PreconditionBase warnings (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/14456

#14455: CMake: remove forgotten GO = statement in tpetra compatibility check 
(proposed by tamiko; merged) https://github.com/dealii/dealii/pull/14455

#14453: Replace DEAL_II_CUDA_HOST_DEV with KOKKOS_FUNCTION (proposed by Rombur; 
merged) https://github.com/dealii/dealii/pull/14453

#14452: GridGenerator: avoid a warning in release mode (proposed by tamiko; 
merged) https://github.com/dealii/dealii/pull/14452

#14451: Sparsity pattern base cleanup (proposed by drwells; merged) 
https://github.com/dealii/dealii/pull/14451

#14450: Bump actions/checkout from 2 to 3 (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/14450

#14449: Move declaration of vtu_stream to where it is used. (proposed by 
bangerth; merged) https://github.com/dealii/dealii/pull/14449

#14448: Use variable-sized std::array instead of fixed number of arguments. 
(proposed by bangerth; merged) 

[deal.II] Re: How to solve 2d elasticity problem for two different materials.

2022-11-28 Thread Bruno Turcksin
Deepika,

You can take a look at step-28. We use material_id to represent different 
materials.

Best,

Bruno

On Monday, November 28, 2022 at 1:56:22 AM UTC-5 deepika...@iitgoa.ac.in 
wrote:

> Hello Everyone,
>
> I have solved the 2d elasticity problem for DBCs and NBCs for a single 
> material. Now I wanted to solve it for different materials (Ex. steel & Al) 
> in the same domain. Can you please suggest which tutorial will help to do 
> the same?
>
>
> Thanks & Regard,
> Deepika
>

-- 
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/791842ec-2591-4db5-8b3c-9801e843ea91n%40googlegroups.com.


Re: [deal.II] Compilation error

2022-11-28 Thread Pedro Rodrigues
Hello
I think I have solved the problem. I had a custom lapack installation in
/usr/local/lib and Trilinos was picking that one instead of the version
from *candi* at /usr/lib. After removing the first one I was able to build
the library with *candi*. I noticed too many warnings related to casting
from *void*. I hope that it doesn't create execution problems.
Another problem I have  found was related to *mpiexec*. It sends a warning
about running as *root* and does not run. I had to add  two *export* lines
in *.bashrc* to bypass it.
My system is UbuntuLTS 22.04 running in VBox.

thanks for your quick reply

regards

On Sun, Nov 27, 2022 at 11:20 PM Wolfgang Bangerth 
wrote:

> On 11/26/22 07:33, Pedro Rodrigues wrote:
> >
> > I am trying to build the library with your candi script but I got an
> error,
> > apparently with a missing function from Trilinos library:
> >
> > [ 75%] Linking CXX executable ../bin/step-1.debug
> > /usr/bin/ld:
> /root/dealii-candi/trilinos-release-13-2-0/lib/libepetra.so:
> > undefined reference to `dggsvd_'
> > /usr/bin/ld:
> /root/dealii-candi/trilinos-release-13-2-0/lib/libepetra.so:
> > undefined reference to `sggsvd_'
> > collect2: error: ld returned 1 exit status
> > make[2]: *** [examples/CMakeFiles/step-1.debug.dir/build.make:249:
> > bin/step-1.debug] Error 1
> > make[1]: *** [CMakeFiles/Makefile2:3998:
> > examples/CMakeFiles/step-1.debug.dir/all] Error 2
> > make: *** [Makefile:136: all] Error 2
>
> Pedro,
> the missing functions are part of BLAS. It is a bug in Trilinos that it
> does
> not make sure the BLAS library appears on the linker line, but I don't
> know
> how to fix this. What kind of system are you on? Can you just install
> deal.II
> from an existing package?
>
> Best
>   W.
>
> --
> 
> Wolfgang Bangerth  email: bange...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
>
> --
> 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/6e6c9529-cec5-133a-7deb-14c8800324fd%40colostate.edu
> .
>


-- 
atenciosamente,

*Pedro Rodrigues*

-- 
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/CA%2BajKM_v2qzurfAz9N3imOr4iM5JVT-8gBM8G83K71q47ErJ-Q%40mail.gmail.com.