[deal.II] Re: which version of trilinos for dealii-8.4.1

[Roc Wang]

Thanks, everybody.

   I installed trilinos-12.0.1, with Muelu enabled. Dealii-8.4.1 can be 
installed successfully. 

   There was an error when I installed trilinos-12.6.3 which stuck me. The 
error info is like:

*Scanning dependencies of target muelu-adapters* 
[ 97%] Building CXX object packages/muelu/adapters/CMakeFiles/
muelu-adapters.dir/epetra/MueLu_EpetraOperator.cpp.o 
[ 97%] Building CXX object packages/muelu/adapters/CMakeFiles/
muelu-adapters.dir/stratimikos/Stratimikos_MueLuHelpers.cpp.o 
g++: internal compiler error: Killed (program cc1plus) 
Please submit a full bug report, 
with preprocessed source if appropriate. 
See ; for instructions. 
packages/muelu/adapters/CMakeFiles/muelu-adapters.dir/build.make:86: recipe 
for target 'packages/muelu/adapters/CMakeFiles/
muelu-adapters.dir/stratimikos/S
tratimikos_MueLuHelpers.cpp.o' failed 
make[2]: *** [packages/muelu/adapters/CMakeFiles/
muelu-adapters.dir/stratimikos/Stratimikos_MueLuHelpers.cpp.o] Error 4 
CMakeFiles/Makefile2:30692: recipe for target 
'packages/muelu/adapters/CMakeFiles/muelu-adapters.dir/all' failed 
make[1]: *** [packages/muelu/adapters/CMakeFiles/muelu-adapters.dir/all] 
Error 2 
Makefile:160: recipe for target 'all' failed 
make: *** [all] Error 2.

 After 12.0.1 installed, dealii can be compiled and built withour errors. 
Thanks a lot again to all of you.

On Tuesday, June 14, 2016 at 11:10:15 AM UTC+8, Roc Wang wrote:
>
> Hi, 
>
>   it seems trilinos-12.6.3 doen't work for dealii.8.4.1. Just wondering 
> which version is good for dealii8.4.1? thanks.
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

[Wolfgang Bangerth]

On 06/14/2016 10:44 AM, Denis Davydov wrote:

As for the solution transfers, i would speculate that you need two
different ones: one for the main field you solve for and another
for quadrature data.


Correct. You need to define a separate DoFHandler for the DG field to 
which you project your quadrature point data. Then you need a separate 
(second) SolutionTransfer object for the second DoFHandler.


Best
 W.
--

Wolfgang Bangerth   email:bange...@math.tamu.edu
www: http://www.math.tamu.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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

[dealii . group]

>
>
> Thank you. Is this class available in the release version 8.3.0?
>
I set separate dofhandler (history_dof_handler)  for FE_DGQ. I just 
> confusing because the history_field vector has different size respect to 
> the main solution vector and it also has different DoFHandler. I attached 
> dof_handler defined for the main problem for FE_Q, to the 
> solution_transfer, so I am confusing how we can use this solution_transfer 
> which defined based on the dof_handler to transfer history_field?!
>

On Tuesday, June 14, 2016 at 7:26:06 PM UTC+4:30, Denis Davydov wrote:
>
> Hi
>
> On 14 Jun 2016, at 16:51, dealii...@gmail.com  wrote:
>
> Hi
>
> I am solving a problem using Trilinos and parallel::distributed 
> 
>  triangulation. 
> It works well. I transfer solution data and old solution data during mesh 
> refinement as follow,
>
>
>   std::vector x(2);
>
>   x[0] = _solution;
>
>   x[1] = _old;
>
>
>   parallel::distributed::SolutionTransfer
>
> solution_transfer(dof_handler);
>
>
>   solution_transfer.prepare_for_coarsening_and_refinement(x);
>
>
>   triangulation.execute_coarsening_and_refinement();
>
>   setup_system();
>
>
>   LA::MPI::BlockVector tmp_v(partition);
>
>   LA::MPI::BlockVector tmp_vv(partition);
>
>   std::vector tmp(2);
>
>   tmp[0] = 
>
>   tmp[1] = _2;
>
>
>   solution_transfer.interpolate(tmp);
>
>   solution_old = tmp_2;
>
>
> Now, I also want to transfer quadrature point's data to the new mesh based 
> on the method explained at the end of step-18 in which we can use a 
> discontinuous field that matches the values in the quadrature points. Now 
> using this method, I can prepare a global field (history_field) and I want 
> to transfer the history_field vector as usual using the SolutionTransfer 
> 
>  class.
>
> Can I attach the history_field with relevant_solution and solution_old 
> during solution transfer method? 
>
> yes, you can. You would need to set up a separate DoFHandler and FE_DGQ 
> etc for this.
> Alternatively, you can use this class to achieve the same: 
> https://www.dealii.org/developer/doxygen/deal.II/classparallel_1_1distributed_1_1ContinuousQuadratureDataTransfer.html
>
>
> Or I must use two solution transfers, the first one for relevant_solution 
> and solution using dof_handler and the second one for history_field using 
> history_dof_handler? I can't understand this step and the relation between 
> FE_DGQ for history field and FE_Q for the main problem during refinement.
>
>
> The filed you solve for in FEM and the auxiliary field for transferring 
> are not related.
>
> Regards,
> Denis.
>
>  
>  
>
> On Monday, June 2, 2014 at 3:57:20 AM UTC+4:30, Wolfgang Bangerth wrote:
>>
>> On 05/26/2014 03:43 AM, Ali Karrech wrote: 
>> > Hi Timo, 
>> > Thank you very much. Your comments helped me. I though it's a good idea 
>> to 
>> > share the solution that I adopted with other users. Apart from 
>> > solution_transfer, I used FESystem and FE_DGQ to create a Finite 
>> element field 
>> > which has the same number of degrees as my integration points. I used 
>> this new 
>> > field (which I called history) to store my variable and recuperate them 
>> when I 
>> > need it. It seems to work properly so far. 
>>
>> This is what I would have done as well. There is a function 
>> FETools::compute_projection_from_quadrature_points_matrix that can help 
>> you 
>> with this, in case you haven't already found it. 
>>
>> Best 
>>   W. 
>>
>> -- 
>>  
>> Wolfgang Bangerth   email:bang...@math.tamu.edu 
>>  www: http://www.math.tamu.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 a topic in the 
> Google Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit 
> https://groups.google.com/d/topic/dealii/OSaxZJR19W8/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to 
> dealii+un...@googlegroups.com .
> For more options, visit https://groups.google.com/d/optout.
>
>
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] which version of trilinos for dealii-8.4.1

[Matthias Maier]

On Mon, Jun 13, 2016, at 22:42 CDT, Matthias Maier  wrote:

> On Mon, Jun 13, 2016, at 22:10 CDT, Roc Wang  wrote:
>
>> Hi, 
>>
>>   it seems trilinos-12.6.3 doen't work for dealii.8.4.1. Just wondering 
>> which version is good for dealii8.4.1? thanks.
>
> That's a bit surprising. I currently do not observe any issues with
> 12.6.2. Let me investigate.

Please have a look at [1] for configuration and installation
instructions.

Best,
Matthias

[1] http://dealii.org/8.4.1/external-libs/trilinos.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.
For more options, visit https://groups.google.com/d/optout.

[deal.II] Re: Elasticity

[Denis Davydov]

It appears that you are solving small strain non-linear elasticity. Then 
you need to study step-44 and linearize your problem.

Regards,
Denis.

On Tuesday, June 14, 2016 at 5:07:05 PM UTC+2, benhour@gmail.com wrote:
>
> Dear all,
> Thanks very much for your response and help. If I can make it simple, in 
> my equation, the lame constants in stress-strain tensor are related to 
> another function that changes with solution vector. For instance, consider 
> a situation in which you have a constant rank-2 tensor for initial strain 
> with this assumption that the indices of the matrix are multiplied by a 
> function of solution vectors. My question is that how can I define such a 
> matrix in deall ii because of the fact that all of my equations are almost 
> related to this interpolation function. Looking forward  to hearing from 
> you. I do really appreciate that.
>
> Best,
> Benhour
>
> On Monday, June 13, 2016 at 8:30:40 AM UTC-5, benhour@gmail.com wrote:
>>
>> Dear all,
>>
>> In step 18, the elasticity problem has been solved with this assumption 
>> that there was not initial strain in the equation. I want to add the 
>> initial strain with consideration of the fact that this tensor is a 
>> function of the solution. In other words, this tensor should be 
>> multiplied by a function that is related to the solution vectors. It  would 
>> be very kind of you if you help what is the best way to introduce this  
>> additional 
>> term in the equation. I do really appreciate your kindness and your time 
>> you can answer to my question. Looking
>> forward to hearing from you.
>>
>> Thanks,
>> Benhour
>>
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

[deal.II] Re: Elasticity

[benhour . amirian66]

Dear all,
Thanks very much for your response and help. If I can make it simple, in my 
equation, the lame constants in stress-strain tensor are related to another 
function that changes with solution vector. For instance, consider a 
situation in which you have a constant rank-2 tensor for initial strain 
with this assumption that the indices of the matrix are multiplied by a 
function of solution vectors. My question is that how can I define such a 
matrix in deall ii because of the fact that all of my equations are almost 
related to this interpolation function. Looking forward  to hearing from 
you. I do really appreciate that.

Best,
Benhour

On Monday, June 13, 2016 at 8:30:40 AM UTC-5, benhour@gmail.com wrote:
>
> Dear all,
>
> In step 18, the elasticity problem has been solved with this assumption 
> that there was not initial strain in the equation. I want to add the 
> initial strain with consideration of the fact that this tensor is a 
> function of the solution. In other words, this tensor should be 
> multiplied by a function that is related to the solution vectors. It  would 
> be very kind of you if you help what is the best way to introduce this  
> additional 
> term in the equation. I do really appreciate your kindness and your time 
> you can answer to my question. Looking
> forward to hearing from you.
>
> Thanks,
> Benhour
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

[dealii . group]

Hi

I am solving a problem using Trilinos and parallel::distributed 

 triangulation. 
It works well. I transfer solution data and old solution data during mesh 
refinement as follow,


  std::vector x(2);

  x[0] = _solution;

  x[1] = _old;


  parallel::distributed::SolutionTransfer

 solution_transfer(dof_handler);


  solution_transfer.prepare_for_coarsening_and_refinement(x);


  triangulation.execute_coarsening_and_refinement();

  setup_system();


  LA::MPI::BlockVector tmp_v(partition);

  LA::MPI::BlockVector tmp_vv(partition);

  std::vector tmp(2);

  tmp[0] = 

  tmp[1] = _2;


  solution_transfer.interpolate(tmp);

  solution_old = tmp_2;


Now, I also want to transfer quadrature point's data to the new mesh based 
on the method explained at the end of step-18 in which we can use a 
discontinuous field that matches the values in the quadrature points. Now 
using this method, I can prepare a global field (history_field) and I want 
to transfer the history_field vector as usual using the SolutionTransfer 
 class.

Can I attach the history_field with relevant_solution and solution_old 
during solution transfer method? Or I must use two solution transfers, the 
first one for relevant_solution and solution using dof_handler and the 
second one for history_field using history_dof_handler? I can't understand 
this step and the relation between FE_DGQ for history field and FE_Q for 
the main problem during refinement.

 

 

On Monday, June 2, 2014 at 3:57:20 AM UTC+4:30, Wolfgang Bangerth wrote:
>
> On 05/26/2014 03:43 AM, Ali Karrech wrote: 
> > Hi Timo, 
> > Thank you very much. Your comments helped me. I though it's a good idea 
> to 
> > share the solution that I adopted with other users. Apart from 
> > solution_transfer, I used FESystem and FE_DGQ to create a Finite element 
> field 
> > which has the same number of degrees as my integration points. I used 
> this new 
> > field (which I called history) to store my variable and recuperate them 
> when I 
> > need it. It seems to work properly so far. 
>
> This is what I would have done as well. There is a function 
> FETools::compute_projection_from_quadrature_points_matrix that can help 
> you 
> with this, in case you haven't already found it. 
>
> Best 
>   W. 
>
> -- 
>  
> Wolfgang Bangerth   email:bang...@math.tamu.edu 
>  
>  www: http://www.math.tamu.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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Transfer of quadrature point history data after mesh refinement.

[Denis Davydov]

Hi

> On 14 Jun 2016, at 16:51, dealii.gr...@gmail.com wrote:
> 
> Hi
> 
> I am solving a problem using Trilinos and parallel::distributed 
> 
>  triangulation. It works well. I transfer solution data and old solution data 
> during mesh refinement as follow,
> 
> 
> 
>   std::vector x(2);
> 
>   x[0] = _solution;
> 
>   x[1] = _old;
> 
> 
> 
>   parallel::distributed::SolutionTransfer
> 
> 
> solution_transfer(dof_handler);
> 
> 
> 
>   solution_transfer.prepare_for_coarsening_and_refinement(x);
> 
> 
> 
>   triangulation.execute_coarsening_and_refinement();
> 
>   setup_system();
> 
> 
> 
>   LA::MPI::BlockVector tmp_v(partition);
> 
>   LA::MPI::BlockVector tmp_vv(partition);
> 
>   std::vector tmp(2);
> 
>   tmp[0] = 
> 
>   tmp[1] = _2;
> 
> 
> 
>   solution_transfer.interpolate(tmp);
> 
>   solution_old = tmp_2;
> 
> 
> 
> Now, I also want to transfer quadrature point's data to the new mesh based on 
> the method explained at the end of step-18 in which we can use a 
> discontinuous field that matches the values in the quadrature points. Now 
> using this method, I can prepare a global field (history_field) and I want to 
> transfer the history_field vector as usual using the SolutionTransfer 
>  class.
> 
> Can I attach the history_field with relevant_solution and solution_old during 
> solution transfer method?
> 
yes, you can. You would need to set up a separate DoFHandler and FE_DGQ etc for 
this.
Alternatively, you can use this class to achieve the same: 
https://www.dealii.org/developer/doxygen/deal.II/classparallel_1_1distributed_1_1ContinuousQuadratureDataTransfer.html
 



> Or I must use two solution transfers, the first one for relevant_solution and 
> solution using dof_handler and the second one for history_field using 
> history_dof_handler? I can't understand this step and the relation between 
> FE_DGQ for history field and FE_Q for the main problem during refinement.
> 
> 

The filed you solve for in FEM and the auxiliary field for transferring are not 
related.

Regards,
Denis.

>  
>  
> 
> On Monday, June 2, 2014 at 3:57:20 AM UTC+4:30, Wolfgang Bangerth wrote:
> On 05/26/2014 03:43 AM, Ali Karrech wrote: 
> > Hi Timo, 
> > Thank you very much. Your comments helped me. I though it's a good idea to 
> > share the solution that I adopted with other users. Apart from 
> > solution_transfer, I used FESystem and FE_DGQ to create a Finite element 
> > field 
> > which has the same number of degrees as my integration points. I used this 
> > new 
> > field (which I called history) to store my variable and recuperate them 
> > when I 
> > need it. It seems to work properly so far. 
> 
> This is what I would have done as well. There is a function 
> FETools::compute_projection_from_quadrature_points_matrix that can help you 
> with this, in case you haven't already found it. 
> 
> Best 
>   W. 
> 
> -- 
>  
> Wolfgang Bangerth   email:bang...@ <>math.tamu.edu 
>  
>  www: http://www.math.tamu.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 a topic in the Google 
> Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit 
> https://groups.google.com/d/topic/dealii/OSaxZJR19W8/unsubscribe 
> .
> To unsubscribe from this group and all its topics, send an email to 
> dealii+unsubscr...@googlegroups.com 
> .
> For more options, visit https://groups.google.com/d/optout 
> .

-- 
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.
For more options, visit https://groups.google.com/d/optout.

Re: [deal.II] Re: casting a 4th order symm tensor to Tensor

[Anup Basak]

Hello Jean,

Thanks for the reply. Yes, I have already done what you suggested. But I
think that casting operation should be incorporated.

Many thanks,

Regards,
Anup.

On Tue, Jun 14, 2016 at 12:22 AM, Jean-Paul Pelteret 
wrote:

> Hi Anup,
>
> It looks like there's no constructor from which one can create a rank 4
> tensor from a rank 4 symmetric tensor, nor an assignment operator. I've
> opened up an issue on GitHub here
> . It looks like you'll just
> have to copy out the entries manually for now.
>
> J-P
>
>
> On Tuesday, June 14, 2016 at 1:07:56 AM UTC+2, Anup Basak wrote:
>>
>> Hello all,
>>
>> I want to cast a 4th order SymmetricTensor<4, dim> to a Tensor<4, dim>
>> and have used the following command:
>>
>> Tensor<4, dim> tens = Tensor<4, dim>(StandardTensors::II);
>>
>> where StandardTensors::II is the fourth order symmetric unit tensor
>> defined in
>> StandardTensors class (same as that defined in Step 44). However, I get
>> error saying that
>>
>> cannot convert ‘const value_type {aka const dealii::Tensor<1, 6,
>> double>}’ to ‘double’ in assignment
>>  t[d][d] = data[d];
>>
>> I shall be thankful if someone can point out where it is going wrong.
>>
>> Thanks and regards,
>> Anup.
>>
> --
> 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.
> For more options, visit https://groups.google.com/d/optout.
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

[deal.II] Re: which version of trilinos for dealii-8.4.1

[Denis Davydov]

Hi Roc,

I just built [1] deal.II 8.4.1. with Trilinos 12.6.3 using 
Spack https://github.com/LLNL/spack which also includes
as tests running Step-40 with Trilinos and (1) iterative solver (2) direct 
Superlu_Dist (3) direct MUMPS. 
So as far as I can tell it works ok.

[1] I will create a PR on Spack to update Trilinos to 12.6.3. soon.

Regards,
Denis.

On Tuesday, June 14, 2016 at 5:10:15 AM UTC+2, Roc Wang wrote:
>
> Hi, 
>
>   it seems trilinos-12.6.3 doen't work for dealii.8.4.1. Just wondering 
> which version is good for dealii8.4.1? thanks.
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.

[deal.II] Re: which version of trilinos for dealii-8.4.1

['Uwe Köcher' via deal . II User Group]

Dear Roc

The version 12.4.2 (shipped by default in candi https://github.com/dealii/candi 
) works.

Can you specify what you mean by "doesn't work"?

Best
  Uwe

On Tuesday, June 14, 2016 at 5:10:15 AM UTC+2, Roc Wang wrote:
>
> Hi, 
>
>   it seems trilinos-12.6.3 doen't work for dealii.8.4.1. Just wondering 
> which version is good for dealii8.4.1? thanks.
>

-- 
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.
For more options, visit https://groups.google.com/d/optout.