appears to be on the "release" branch as well as "main" - i.e.
I can see it if I look at the commits for the "release" branch. So it
looked to me like a minor release. How are users supposed to interpret
those tags?
- Adrian
--
Dr Adrian Croucher
Senior Research Fel
ains that the latest patch update 3.21.4 doesn't
have PetscObjectIsNull in it.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email:a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
On 31/07/24 1:47 pm, Adrian Croucher wrote:
I just checked the file include/petsc/finclude/petscsys.h, and it
contained PetscObjectIsNull() until yesterday, but then it disappeared.
It is back in there again in commit 96fb8409 immediately after the
3.21.4 tag. But if you check out 3.21.4
I just checked the file include/petsc/finclude/petscsys.h, and it
contained PetscObjectIsNull() until yesterday, but then it disappeared.
- Adrian
On 31/07/24 1:27 pm, Adrian Croucher wrote:
On 31/07/24 1:25 pm, Barry Smith wrote:
You need to have
#include
in the file above the
On 31/07/24 1:25 pm, Barry Smith wrote:
You need to have
#include
in the file above the usage. That is where it is defined.
I have:
#include
Should that not still work?
- Adrian
Barry
On Jul 30, 2024, at 8:03 PM, Adrian Croucher
wrote:
This Message Is From an External
t week.
I deleted the build dirs for both PETSc and my code and did completely
clean builds of both, but got the same error.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email: a.crouc...@auck
6f90.F90 for the routine.
>
> Barry
>
>> On Jun 26, 2024, at 10:52 PM, Adrian Croucher wrote:
>>
>> PS, I also just tried adding that line to src/dm/imps/plex/tests/ex26f90.F90 using an earlier version of PETSc (before the Fortran interface changes) and that did run
PS, I also just tried adding that line to src/dm/imps/plex/tests/ex26f90. F90 using an earlier version of PETSc (before the Fortran interface changes) and that did run ok. - Adrian -- Dr Adrian Croucher Senior Research Fellow Department of Engineering
ZjQcmQRYFpfptBannerStart
--download-ptscotch --download-chaco --download-hypre
[0]PETSC ERROR: #1 PetscSectionGetValueLayout() at
/home/acro018/software/PETSc/code/src/vec/is/section/interface/section.c:1777
[0]PETSC ERROR: #2 ex26f90.F90:278
Abort(85) on node 0 (rank 0 in comm 16): application called
MPI_Abort(MPI_C
_utils.F90:401
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email:a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
;t see. You can try replacing SNESSetJacobian() in your code with SNESSetJacobianNoInterface() and see what happens.
>
>Barry
>
>
>
>
>> On Jun 21, 2024, at 12:56 AM, Adrian Croucher wrote:
>>
>> hi Barry,
>>
>> On 21/06/24 10:54 am, Barry Smith wrote:
>&g
x1f
and run it with -snes_fd_coloring ok.
However when I rebuilt my code in a new build dir against your PETSc
branch, it built ok but I still got the runtime error "Wrong type of
object: Parameter # 5" in SNESComputeJacobianDefaultColor().
- Adrian
--
Dr Adrian Croucher
Senior Rese
/2024-06-14/fix-fortransnessetjacobianmatfdcoloring)
Thanks, yes, I get the same result from git show -q.
So it looks a bit like the latest fix broke the previous one.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of
get the changes
Barry
On Jun 16, 2024, at 7:30 PM, Adrian Croucher
wrote:
hi Barry,
On 15/06/24 7:56 am, Barry Smith wrote:
Thanks for reporting the problem. It was an oversight on my part
when working on the Fortran interfaces.
It should work now in
https://urldefense.us/v3/__htt
024, at 1:29 AM, Adrian Croucher
wrote:
This Message Is From an External Sender
This message came from outside your organization.
hi,
I just tried updating to the latest PETSc 'main' branch. Since I'm using
Fortran I had to modify a few bits of my code for the new type checking
8/software/PETSc/code/src/snes/interface/snes.c:4755
[0]PETSC ERROR: #7 ../src/timestepper.F90:2295
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email: a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
also works even without the modification you just made to avoid
trouble with the depth label. So that was probably a red herring after all.
Thanks!
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
em
tached.
3. I ran stratification on this mesh, but I did not get an error:
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email:a.crouc...@auckland.ac.nz
tel: +64 (0)9 923 4611
0]PETSC ERROR: Petsc has generated inconsistent data
[0]PETSC ERROR: New depth 0 range [10,54) overlaps with depth 0 range
[10,54)
Would you expect that to work?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, Ne
st needs to be redistributed as well? But I guess that might not
always be the case.
Has anything changed in this part of the PETSc code recently that might
make this stop working?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / U
pointSF because the original
DMPlex points have been shifted to make space for the new cells, faces etc.
But if there is something not quite right in there, perhaps that could
cause trouble during the subsequent redistribution?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department
when computing flow
across it. (Not sure if that is the kind of topological operation you
meant.)
It has all been working fine with PETSc 3.15 since 2020! But because of
the low-level surgery I guess it is vulnerable to changes.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department
On 4/03/24 2: 36 pm, Adrian Croucher wrote: > I just re-tested my code with this "no cells in the interface" check > commented out of the PETSc code, and it appears to run fine. Whoops, no it doesn't. . . I checked more closely and there
ZjQcmQRYFpfptBannerStart
hi again, On 1/03/24 3: 41 pm, Adrian Croucher wrote: > > I am still getting problems with redistribution of MINC (dual > porosity) meshes. The function DMPlexCheckPointSF() is now raising > errors like this: > > [0]PETSC ERROR:
ZjQcmQRYFpfptBannerStart
commented "Check there
are no cells in interface", which is what is raising the error.
My understanding of the point SF is a bit rusty now but I would have
thought that for finite volume meshes with overlap 1 (which is what I
have), there should be cells in the interface?
Cheers,
ine but I think whatever you did
has probably fixed the issue.
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
;s trying these
URLs in sequence until it finds one that works?
I am using an older version of PETSc (3.15.5) as there is a PETSc bug in
more recent versions that affects my project and means I can't upgrade
until it is resolved.
Regards, Adrian
--
Dr Adrian Croucher
Senior Research
again.
Cheers, Adrian
On 1/12/23 9:00 am, Matthew Knepley wrote:
On Tue, Mar 28, 2023 at 6:12 PM Adrian Croucher
wrote:
> Yes, I think so. The partitioning groups cells, and does not
expect cells
> in the overlap. We could generalize that, which is probably
> the right th
.
Yes, and it might be useful to be able to disable that check for cells
in the overlap, for cases like this.
- 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
ap 1, then adding the
MINC cells etc. if needed, then redistributing (again with overlap 1).
It seems to have worked ok but could the fact that it's already been
distributed with overlap 1 be causing the error with redistribution?
- Adrian
--
Dr Adrian Croucher
Senior Research Fello
waiwera/waiwera/blob/55d26483b3392b70985fb25f998fbda3297c6154/src/mesh.F90#L3408
- 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
might be showing up a problem that was already
there? Any other clues as to how to approach 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
",PETSC_BOOL,&hdf5->defTimestepping,×tepping)
Is that not right?
Regards, Adrian
On 4/11/21 00:39, Matthew Knepley wrote:
On Tue, Nov 2, 2021 at 9:41 PM Adrian Croucher
wrote:
hi again,
I found the bug. At hdf5io.c:24, the arguments in the call to
PetscViewerHDF5ReadA
dire problem.
OK, I'll stop worrying and learn to love the bomb :-)
- 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
tascale systems, so MPI_COMM_CREATE_GROUP was introduced in
MPI-3 to alleviate this problem."
- 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
guess mpi_comm_split() might also have the advantage
that each rank only needs to know if it is in the group or not, rather
than needing an array of all participating ranks, as mpi_comm_create()
and mpi_comm_create_group() do.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Departme
a bit
misleading then.
Thanks!
- 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
27;ve read suggests you only
need to call it on processes in the group. It seems to work either way,
but you have to use the communicator a little differently.
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
emai
May 16, 2019, at 5:59 AM, Adrian Croucher via petsc-dev
mailto:petsc-dev@mcs.anl.gov>> wrote:
hi
A couple of years ago (see below) I was having problems reading in
GMSH files using DMPlexCreateFromFile(), if the MSH file didn't have
a newline on the end.
It was fixed at the ti
reate.c
[0]PETSC ERROR: #6 User provided function() line 0 in User file
On 17/01/17 4:32 PM, Matthew Knepley wrote:
On Mon, Jan 16, 2017 at 9:15 PM, Adrian Croucher
mailto:a.crouc...@auckland.ac.nz>> wrote:
hi Matt
On 17/01/17 15:52, Matthew Knepley wrote:
Its a DOS file an
es setting up a data layout, but if
so I am unclear as to what they are.
- Adrian
On 15/05/19 2:37 PM, Adrian Croucher wrote:
hi
I just tried upgrading from PETSc 3.10 to 3.11 and working my way
through all the changes that have broken things in my code.
I think I have most of them wor
looks like I am now supposed to use PetscFVSetComponentName() to set
field component names. Unfortunately there does not appear to be a
Fortran interface for this function yet?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New
test the new approach there also.
many thanks
Martin
--
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
module context_module
#include
use petsc
implicit none
pr
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
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
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
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
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
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
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
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
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
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
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
--
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
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,
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
/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
--
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
(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
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
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 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
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
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
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
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
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:
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
- 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
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
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
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
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
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
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
--
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
.
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..
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
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
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
.
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
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
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
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
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
1 - 100 of 166 matches
Mail list logo