[petsc-users] QCG and KSPSetPCSide

All, I am trying to use KSPSolve for a QCG solver using petsc 3.4.3 and it gives me the error of "PC does not have left symmetric apply!". Further reading on the web pages suggest I need to use KSPSetPCSide to set up a symmetric side for the solver. However, this does not seem to work. I am doi

[petsc-users] KSPSetComputeEigenvalues memory corruption

All, I am using KSPSetComputeEigenvalues (in petsc 3.4.3) and getting a memory corruption problem. In the manual it says: "Currently this option is not valid for all iterative methods" Which methods and could this be why there is a memory corruption error message? Carol Dr Carol Brickley BSc,P

Re: [petsc-users] EXTERNAL: Re: “PETSC ERROR: Unable to find requested PCtype hypre”

Hi, I thought I had but looking at the log summary it isn’t coming up that petsc was built properly with hypre. My script had a bug which caused it to ignore everything after the blas library entries. Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applied Computer Scie

[petsc-users] “PETSC ERROR: Unable to find requested PCtype hypre”

Hi, I am trying to run an executable built with petsc 3.5.3 and hypre 2.9.0b with flags “-pc_type hypre –ksp_type cg” but keep getting the error “PETSC ERROR: Unable to find requested PCtype hypre” Any ideas why this might be happening? Thanks Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS S

[petsc-users] petsc 3.4.3 to 3.5.3

All, I am switching to petsc 3.5.3 and finding some differences with 3.4.3. KSPSetOperators no longer has the MatStructure argument. I see for SAME_PRECONDITIONER that you now use KSPSetReusePreconditioner with flag set to tru before the call to set to true KSPSetOperators. Is there any need to

Re: [petsc-users] EXTERNAL: Re: Problems trying to use -log_summary

Thanks Couldn't get it to recognise the flag, -log_summary, but wrote to file using PetscLogPrintDetailed(PETSC_COMM_WORLD, 'petsc_log_summary.txt', petsc_err). Thanks Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applied Computer Science DS+T, AWE Aldermaston Readi

[petsc-users] Problems trying to use -log_summary

Hi, I am having problems trying to use -log_summary. I have added PetscLogBegin,PetscLogInitSummary and compiled with PETSC_USE_LOG. I then run with -log_summary but the code will not accept this flag, Any ideas? Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applie

[petsc-users] petsc 3.4.3

All I am using petsc 3.4.3 but finding it not working as quickly as it should. Any ideas? Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applied Computer Science DS+T, AWE Aldermaston Reading Berkshire RG7 4PR Direct: 0118 9855035 ___

Re: [petsc-users] EXTERNAL: Re: DMDAVecGetArrayF90

Excellent - just remembered that. I'll see if it works. Had taken it out to do my own thing. Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applied Computer Science DS+T, AWE Aldermaston Reading Berkshire RG7 4PR Direct: 0118 9855035 -Original Message- Fro

[petsc-users] DMDAVecGetArrayF90

Hi, I am using DMDAVecGetArrayF90 with a pointer in petsc 3.4.3. Is there anyway of finding out what the start and finish indices of the pointer coming back from this routine are per processor? I have declared: DM :: da Vec :: X Pets

[petsc-users] FW: EXTERNAL: Re: EXTERNAL: Re: F90 to petsc 3.4.3 interface - loss of accuracy

All, I am still having problems with data inconsistencies. We think it is loss of accuracy but am not sure. I have tried using 1 tolerance of 1e-12 for the ksp tolerance (ksp_rtol) but that doesn't seem to change matters. Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engin

[petsc-users] F90 to petsc 3.4.3 interface - loss of accuracy

Hi, Does anyone find loss of accuracy between F90 code and Petsc 3.4.3 C codes? I am noticing that some of my data is returning 9.9997e-13 instead of 1.e-12 for example. Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS Senior Software Engineer Applied Computer Science DS+T, AWE Alder

[petsc-users] VecGetArray issue

All, When I use VecGetArray it only seems to return half of the array? I am using it: VecGetArray (ptr->sol_x_loc, sol_v, iss, err) And sol_v only seems to come back with half the vector? I am using a F90 interface to the Petsc 3.4.3 routines. Carol Dr Carol Brickley BSc,PhD,ARCS,DIC,MBCS S

[petsc-users] F90 program interfaced with petsc 3.4.3 - inconsitencies

Hi, I am trying to runs a F90 program interfaced with petsc 3.4.3. When the routine DMDACreate2D is called, if I run with a debugger such as ddt, it shows me that in the petsc routines the DM array (inra) is populated, however when I return to the F90 code the DM array has one value and no long

[petsc-users] Error compiling f90 code with petsc 3.5.2

Hi, When I try to compile an F90 code and include the 3.5.2 petsc built under intelmpi-4.1.3 and intel13.1, I get: ...error #5082: Syntax error, found '::' when expecting one of: ( % [ : . = => DMDABoundaryType:: bx, by, bz ! boundary type .