On 12/08/2022 18:15, Matthew Knepley wrote:
This email was sent to you by someone outside the University.
You should only click on links or attachments if you are certain that
the email is genuine and the content is safe.
On Fri, Aug 12, 2022 at 9:27 AM David Scott wrote:
On 12/08/2022
On 12/08/2022 11:26, Matthew Knepley wrote:
This email was sent to you by someone outside the University.
You should only click on links or attachments if you are certain that the email
is genuine and the content is safe.
On Fri, Aug 12, 2022 at 2:21 AM Quynh Nguyen
mailto:nguyenquynh9...@gmai
document suggests 1MB is a good option.
Lawrence
On Fri, 26 Nov 2021 at 15:10, David Scott wrote:
Hello,
I am trying to write out HDF5 files. The program I have works all
right
up to a point but then fails when I double the size of the file that I
am trying to write out. I
Hello,
I am trying to write out HDF5 files. The program I have works all right
up to a point but then fails when I double the size of the file that I
am trying to write out. I have attached a file containing the error
message and another file containing the short program that I used to
generate i
Hello,
I am specifying an options file on the command line with '-options_file
'. I would like the program to write out to
stdout. Is there a simple way to do this?
I thought that I would be able interrogate the options database to
obtain the file name but, whilst the file is definitely being r
Matt,
I'll have a go at writing Xdmf files.
Am I right in thinking that VTK files are written sequentially?
Thanks,
David
On 09/07/2020 20:17, Matthew Knepley wrote:
On Thu, Jul 9, 2020 at 2:51 PM David Scott <mailto:d.sc...@epcc.ed.ac.uk>> wrote:
Hello,
I have wri
I should have said that I have created files with the suffix '.h5'. Was
that the right thing to do?
David
On 09/07/2020 19:50, David Scott wrote:
Hello,
I have written out a Vec using the (Fortran) HDF5 routines provided by
PETSc.
I now want to import the data into ParaView 5.8.0 a
Hello,
I have written out a Vec using the (Fortran) HDF5 routines provided by
PETSc.
I now want to import the data into ParaView 5.8.0 and display a contour.
I can read the data in using VisItPixieReader but the contour tool is
not available. Can I use the HDF5 files produced by PETSc 13.3.3 with
ding something?
Thanks,
Matt
Thanks,
David
On 18/02/2020 12:42, Matthew Knepley wrote:
On Tue, Feb 18, 2020 at 6:03 AM David Scott
mailto:d.sc...@epcc.ed.ac.uk>> wrote:
Hello,
I wish to solve a channel flow problem with different boundary
Hello,
I wish to solve a channel flow problem with different boundary
conditions. In the streamwise direction I may have periodic or
inlet/outlet BCs. I would like to make my code for the two cases as
similar as possible. If I use DM_BOUNDARY_PERIODIC then when performing
a linear solve the ghost
Hello,
In my code I have
call KSPSetComputeOperators(ksp_u, compute_matrix_u,
PETSC_NULL_OBJECT, ierr)
later I have something like
call KSPSolve(ksp_u, RHS_u_vec, u3_vec, ierr)
call DMCreateMatrix(da_u, B, ierr)
call compute_matrix_u(ksp_u, B, B, PETSC_NULL_INTEGER, ierr)
call KSPG
n 11, 2015 at 6:48 AM, David Scott <mailto:d.sc...@ed.ac.uk>> wrote:
Hello,
I am using MINRES with GAMG and have supplied various options
#PETSc Option Table entries:
-ksp_max_it 500
-ksp_monitor_true_residual
-log_summary
-mg_levels_ksp_max_it 2
Hello,
I am using MINRES with GAMG and have supplied various options
#PETSc Option Table entries:
-ksp_max_it 500
-ksp_monitor_true_residual
-log_summary
-mg_levels_ksp_max_it 2
-mg_levels_ksp_type chebyshev
-mg_levels_pc_type sor
-options_left
-pc_gamg_agg_nsmooths 1
-pc_gamg_threshold 0.03
-pc
On 12/12/2014 16:07, Matthew Knepley wrote:
On Fri, Dec 12, 2014 at 8:58 AM, David Scott <mailto:d.sc...@ed.ac.uk>> wrote:
How can I get GAMG to take into account the fact that I have more
than one degree of freedom per node when it is producing a coarser
grid? I am deali
How can I get GAMG to take into account the fact that I have more than
one degree of freedom per node when it is producing a coarser grid? I am
dealing with a problem in structural mechanics and have an unstructued
grid. The matrix is read from a file produced by another application.
The compon
On 03/05/2013 15:16, David Scott wrote:
> I have inserted
> call PetscViewerBinaryOpen(PETSC_COMM_WORLD, 'DifficultSystem', &
>FILE_MODE_WRITE, viewer, ierr)
> call KSPGetOperators(ksp, A, PETSC_NULL_OBJECT, &
>PETSC_NULL_INTEGER, ie
Barry,
On 23/04/2013 02:36, Barry Smith wrote:
>
> On Apr 22, 2013, at 11:26 AM, David Scott wrote:
>
>> Thanks for the suggestion.
>>
>> I had tried '-pc_type gamg -pc_gamg_agg_nsmooths 1' with an earlier version
>> of the code without success.
Hello Barry,
On 23/04/2013 02:36, Barry Smith wrote:
>David,
>
> We would be very interested in determining what is "going wrong" with
> the solver here since we hope to make it robust. Would it be possible for
> you to use a MatView() and VecView() on the matrix and the right hand sid
uces physically sensible results.
David
On 22/04/2013 13:12, Matthew Knepley wrote:
> On Mon, Apr 22, 2013 at 7:56 AM, David Scott <mailto:d.scott at ed.ac.uk>> wrote:
>
> Hello,
>
> I am working on a fluid-mechanical code to solve the two-phase
> Navier?Stokes e
.
If I use GMRES for the first few thousand time steps and then swap to CG
it does converge but the speed of execution is somewhat reduced.
I have attached relevant excerpts from the code.
Yours sincerely,
David Scott
--
Dr. D. M. Scott
Applications Consultant
Edinburgh Parallel Computing Centre
d. (Indeed, consider the family of problems in which you
>> extend the domain in that patch and apply _any_ boundary conditions in
>> the extended domain. All of those solutions are also solutions of your
>> problem with "second derivative is zero" on your "boundary&
ith implementation
> details.
>
>
> On Thu, Jan 10, 2013 at 9:03 AM, David Scott <mailto:d.scott at ed.ac.uk>> wrote:
>
> Hello,
>
> I am solving Poisson's equation (actually Laplace's equation in this
> simple test case) on a 3D structured gri
Hello,
I am solving Poisson's equation (actually Laplace's equation in this
simple test case) on a 3D structured grid. The boundary condition in the
first dimension is periodic. In the others there are Von Neumann
conditions except for one surface where the second derivative is zero. I
have sp
I wish to set boundary values in a structured grid. How do I do this if
I am using DMDA_BOUNDARY_GHOSTED? I thought that I should set the
boundary values in the function specified by a call to DMSetInitialGuess
but it appears to be passed a global vector.
I would prefer any example code to be i
24 matches
Mail list logo