Re: [deal.II] Error on running file with mounted volume

2020-05-29 Thread Shanmukh Sripada
Yes, the 52nd line in the parameters.in file says "set Print timing 
information with output = true".
Deleting this line solves the problem! 

Thank you for bringing this to my notice.

Kind regards,
Shanmukh 

On Friday, May 29, 2020 at 10:29:39 PM UTC+5:30, Wolfgang Bangerth wrote:
>
> On 5/29/20 5:15 AM, Shanmukh Sripada wrote: 
> > I get the following error when I run the "mpirun -n 1 main" command when 
> I 
> > open a container on docker using host volume (like this: docker run -it 
> -v 
> > (host/directory):(container directory) image. 
> > However, I don't get this error when I don't use the volume feature 
> (like 
> > this: docker run -it image). 
> > What could be the reason for this? 
>
> I have no idea what the reason for this might be, but what *is* in line 52 
> of 
> your input file? Does it say the "Print timing..." text the error message 
> is 
> talking about? If not, is there another file anywhere that has that text 
> in 
> line 52 and that would suggest what is going on? 
>
> 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/dae81cb6-7e56-4c0d-8c90-cc176cb983d5%40googlegroups.com.


Re: [deal.II] Error on running file with mounted volume

2020-05-29 Thread Wolfgang Bangerth

On 5/29/20 5:15 AM, Shanmukh Sripada wrote:
I get the following error when I run the "mpirun -n 1 main" command when I 
open a container on docker using host volume (like this: docker run -it -v 
(host/directory):(container directory) image.
However, I don't get this error when I don't use the volume feature (like 
this: docker run -it image).

What could be the reason for this?


I have no idea what the reason for this might be, but what *is* in line 52 of 
your input file? Does it say the "Print timing..." text the error message is 
talking about? If not, is there another file anywhere that has that text in 
line 52 and that would suggest what is going on?


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/1c8854f1-c593-535f-0931-e8df2f26a39c%40colostate.edu.


[deal.II] Please explain the projection method (tutorial step-35)

2020-05-29 Thread Alexey Ozeritskiy
Hello,

In step-35 we find u^{k+1} (diffusion step)
Then we find \phi^{k+1} (projection step) and correct p^{k+1} pressure. 
But why we don't correct u^{k+1} velocity? I think u^{k+1} must be 
solenoidal but it is not.

-- 
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/88d0faab-9653-4d91-9ace-cc8a13b5e79b%40googlegroups.com.


Re: [deal.II] problem measure() for quads

2020-05-29 Thread Jean-Paul Pelteret
Hi Andreas,

Thanks for being willing to make a pull request of this. The basic instructions 
for how to do so are listed in our FAQ, 
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#i-found-a-typo-or-a-bug-and-fixed-it-on-my-machine-how-do-i-get-it-included-in-dealii
 

but of course if you have any troubles just reach out to us and we’d be happy 
to help further.

Best,
Jean-Paul

> On 29 May 2020, at 08:55, Andreas Kyritsakis  wrote:
> 
> Thank you very much Wolfgang. Yes I would like to propose a patch. How is it 
> done? I am not familiar with the process.
> 
> On Thursday, May 28, 2020 at 7:31:04 PM UTC+3, Wolfgang Bangerth wrote:
> On 5/28/20 3:08 AM, Andreas Kyritsakis wrote: 
> > 
> > I have a problem with using the measure() function for quads in 3D. I am 
> > currently using version 9.1.1 and the function seems to check whether the 
> > face 
> > is planar and if not, it returns nan. The problem is that the tolerance for 
> > considering it planar is extremely low (1.e-24), which causes problems. If 
> > for 
> > example the mesh is read in from a standard .msh file, written by the 
> > GridOut 
> > class, we cannot expect to get such a precision in the vertices. I saw that 
> > in 
> > the 9.2 version this is partially fixed by handling the non-planar case by 
> > numerical intergration, but the tolerance has remained the same extremely 
> > low, 
> > which might result in unnecessary computational effort. 
> > 
> > I wonder whether the tolerance can be increased to something more 
> > reasonable 
> > (e.g. 1.e-10) or there is something I am missing that suggests such a low 
> > tolerance. 
> 
> It's a good question. I checked that the expression in question is 
> dimensionless and scaled correctly so that the mesh size h does not appear in 
> the size of the expression. I think it would be quite ok if you changed the 
> tolerance to 1e-10. 
> 
> Want to propose a patch for this? 
> 
> 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/231d1ecd-03a2-4cce-9590-ed1470fd5278%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/DB0C30AC-B940-4D80-9250-4C7AAA7E91CF%40gmail.com.


Re: [deal.II] problem measure() for quads

2020-05-29 Thread Andreas Kyritsakis
Thank you very much Wolfgang. Yes I would like to propose a patch. How is 
it done? I am not familiar with the process.

On Thursday, May 28, 2020 at 7:31:04 PM UTC+3, Wolfgang Bangerth wrote:
>
> On 5/28/20 3:08 AM, Andreas Kyritsakis wrote: 
> > 
> > I have a problem with using the measure() function for quads in 3D. I am 
> > currently using version 9.1.1 and the function seems to check whether 
> the face 
> > is planar and if not, it returns nan. The problem is that the tolerance 
> for 
> > considering it planar is extremely low (1.e-24), which causes problems. 
> If for 
> > example the mesh is read in from a standard .msh file, written by the 
> GridOut 
> > class, we cannot expect to get such a precision in the vertices. I saw 
> that in 
> > the 9.2 version this is partially fixed by handling the non-planar case 
> by 
> > numerical intergration, but the tolerance has remained the same 
> extremely low, 
> > which might result in unnecessary computational effort. 
> > 
> > I wonder whether the tolerance can be increased to something more 
> reasonable 
> > (e.g. 1.e-10) or there is something I am missing that suggests such a 
> low 
> > tolerance. 
>
> It's a good question. I checked that the expression in question is 
> dimensionless and scaled correctly so that the mesh size h does not appear 
> in 
> the size of the expression. I think it would be quite ok if you changed 
> the 
> tolerance to 1e-10. 
>
> Want to propose a patch for this? 
>
> 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/231d1ecd-03a2-4cce-9590-ed1470fd5278%40googlegroups.com.