say "there
cannot be any "holes" in the indices. Use AOCreateMapping() or
AOCreateMappingIS() if you wish to have "holes" in the indices." That's
why I am using AOCreateMapping in my code.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of En
ype-bound procedures. Previously this didn't bother it, but it looks
like it does now.
Is its complaint legitimate? or perhaps a compiler bug? (this is using
gcc 6.3.0)
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, N
nceTest(snes, convergence, context, &
1
Error: Actual argument at (1) to assumed-type dummy is of derived type
with type-bound or FINAL procedures
Cheers, Adrian
On Oct 17, 2018, at 5:21 PM, Adrian Croucher wrote:
hi
A colleague has just re
intel.com/en-us/forums/intel-fortran-compiler-for-linux-and-mac-os-x/topic/798744
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
name to return type***
***
Any ideas?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
a bug?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
ct.
This does compile, but doesn't look like it is actually doing anything-
i.e it's not giving the same behaviour as if I set "-mat_fd_type ds" on
the command line. It just gives the default wp behaviour.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department o
ian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
alling
MatMFFDSetType() with argument MATMFFD_DS instead of "ds"?
I will try debugging into MatMFFDSetType() and seeing what it's doing
when I pass in "ds".
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
On 24/06/16 15:58, Adrian Croucher wrote:
I will try debugging into MatMFFDSetType() and seeing what it's doing
when I pass in "ds".
I did this. In MatMFFDSetType() (mffd.c: 116), it executes
ierr =
PetscObjectTypeCompare((PetscObject)mat,MATMFFD,&match);CHKERRQ(ier
TMFFD (=
"mffd").
Does that mean I should be setting the matrix type to MATMFFD to get
this to work?
- Adrian
On 24/06/16 16:35, Adrian Croucher wrote:
On 24/06/16 15:58, Adrian Croucher wrote:
I will try debugging into MatMFFDSetType() and seeing what it's doing
when I pa
I started going astray.
How do I do the equivalent of setting the option "-mat_fd_type ds" in my
code? There is MatMFFDSetType() for MATMFFD, but doesn't seem to be
anything like a plain MatFDSetType().
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engin
t may possibly not matter any
more about wp or ds.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
.
Thanks, that appears to work fine for me.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 84611
missing?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
On 11/11/16 15:10, Barry Smith wrote:
Fixed in master and next. We forgot to add it for fortran
Thanks
Barry
You should just be able to do a git pull and have it working, no need to
reconfigure or recompile.
Yes, that's done the trick- thanks.
- Adrian
--
Dr A
appears to be complaining about the $EndElements line not matching
properly, but it looks ok to me.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
testmesh.msh
Description: Mesh model
.
Would it be possible to make DMPlexCreateGmsh() a bit more robust about
the line ending at the end of the file? Gmsh itself isn't bothered by it.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc..
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
dge+1)%(dim+1))*dim+d]);}
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
moment. But at some point I will need it.
Okay, when you need it, we will make it work.
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
hi Matt,
On 27/01/17 03:04, Matthew Knepley wrote:
On Wed, Jan 25, 2017 at 11:12 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi Matt,
A couple of years ago we discussed adding support for 6-node wedge
cells in DMPlex- the hurdle being that at present
hi
Currently the function MatFDColoringGetPerturbedColumns() doesn't apear
to be accessible from Fortran. I get an 'undefined reference' error when
I try to compile.
Is there just a missing interface?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department
which return arrays like this, e.g. ISGetIndicesF90().
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
PETSc 3.7.x
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
- Error Message
--
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
hanks very much, that seems to have fixed it.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
hi,
If I try to use the function DMGetLabelName() from Fortran I get an
'undefined reference' error while compiling.
Is there a missing interface?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
emai
On 29/06/17 05:59, Matthew Knepley wrote:
On Tue, Jun 27, 2017 at 6:39 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi,
If I try to use the function DMGetLabelName() from Fortran I get
an 'undefined reference' error while compiling.
Is
need to include petscdmlabel.h
as well.
Is that a bug or a feature?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
ves access to everything.
Also I'm not sure of exactly how to set up the ilocal and iremote array
parameters for this function in Fortran. How should they be declared- as
pointer arrays?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Unive
age in doing it that way?
However, at the moment I can't, because it appears the Fortran interface
for DMCreateDefaultSF() is missing. I would have thought the interface
for that should be simple enough, as it just takes a DM and a couple of
PetscSections.
- Adrian
Adrian Croucher w
On 25/09/17 13:18, Adrian Croucher wrote:
Alternatively, it occurred to me I could probably just use
DMCreateDefaultSF() to create the point SF on the new DM- would there
be any disadvantage in doing it that way?
However, at the moment I can't, because it appears the Fortran
inte
On 25/09/17 13:37, Adrian Croucher wrote:
Actually it looks like DMGetDefaultSF() does already work in Fortran
though, and could be used to do the same thing.
Do you think that would be a reasonable way to do it?
Hmm, I just tried it and it looks like the default SF produced by this
s so that the depth stratum bounds are
updated correctly. I also shift the end_interior value for each depth
stratum, so that code relying on DMPlexGetHybridBounds() should still work.
Would you expect copying the SF to work in this case?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
ndices, and also
update the remote root index for each leaf (presumably I can use the
original SF and PetscSFBcastBegin/ PetscSFBcastEnd to do that).
If you agree, then I will need working Fortran interfaces to the
PetscSFGetGraph/ PetscSFSetGraph functions, which are missing at
present. Are
hi
On 28/09/17 15:34, Adrian Croucher wrote:
So I think I will need to use PetscSFSetGraph() after all, so I can
increase the number of roots, update the leaf point indices, and also
update the remote root index for each leaf (presumably I can use the
original SF and PetscSFBcastBegin
hi
On 03/08/17 13:10, Matthew Knepley wrote:
On Wed, Aug 2, 2017 at 8:06 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
I haven't tried this myself, but two of my colleagues are
reporting problems accessing DMLabel from Fortran, when using the
latest &#
st don't understand how this
should be done.
- Adrian
On 28/09/17 15:34, Adrian Croucher wrote:
hi
On 28/09/17 04:18, Matthew Knepley wrote:
Okay, I think this should be easy to solve.
First a little bit about SF. There are two parts to the
specification. You have the communication part
ctory and rebuilding everything
(including make allfortranstubs) but again it didn't help.
I'm a bit surprised that custom Fortran bindings should produce anything
in include/petsc/finclude/ftn-auto/, I thought that would only be for
the auto-generated stuff that comes out of b
imilar "petsc/finclude/petscsf.h"?
It looks like PetscSF is already defined inside
petsc/finclude/petscis.h. It's not clear to me which bits of the SF
interface need to be in their own files and which are lumped in with the
IS interface.
- Adrian
--
Dr Adrian Croucher
Senior Rese
ms to me that it would be cleaner to do so.
So do you think the SF stuff in petsc/finclude/petscis.h should be taken
out, and put into a new petsc/finclude/petscsf.h ?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email:
hi
On 16/10/17 15:16, Jed Brown wrote:
Adrian Croucher writes:
So do you think the SF stuff in petsc/finclude/petscis.h should be taken
out, and put into a new petsc/finclude/petscsf.h ?
I think that's desirable for symmetry with include/petscsf.h, but it
isn't importan
hi
On 19/10/17 06:45, Matthew Knepley wrote:
On Tue, Oct 17, 2017 at 11:35 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
So, now I'm trying to add Fortran bindings for PetscSFBcastBegin()
and PetscSFBcastEnd().
From the C side I have added the fol
rray1dCreate((void*) ilocal, PETSC_INT, 1, *nleaves, lptr
PETSC_F90_2PTR_PARAM(lptrd)); if (*ierr) return;
+ *ierr = F90Array2dCreate((void*) iremote, PETSC_INT, 1, 2, 1,
*nleaves, rptr PETSC_F90_2PTR_PARAM(rptrd));
+}
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
it in /vec/f90-mod/ftn-auto-interfaces/petscpetscsf.h90), but the
one for PetscSFGetGraph() is missing for some reason. Any clues?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +6
hi
On 24/10/17 13:55, Adrian Croucher wrote:
The stub for PetscSFSetGraph() is getting automatically generated (I
can see it in /vec/f90-mod/ftn-auto-interfaces/petscpetscsf.h90), but
the one for PetscSFGetGraph() is missing for some reason. Any clues?
Oh, I see, all it needs to turn on
n't really seem to do that.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
hi,
On 11/10/17 22:18, Matthew Knepley wrote:
On Tue, Oct 10, 2017 at 9:07 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi
On 03/08/17 13:10, Matthew Knepley wrote:
On Wed, Aug 2, 2017 at 8:06 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz&g
(pilocal));
Other than that I think it is all go. (It even rightfully complains now
if I try to pass MPI_INT instead of MPI_INTEGER into PetscSFBcastBegin()
in my Fortran calling code.)
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
/impls/basic/sfbasic.c:497, which is:
DEF_PackCmp(int)
Dunno what that means though...
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
0)
at /home/acro018/software/PETSc/code/src/vec/is/sf/impls/basic/sfbasic.c:497
497 DEF_PackCmp(int)
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
On 01/11/17 16:20, Jed Brown wrote:
Adrian Croucher writes:
On 01/11/17 09:50, Smith, Barry F. wrote:
Please send the full traceback. Cut and paste
It's not really giving me much, only:
Program received signal SIGSEGV, Segmentation fault.
0x7efe8574d6a9 in Pack_int_1 (n=2,
hi
On 31/10/17 11:21, Adrian Croucher wrote:
I just pulled the latest next branch and found that the PetscSF stuff
doesn't appear to work in Fortran anymore for me.
I think I've found the trouble- after you put in the check to see if
Fortran has type(*) support, I needed to r
e do the following (in all your git clones):
git checkout master
git fetch
git branch -D next
git checkout next
I just did this, and now the PetscSFNode type is unavailable from
Fortran again. Barry added Fortran support for this a couple of weeks
ago, but now it seems to be gone.
- Adrian
--
x27;s fixed it.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
for the block
version the bounds should not be scaled.
I've just created a pull request
(acroucher/fix-IS-global-to-local-mapping-block) with a suggested fix.
- Adrian
On 16/11/17 11:52, Adrian Croucher wrote:
I actually attached the wrong test program last time- I've attached
sEnd = 0,0 and a default (or maybe null)
value for ft?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
n
/home/acro018/software/PETSc/code/src/vec/vec/interface/vector.c
[1]PETSC ERROR: #13 User provided function() line 0 in User file
#11 User provided function() line 0 in User file
On 13/02/18 17:13, Matthew Knepley wrote:
On Mon, Feb 12, 2018 at 10:57 PM, Adrian Croucher
mailto:a.crouc...@auckl
hi Matt,
On 20/02/18 12:32, Matthew Knepley wrote:
On Mon, Feb 19, 2018 at 4:33 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi Matt
I tried what you suggested and modified
DMPlexGetFieldType_Internal() so that all processes have the same
field ty
hi,
I have created a pull request (acroucher/dmplex-hdf5-vecview) with the
two fixes I came up with to prevent HDF5 DMPlex VecView() crashing when
one process has no cells on it.
Cheers, Adrian
On 22/02/18 17:12, Adrian Croucher wrote:
hi Matt,
On 20/02/18 12:32, Matthew Knepley wrote
and configured with --with-fortran-interfaces=1. At the
top of my code I have #include .
I noticed the fortran interfaces to some other functions (e.g.
DMPlexCreateFromDAG) also seemed to be missing in older versions (e.g.
the 'maint' branch)- is this one perhaps still to be added too?
troid(2), centroid(3)
end do
call DMDestroy(dm, ierr); CHKERRQ(ierr)
call PetscFinalize(ierr); CHKERRQ(ierr)
deallocate(numPoints, coneSize, cones, coneOrientations, vertexCoords)
end program main
--
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611
27;(a, i2, a, f8.4, a, 3(f8.4, 1x))'), &
'cell: ', i, ' volume: ', vol, ' centroid: ', &
pcentroid(1), pcentroid(2), pcentroid(3)
end do
call DMDestroy(dm, ierr); CHKERRQ(ierr)
call PetscFinalize(ierr); CHKERRQ(ierr)
deallocate
The interface needed to be declared. I pushed the fix, and your
example to PETSc (ex3f90).
Great, that's fixed it. The weird behaviour concerning initializing the
centroid and normal variables is gone now too. Thanks very much for your
help- much appreciated.
- Adrian
--
Dr A
a DMPlex that's been
distributed using DMPlexDistribute()? DMView just seems to give me the
cell counts as if the DMPlex was still all on rank 0. There doesn't seem
to be anything in the docs about DMViewFromOptions() yet, so I'm just
guessing what it's for from the examples.
, coneOrientations, vertexCoords)
end program main
--
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611
Scalar, pointer:: p(:)
However the compiler complains:
Error: Rank mismatch in argument 'array' at (1) (scalar and rank-1)
It's not liking the 'part' argument, but I'm not sure how it's supposed
to be passed in. Any clues?
Cheers, Adrian
--
Dr Adrian Croucher
De
)
call MPI_Comm_rank(comm, rank); CHKERRQ(ierr)
do c = cStart, cEnd-1
call PetscSectionGetOffset(cell_section, c, off, ierr); CHKERRQ(ierr)
part(off) = rank
end do
call VecRestoreArrayF90(partition, part, ierr); CHKERRQ(ierr)
end subroutine create_partition_vector
--
Dr Adrian
have I just
messed it up?
I also came across another function that seems to be lacking a fortran
interface: DMPlexConstructGhostCells(). Not totally sure if I need that,
but again was just trying to follow what is done in TS ex 33.
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering
that works out
I'll use it as a basis for solving my real problem.
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611
1
call DMPlexComputeCellGeometryFVM(dm, f, facearea, pcentroid,
pnormal, ierr)
write(*,'(i2, 1x, 4f10.4)'), f, facearea, centroid
end do
call DMDestroy(dm, ierr); CHKERRQ(ierr)
call PetscFinalize(ierr); CHKERRQ(ierr)
end program main
--
Dr Adrian Croucher
Department of
fixed it.
The problems I was seeing last week with VecView() writing to VTK seem
to have gone now too. Wunderbar!
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611
d cell types (like 6-node wedge) that have faces
with different numbers of nodes? We use wedge cells quite a bit, so it
would limit the usefulness of DMPlex if it can't support them. Unless
there is some kind of workaround?
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering
works first, and it probably
isn't top of my priority list either at present- I can continue getting
my code together without support for wedge cells at the moment. But at
some point I will need it.
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering Science
University of Au
cross
any more missing ones. I couldn't spot anything obvious in the git log
about how previous ones have been added though.
Cheers, Adrian
--
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611
working as yet. So I'm
a bit stuck. Should this approach work with DMPLex, or not?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
orf.o
CC linux-gnu-c-opt/obj/src/ksp/pc/impls/pbjacobi/pbjacobi.o
CC linux-gnu-c-opt/obj/src/ksp/pc/impls/mat/pcmat.o
CC linux-gnu-c-opt/obj/src/ksp/pc/impls/cp/cp.o
/home/acro018/software/PETSc/code/src/dm/f90-mod/petscdmmod.F:35.11:
use petscdm
1
--
Dr A
hi Matt,
On 13/05/14 00:54, Matthew Knepley wrote:
On Sun, May 11, 2014 at 8:46 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi, I just checkout the latest updates to the next branch, but am
getting an error on compiling PETSc- see below. Looks like it'
hi again
Oh, it might be just a typo in petscdmdef.h, line 21- should be
PetscLimiter, not PetscFLimiter ?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
On 13/05/14 11:39, Adrian Croucher wrote:
hi again
Oh, it might be just a typo in petscdmdef.h, line 21- should be
PetscLimiter, not PetscFLimiter ?
That fixes the problem with compiling PETSc. I am getting some warnings
now compiling TS ex11 though:
/home/acro018/software/PETSc/code/src
l *, const PetscReal *, const Petsc
Scalar *, const PetscScalar *, PetscScalar *, void *)’ but argument is
of type ‘RiemannFunction’
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.n
W() the IDs are
{100,101,200,300}. It looks like these IDs define where the BCs are
applied, but what do these numbers actually refer to?
Thanks!
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
te
hi
On 21/05/14 22:43, Matthew Knepley wrote:
On Wed, May 21, 2014 at 12:46 AM, Jed Brown <mailto:j...@jedbrown.org>> wrote:
Adrian Croucher mailto:a.crouc...@auckland.ac.nz>> writes:
> In the 'PhysicsCreate' functions there are calls to
DMPlexAddBoun
) function?
Such a thing doesn't appear to exist at the moment, so I'm just
wondering if it is in the pipeline, or if there is already some other
way to do this. I see there is a DMTSSetIFunctionLocal(), would that work?
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Depar
by signals sent by mpirun (as reported here).
--
[des108:00604] 1 more process has sent help message help-mpi-api.txt /
mpi-abort
[des108:00604] Set MCA parameter "orte_base_help_aggregate" to 0 to see
all help / error
code in C to avoid this issue, but will need to get it going in
fortran eventually.
Okay, let me know if it works. I will finish going through the ex11
changes by the end of the week.
Yes, that's fixed the problem. Thanks!
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Dep
nload-chaco --download-ctetgen
[0]PETSC ERROR: #1 PetscFVCreate() line 1356 in
/home/acro018/software/PETSc/code/src/dm/dt/interface/dtfv.c
--
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engine
$PETSC_DIR/$PETSC_ARCH/conf/sowing
and reconfigure using
./$PETSC_ARCH/conf/reconfigure-$PETSC_ARCH.py
OK thanks, I just tried that, but no luck unfortunately- am getting the
exact same error message.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Uni
s the same include as in my little code snippet. The
only thing different I can see is the format- I used free format and
ex3f90 seems to be fixed format?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auck
heers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
Looks like there are other missing Fortran interfaces:
DMSetNumFields()
DMSetField()
PetscFVSetNumComponents()
PetscFVSetSpatialDimension()
I have a feeling I might find more yet...
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of
dices()
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
ne get around this problem? I couldn't find any example
problems using Fortran modules.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
ian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
On 05/06/14 23:14, Matthew Knepley wrote:
On Thu, Jun 5, 2014 at 12:56 AM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
OK, I found one way around this- if in test.F90 I replace
Can you try just using the finclude/petscdef.h90 headers in the module?
Aha, yes tha
Riemann
function, to get access to the cells supporting the face and their
centroid positions- but the way the Riemann function interface is
currently set up doesn't appear to allow for that.
Or have I got the wrong idea about how one would go about implementing this?
- Adrian
--
Dr A
uot;y", ierr); CHKERRQ(ierr)
call PetscSectionSetUp(section, ierr); CHKERRQ(ierr)
call DMSetDefaultSection(dm, section); CHKERRQ(ierr)
call PetscSectionDestroy(section, ierr); CHKERRQ(ierr)
end subroutine setup_dm_section
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engi
On 09/08/14 06:43, Matthew Knepley wrote:
In 'next' I have changed the signature:
Aha, thanks, I hadn't noticed that.
Cheers, Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.crouc...@auckland
1 - 100 of 166 matches
Mail list logo
