Re: [deal.II] Output of Gauss point stress tensor

2020-07-30 Thread Muhammad Mashhood
Dear Prof. Wolfganag,
   Thank you for the suggestion. I ran 
debugger with the ".prm" file and was successful to backtrace the error and 
resolved it. After this, using the step-19 particles tutorial refered by 
you, so far I am also successful in creating and exporting out the 
particles in ".vtu" file to visualize their ids and location (which is 
quadrature points location in my case).

As a next step, After creating particle on quadrature point, I am now 
trying to assign it the quadrature point value which is "*SymmetricTensor<2, 
dim>*" in my case. *I am using following simplified function similar to the 
one in particles tutorial file of step-19:*



































*template void  Problem::create_particles(){  
next_unused_particle_id = 1;FEValues fe_values (fe, 
quadrature_formula, update_values | 
update_gradients | 
update_quadrature_points);for (const auto  : 
dof_handler.active_cell_iterators())  {
PointHistory *local_quadrature_points_history
   = reinterpret_cast 
*>(cell->user_pointer());  Assert 
(local_quadrature_points_history >=   
   _point_history.front(),
  ExcInternalError());   
Assert (local_quadrature_points_history <
   _point_history.back(),
  ExcInternalError());
   fe_values.reinit(cell);for (const unsigned int 
q_point :fe_values.quadrature_point_indices())
  {const Point location =
fe_values.quadrature_point(q_point);
Particles::Particle new_particle;
new_particle.set_location(location);
new_particle.set_reference_location(
  mapping.transform_real_to_unit_cell(cell, 
location));
new_particle.set_id(next_unused_particle_id);
SymmetricTensor<2, dim> total_strain = 
local_quadrature_points_history[q_point].old_strain;*









*
new_particle.set_properties(make_array_view(total_strain));
particle_handler.insert_particle(new_particle, 
cell);++next_unused_particle_id;  
}  }particle_handler.update_cached_numbers();}*

*But it gives the following error on running:*










*An error occurred in line <298> of file 
 in 
functionvoid dealii::Particles::Particle::set_properties(const dealii::ArrayView&) [with int 
dim = 3; int spacedim = 3]The violated condition was: property_pool != 
nullptrAdditional information: This exception -- which is used in many 
places in the library -- usually indicates that some condition which the 
author of the code thought must be satisfied at a certain point in an 
algorithm, is not fulfilled. An example would be that the first part of an 
algorithm sorts elements of an array in ascending order, and a second part 
of the algorithm later encounters an element that is not larger than the 
previous one.There is usually not very much you can do if you encounter 
such an exception since it indicates an error in deal.II, not in your own 
program. Try to come up with the smallest possible program that still 
demonstrates the error and contact the deal.II mailing lists with it to 
obtain help.*I tried to fix it by initializing and setting the property 
pool to the new_particle in my current function :

*Particles::PropertyPool propertypool(6); // six entries present in the 
symmetric tensor of second order*
*new_particle.set_property_pool(propertypool);*

*But still it stops with the following error:*







*An error occurred in line <531> of file 
 in 
functiondealii::ArrayView::value_type& 
dealii::ArrayView::operator[](std::size_t) 
const [with ElementType = double; MemorySpaceType = 
dealii::MemorySpace::Host; dealii::ArrayView::value_type = double; std::size_t = long unsigned int]The 
violated condition was: static_cast::type)>:: type>(i) < 
static_cast::type)>:: 
type>(n_elements)Additional information: Index 3 is not in the 
half-open range [0,3).*
I explored the mailing list and found someone from the community already 
had similar problem of assigning the properties to the particles where it 
was suggested to use following way of assigning the properties:



*for( auto iter=particleHandler.begin(); iter!=particleHandler.end(); 
++iter, ++part ) *





* iter->set_properties 
(make_array_view(interpolatedParticleQuantities.begin()+ part*ncomponents, 
interpolatedParticleQuantities.begin()+ 

Re: [deal.II] Unexpected behavior when using GridGenerator::subdivided_hyper_rectangle in parallel

2020-07-30 Thread Jimmy Ho
Hi Dr. Bangerth,

Thanks a lot for the clarification! They are really helpful!

Best,
Jimmy

On Thursday, July 30, 2020 at 11:47:21 AM UTC-5, Wolfgang Bangerth wrote:
>
> On 7/30/20 10:11 AM, Jimmy Ho wrote: 
> > 
> > As a follow-up question, upon calling compress(), will the local copy of 
> the 
> > system matrix on a specific processor get updated to contain information 
> from 
> > all other processors? In other words, if I print out the system matrix 
> from a 
> > particular processor after calling compress(), is that the same global 
> system 
> > matrix that the linear solver is solving? 
>
> Each processor only stores a subset of the rows of the matrix. During 
> assembly, each processor computes a number of matrix entries, but it can 
> not 
> compute all entries for the rows of the matrix it owns -- some need 
> contributions from other processes. The call to compress() makes sure the 
> contributions from these other processes are sent to the one that owns 
> these 
> rows of the matrix. 
>
> In any case, after compress(), the rows each processor owns are correct -- 
> but 
> each processor doesn't know anything about the rows of the matrix it 
> doesn't own. 
>
> Best 
>   W. 
>
>
> -- 
>  
> 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/b424e816-cd22-4854-b432-d9d9d745ae8do%40googlegroups.com.


Re: [deal.II] Unexpected behavior when using GridGenerator::subdivided_hyper_rectangle in parallel

2020-07-30 Thread Wolfgang Bangerth

On 7/30/20 10:11 AM, Jimmy Ho wrote:


As a follow-up question, upon calling compress(), will the local copy of the 
system matrix on a specific processor get updated to contain information from 
all other processors? In other words, if I print out the system matrix from a 
particular processor after calling compress(), is that the same global system 
matrix that the linear solver is solving?


Each processor only stores a subset of the rows of the matrix. During 
assembly, each processor computes a number of matrix entries, but it can not 
compute all entries for the rows of the matrix it owns -- some need 
contributions from other processes. The call to compress() makes sure the 
contributions from these other processes are sent to the one that owns these 
rows of the matrix.


In any case, after compress(), the rows each processor owns are correct -- but 
each processor doesn't know anything about the rows of the matrix it doesn't own.


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/9a27440e-2864-afb2-d21c-cecdcc229305%40colostate.edu.


Re: [deal.II] Unexpected behavior when using GridGenerator::subdivided_hyper_rectangle in parallel

2020-07-30 Thread Jimmy Ho
Hi Dr. Bangerth,

As a follow-up question, upon calling compress(), will the local copy of 
the system matrix on a specific processor get updated to contain 
information from all other processors? In other words, if I print out the 
system matrix from a particular processor after calling compress(), is that 
the same global system matrix that the linear solver is solving?

Thanks a lot for clarifying!

Best,
Jimmy

On Wednesday, July 29, 2020 at 9:55:38 AM UTC-5, Wolfgang Bangerth wrote:
>
>
> Jimmy, 
>
> > A minimum example to reproduce this is attached. When the mesh is built 
> using 
> > GridGenerator::hyper_cube or GridGenerator::subdivided_hyper_rectangle 
> with 
> > subsequent refinement, the program works as expected. When the same mesh 
> is 
> > generated using GridGenerator::subdivided_hyper_rectangle without any 
> > subsequent refinement, some entries in the global stiffness matrix (the 
> nodes 
> > on the right hand side of the mesh, in this case, using 2 processors) do 
> not 
> > get updated. The example output of stiffness matrices when using one and 
> two 
> > processors are also attached for reference. 
> > 
> > So my question is, is this the expected behavior of the function? If so, 
> why 
> > is that the case? 
>
> There is still a lot of stuff in the program that could be remove, 
> including 
> all of the comments, to make it substantially smaller and easier to 
> understand. 
>
> I don't know whether there is a bug, but here is a suggestion: The finite 
> element solution u_h(x) that results from the linear system should be the 
> same 
> independent of the partitioning. But the order of degrees of freedom may 
> be 
> different, and consequently the matrix may not be the same -- it should 
> only 
> be the same up to some column and row permutation. Have you verified that 
> the 
> *solution function* (not the solution vector) that results is the same 
> independent of the number of processors? 
>
> Best 
>   W. 
>
> -- 
>  
> 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/39a0b8f8-c0f4-437d-b70c-7eed4295b13co%40googlegroups.com.


Re: [deal.II] Unexpected behavior when using GridGenerator::subdivided_hyper_rectangle in parallel

2020-07-30 Thread Jimmy Ho
Hi Dr. Bangerth,

Thanks a lot for your guidance! I compared the solution in the vtu files 
using the minimal example above, they are nearly identical. Looking back 
into the code, I am outputting the system matrix from processor 0, which 
probably only printed the part that it locally owns, hence there's a 
difference in the matrices when running in serial and parallel. I guess I 
rushed to conclusion too quickly.

Best,
Jimmy

On Wednesday, July 29, 2020 at 9:55:38 AM UTC-5, Wolfgang Bangerth wrote:
>
>
> Jimmy, 
>
> > A minimum example to reproduce this is attached. When the mesh is built 
> using 
> > GridGenerator::hyper_cube or GridGenerator::subdivided_hyper_rectangle 
> with 
> > subsequent refinement, the program works as expected. When the same mesh 
> is 
> > generated using GridGenerator::subdivided_hyper_rectangle without any 
> > subsequent refinement, some entries in the global stiffness matrix (the 
> nodes 
> > on the right hand side of the mesh, in this case, using 2 processors) do 
> not 
> > get updated. The example output of stiffness matrices when using one and 
> two 
> > processors are also attached for reference. 
> > 
> > So my question is, is this the expected behavior of the function? If so, 
> why 
> > is that the case? 
>
> There is still a lot of stuff in the program that could be remove, 
> including 
> all of the comments, to make it substantially smaller and easier to 
> understand. 
>
> I don't know whether there is a bug, but here is a suggestion: The finite 
> element solution u_h(x) that results from the linear system should be the 
> same 
> independent of the partitioning. But the order of degrees of freedom may 
> be 
> different, and consequently the matrix may not be the same -- it should 
> only 
> be the same up to some column and row permutation. Have you verified that 
> the 
> *solution function* (not the solution vector) that results is the same 
> independent of the number of processors? 
>
> Best 
>   W. 
>
> -- 
>  
> 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/285f4af3-3e01-4cd8-bcf0-9e6062a53967o%40googlegroups.com.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread Wolfgang Bangerth

On 7/30/20 3:59 AM, kaleem iqbal wrote:
I want to simulate FSI with elastic walls in bifurcation carotid arteries in 
2D. is it possible to use that geometry?


It is possible to use this geometry, but you need to generate a mesh for it. 
(See step-49 for mesh generation.)


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/6935fa53-0f02-121f-53dc-3489962b762f%40colostate.edu.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread kaleem iqbal
I want to simulate FSI with elastic walls in bifurcation carotid arteries
in 2D. is it possible to use that geometry?


On Thu, Jul 30, 2020 at 2:54 PM kaleem iqbal  wrote:

> thanks for the reply. simulation run
>
>
> On Thu, Jul 30, 2020 at 2:18 PM Luca Heltai  wrote:
>
>> Take a look at your parameter file. It is probably trying to write output
>> file to a directory that does not exist.
>>
>> Luca
>>
>> Il giorno 30 lug 2020, alle ore 10:58, kaleem iqbal <
>> kaleemiqb...@gmail.com> ha scritto:
>>
>> 
>> Dear Prof. Wolfgang;
>> During running step-70. I found the following error
>> Exception on processing:
>>
>> 
>> An error occurred in line <1360> of file
>> 
>> in function
>> void dealii::ParameterHandler::print_parameters(const string&,
>> dealii::ParameterHandler::OutputStyle) const
>> The violated condition was:
>> out
>> Additional information:
>> An input/output error has occurred. There are a number of reasons why
>> this may be happening, both for reading and writing operations.
>>
>> If this happens during an operation that tries to read data: First, you
>> may be trying to read from a file that doesn't exist or that is not
>> readable given its file permissions. Second, deal.II uses this error at
>> times if it tries to read information from a file but where the information
>> in the file does not correspond to the expected format. An example would be
>> a truncated file, or a mesh file that contains not only sections that
>> describe the vertices and cells, but also sections for additional data that
>> deal.II does not understand.
>>
>> If this happens during an operation that tries to write data: you may be
>> trying to write to a file to which file or directory permissions do not
>> allow you to write. A typical example is where you specify an output file
>> in a directory that does not exist.
>> 
>>
>> Aborting!
>> 
>> CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run'
>> failed
>> make[3]: *** [CMakeFiles/run] Error 1
>> CMakeFiles/Makefile2:131: recipe for target 'CMakeFiles/run.dir/all'
>> failed
>> make[2]: *** [CMakeFiles/run.dir/all] Error 2
>> CMakeFiles/Makefile2:138: recipe for target 'CMakeFiles/run.dir/rule'
>> failed
>> make[1]: *** [CMakeFiles/run.dir/rule] Error 2
>> Makefile:144: recipe for target 'run' failed
>> make: *** [run] Error 2
>> Please guide me
>> Regards
>> Kaleem iqbal
>>
>> --
>> 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/168928f6-65e2-40e2-924a-84bfe514e2bco%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/2E791A8D-2E6F-4FF2-A232-0ABDD68618BF%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/CAOHBfT1vFJOOp328GXce-TjH2tGFfTjnBPQF4c3P3C1%2BVVhizA%40mail.gmail.com.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread kaleem iqbal
thanks for the reply. simulation run


On Thu, Jul 30, 2020 at 2:18 PM Luca Heltai  wrote:

> Take a look at your parameter file. It is probably trying to write output
> file to a directory that does not exist.
>
> Luca
>
> Il giorno 30 lug 2020, alle ore 10:58, kaleem iqbal <
> kaleemiqb...@gmail.com> ha scritto:
>
> 
> Dear Prof. Wolfgang;
> During running step-70. I found the following error
> Exception on processing:
>
> 
> An error occurred in line <1360> of file
> 
> in function
> void dealii::ParameterHandler::print_parameters(const string&,
> dealii::ParameterHandler::OutputStyle) const
> The violated condition was:
> out
> Additional information:
> An input/output error has occurred. There are a number of reasons why
> this may be happening, both for reading and writing operations.
>
> If this happens during an operation that tries to read data: First, you
> may be trying to read from a file that doesn't exist or that is not
> readable given its file permissions. Second, deal.II uses this error at
> times if it tries to read information from a file but where the information
> in the file does not correspond to the expected format. An example would be
> a truncated file, or a mesh file that contains not only sections that
> describe the vertices and cells, but also sections for additional data that
> deal.II does not understand.
>
> If this happens during an operation that tries to write data: you may be
> trying to write to a file to which file or directory permissions do not
> allow you to write. A typical example is where you specify an output file
> in a directory that does not exist.
> 
>
> Aborting!
> 
> CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run' failed
> make[3]: *** [CMakeFiles/run] Error 1
> CMakeFiles/Makefile2:131: recipe for target 'CMakeFiles/run.dir/all' failed
> make[2]: *** [CMakeFiles/run.dir/all] Error 2
> CMakeFiles/Makefile2:138: recipe for target 'CMakeFiles/run.dir/rule'
> failed
> make[1]: *** [CMakeFiles/run.dir/rule] Error 2
> Makefile:144: recipe for target 'run' failed
> make: *** [run] Error 2
> Please guide me
> Regards
> Kaleem iqbal
>
> --
> 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/168928f6-65e2-40e2-924a-84bfe514e2bco%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/2E791A8D-2E6F-4FF2-A232-0ABDD68618BF%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/CAOHBfT2s7cTd5dMnJ%2BDYxS5TCv%3Dw2NmrODudCGMYt6VVkU4%3DCw%40mail.gmail.com.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread kaleem iqbal
Thanks heena. Simulation run now.

Regards
K Iqbal

-- 
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/CAOHBfT2DnbiHuJ%2BhTEqreHtbzcxM3%2B9vT3dXQDc7wANOJk1KBg%40mail.gmail.com.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread Luca Heltai
Take a look at your parameter file. It is probably trying to write output file 
to a directory that does not exist. 

Luca

> Il giorno 30 lug 2020, alle ore 10:58, kaleem iqbal  
> ha scritto:
> 
> 
> Dear Prof. Wolfgang; 
> During running step-70. I found the following error  
> Exception on processing: 
> 
> 
> An error occurred in line <1360> of file 
> 
>  in function
> void dealii::ParameterHandler::print_parameters(const string&, 
> dealii::ParameterHandler::OutputStyle) const
> The violated condition was: 
> out
> Additional information: 
> An input/output error has occurred. There are a number of reasons why 
> this may be happening, both for reading and writing operations.
> 
> If this happens during an operation that tries to read data: First, you may 
> be trying to read from a file that doesn't exist or that is not readable 
> given its file permissions. Second, deal.II uses this error at times if it 
> tries to read information from a file but where the information in the file 
> does not correspond to the expected format. An example would be a truncated 
> file, or a mesh file that contains not only sections that describe the 
> vertices and cells, but also sections for additional data that deal.II does 
> not understand.
> 
> If this happens during an operation that tries to write data: you may be 
> trying to write to a file to which file or directory permissions do not allow 
> you to write. A typical example is where you specify an output file in a 
> directory that does not exist.
> 
> 
> Aborting!
> 
> CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run' failed
> make[3]: *** [CMakeFiles/run] Error 1
> CMakeFiles/Makefile2:131: recipe for target 'CMakeFiles/run.dir/all' failed
> make[2]: *** [CMakeFiles/run.dir/all] Error 2
> CMakeFiles/Makefile2:138: recipe for target 'CMakeFiles/run.dir/rule' failed
> make[1]: *** [CMakeFiles/run.dir/rule] Error 2
> Makefile:144: recipe for target 'run' failed
> make: *** [run] Error 2
> Please guide me
> Regards
> Kaleem iqbal
> -- 
> 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/168928f6-65e2-40e2-924a-84bfe514e2bco%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/2E791A8D-2E6F-4FF2-A232-0ABDD68618BF%40gmail.com.


Re: [deal.II] Mesh-induced elastic anisotropy and distorting the quad. points as a way to palliate it

2020-07-30 Thread David F
Hi Sebastian,

first of all, sorry for my late reply. Thank you very much for comment, it 
certainly raises some very interesting points. I think the only thing which 
is still left to be explained is the fact that a triangular grid, which 
yields the right isotropic result, becomes anisotropic just by rearranging 
the elastic properties to form squares-like clusters. I think this case 
does not correspond to any of your 4 points since I do this without h->inf. 
In fact, each of those clusters are just 2 FE.

Best,
David.


On Saturday, 11 July 2020 06:19:50 UTC+2, Sebastian Stark wrote:
>
> Hi David,
>
> I'm neither an expert, nor do I know the literature well, but looking on 
> your pictures, I think, the situations you are studying are geometrically 
> anisotropic. Just plot the distribution of angles the faces of your 
> inhomogenities make with the x-axis. For the quad-case, you'll get two 
> discrete peaks at 0 and 90 degree. for the triangular case, you get 0, 45 
> and 90 degree. So, from this, the results do not seem surprising to me 
> (just consider the extreme case of cracks - if you have them only at 0 and 
> 90 degree oriented, this is unlikely to be isotropic). The fact that you 
> have randomly assigned elastic properties won't help to fix that.
>
> A few examples (in 2d):
>
> (1) equidistant circular inclusions in a matrix (matrix and inclusions two 
> different isotropic linearly elastic materials) -> this should be isotropic 
> if mesh size h->0
>
> (2) equidistant square inclusions in a matrix, all aligned with x-axis 
> (matrix and inclusions two different isotropic linearly elastic materials) 
> -> probably anisotropic if mesh size h->0
>
> (3) equidistant square inclusions in a matrix, random orientation of 
> inclusions (matrix and inclusions two different isotropic linearly elastic 
> materials) -> should be isotropic if mesh size h->0 and number of 
> orientations->infinity
>
> (4) equidistant square inclusions in a matrix, random orientation of 
> inclusions (random isotropic linearly elastic materials) -> should be 
> isotropic if mesh size h->0 and number of orientations->infinity
>
> Also to consider: In 3d, there are 21 elastic constants for a linearly 
> elastic material. In a mathematical 2d scenario, it should be 6. This 
> suggests that, in example (3), it is not strictly necessary to have random 
> orientation. Rather, a few (equally spaced) discrete orientations might be 
> good enough. If that's the case, how many does one need? I'm betting on 6, 
> not sure though. Alternatively, one could replace the square inclusions in 
> (2) by regular polygons and ask how many vertices the polygon needs for 
> isotropy. Again, I'm betting on 6.
>
> Related: Are there crystal structures with such a high degree of symmetry, 
> that they are elastically isotropic? For dielectric properties, a cubic 
> crystal is good enough already. But the dielectric tensor is rank 2 and the 
> elastic one rank 4. So you'll need much more symmetry in the crystal; and 
> considering that a crystal can have 6-fold rotational symmetry at most that 
> might be impossible.
>
> What I'm just noticing: Hexagonal crystals are elastically isotropic 
> perpendicular to the hexagonal axis. So, my bet on 6 might be good. And it 
> might explain your observation that the triangular elements are relatively 
> isotropic (though maybe not perfectly).
>
> I hope that gives you some input. If you have definitive answers to any of 
> the questions, I'm curious.
>
> Regards,
> Sebastian
>
>
> Am 11.07.20 um 00:12 schrieb David F:
>
> I have made a somewhat extensive study on his issue and prepared some 
> plots that will hopefully answer your questions, and also includes Bruno's 
> suggestion about distorting the mesh. The basic setup is: I sheared the 
> mesh along different orientations (see x-axis on the plots) and measured 
> the shear modulus (y-axis). I have repeated the random process of setting 
> the elastic properties many times to have good statistics (see errorbars on 
> the plots). Each element has an isotropic stiffness tensor with a Poisson 
> ratio of 1/3 and a shear modulus which is exponentially distributed with an 
> average of 10. I use linear shape functions unless otherwise stated. If the 
> picture are not big enough, you can find them in the links beneath them.
>
>
> *1) I change the resolution. *By this I don't mean just a mesh with a 
> bigger number of elements, but importantly each inhomogeneities is 
> represented by a bigger number of elements. Therefore, we solve problems 
> with exactly the same physical domain but with different resolution. In the 
> legend, n means the resolution of the inhomogeneities. For n=1 each 
> inhomogeneity is described by 1 FE. For n=2, by 2^2, and for n=4 by 4^2. We 
> can see that for shearing with angle 0 (see pictures on the bottom for 
> clarity) the shear modulus is minimum, while it is maximum for 45 degrees, 
> when the principal axes are aligned 

Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread heena patel
Hi Iqbal,
   Check this link below.
https://groups.google.com/forum/?utm_medium=email_source=footer#!searchin/dealii/step$2070$209.2$20version%7Csort:date/dealii/Y5Lgt9Gn0Zs/vNHTNDdaAwAJ


Regards,
Heena

On Thu, Jul 30, 2020 at 10:58 AM kaleem iqbal 
wrote:

> Dear Prof. Wolfgang;
> During running step-70. I found the following error
> Exception on processing:
>
> 
> An error occurred in line <1360> of file
> 
> in function
> void dealii::ParameterHandler::print_parameters(const string&,
> dealii::ParameterHandler::OutputStyle) const
> The violated condition was:
> out
> Additional information:
> An input/output error has occurred. There are a number of reasons why
> this may be happening, both for reading and writing operations.
>
> If this happens during an operation that tries to read data: First, you
> may be trying to read from a file that doesn't exist or that is not
> readable given its file permissions. Second, deal.II uses this error at
> times if it tries to read information from a file but where the information
> in the file does not correspond to the expected format. An example would be
> a truncated file, or a mesh file that contains not only sections that
> describe the vertices and cells, but also sections for additional data that
> deal.II does not understand.
>
> If this happens during an operation that tries to write data: you may be
> trying to write to a file to which file or directory permissions do not
> allow you to write. A typical example is where you specify an output file
> in a directory that does not exist.
> 
>
> Aborting!
> 
> CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run' failed
> make[3]: *** [CMakeFiles/run] Error 1
> CMakeFiles/Makefile2:131: recipe for target 'CMakeFiles/run.dir/all' failed
> make[2]: *** [CMakeFiles/run.dir/all] Error 2
> CMakeFiles/Makefile2:138: recipe for target 'CMakeFiles/run.dir/rule'
> failed
> make[1]: *** [CMakeFiles/run.dir/rule] Error 2
> Makefile:144: recipe for target 'run' failed
> make: *** [run] Error 2
> Please guide me
> Regards
> Kaleem iqbal
>
> --
> 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/168928f6-65e2-40e2-924a-84bfe514e2bco%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/CAG_OxbgHXT4P448t%3DZo24ncfAjd4XfaK5MtZ-esO6_LEJr7Hbg%40mail.gmail.com.


Re: [deal.II] Installation didn't give any errors but when I tried make test, it failed all tests

2020-07-30 Thread kaleem iqbal
Dear Prof. Wolfgang; 
During running step-70. I found the following error  
Exception on processing: 


An error occurred in line <1360> of file 

 
in function
void dealii::ParameterHandler::print_parameters(const string&, 
dealii::ParameterHandler::OutputStyle) const
The violated condition was: 
out
Additional information: 
An input/output error has occurred. There are a number of reasons why 
this may be happening, both for reading and writing operations.

If this happens during an operation that tries to read data: First, you may 
be trying to read from a file that doesn't exist or that is not readable 
given its file permissions. Second, deal.II uses this error at times if it 
tries to read information from a file but where the information in the file 
does not correspond to the expected format. An example would be a truncated 
file, or a mesh file that contains not only sections that describe the 
vertices and cells, but also sections for additional data that deal.II does 
not understand.

If this happens during an operation that tries to write data: you may be 
trying to write to a file to which file or directory permissions do not 
allow you to write. A typical example is where you specify an output file 
in a directory that does not exist.


Aborting!

CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run' failed
make[3]: *** [CMakeFiles/run] Error 1
CMakeFiles/Makefile2:131: recipe for target 'CMakeFiles/run.dir/all' failed
make[2]: *** [CMakeFiles/run.dir/all] Error 2
CMakeFiles/Makefile2:138: recipe for target 'CMakeFiles/run.dir/rule' failed
make[1]: *** [CMakeFiles/run.dir/rule] Error 2
Makefile:144: recipe for target 'run' failed
make: *** [run] Error 2
Please guide me
Regards
Kaleem iqbal

-- 
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/168928f6-65e2-40e2-924a-84bfe514e2bco%40googlegroups.com.