scaling in 4-core machine: unassembled structured

2009-11-20 Thread jaru...@ascomp.ch

 Getting better multicore performance requires changing the algorithms to
 better reuse L1 cache.  This means moving away from assembled matrices
 where possible and of course finding good preconditioners.

I do not know how to move away from assembled matrix. As I have to  
reset values to the matrix in each iteration, I oblige to call  
MatAssemblyBegin() and MatAssemblyEnd(). Is there other option to  
create and set values the matrix??

 High-order
 and fast multipole methods are good for this.

For example, please?


Jarunan






-- 
Jarunan Panyasantisuk
Development Engineer
ASCOMP GmbH, Technoparkstr. 1
CH-8005 Zurich, Switzerland
Phone : +41 44 445 4072
Fax   : +41 44 445 4075
E-mail: jarunan at ascomp.ch
www.ascomp.ch


suprelu

2009-11-20 Thread Satish Balay
works fine for me.

Perhaps you can retry after 'rm -rf externalpackages'

If the problem persists - send configure.log to petsc-maint at mcs.anl.gov

Satish

On Fri, 20 Nov 2009, Lukasz Kaczmarczyk wrote:

 Hello
 I try to configure petsc (p9) with SuperLU and I get error message 
 
 *
  UNABLE to CONFIGURE with GIVEN OPTIONS(see configure.log for 
 details):
 ---
 Error unzipping _d_SuperLU_DIST.tar.gz: Could not execute 'cd 
 /Users/likask/MyBuild/src/petsc-3.0.0-p9/externalpackages; gunzip 
 _d_SuperLU_DIST.tar.gz':
 
 gzip: _d_SuperLU_DIST.tar.gz: not in gzip format
 *
 
 Regards,
 Lukasz



Change orthogonalization option in fortran?

2009-11-20 Thread Barry Smith

   Sorry, we don't have the Fortran interfaces for these operations.  
You can use

   call PetscOptionsSet(- 
ksp_gmres_classicalgramschmidt,PETSC_NULL_CHARACTER,ierr)
   call PetscOptionsSet(- 
ksp_gmres_cgs_refinement_type,REFINE_IFNEEDED,ierr)

  before creating the KSP object.

Barry


On Nov 19, 2009, at 1:14 PM, hxie at umn.edu wrote:

 Hi,

 I want to change the method for orthogonalization for the default  
 ksp solver in fortran. I added the following in my code:
 --
   call  
 KSPGMRESSetOrthogonalization 
 (ksp,KSPGMRESClassicalGramSchmidtOrthogonalization,pterr)
   call  
 KSPGMRESSetCGSRefinementType 
 (ksp,KSP_GMRES_CGS_REFINEMENT_IFNEEDED,pterr)
 --

 And I got the following error message when compiling.
 --
 error #6404: This name does not have a type, and must have an  
 explicit type. [KSPGMRESCLASSICALGRAMSCHMIDTORTHOGONALIZATI]
  call  
 KSPGMRESSetOrthogonalization 
 (ksp,KSPGMRESClassicalGramSchmidtOrthogonalization,pterr)

 error #6404: This name does not have a type, and must have an  
 explicit type. [KSP_GMRES_CGS_REFINEMENT_IFNEEDED]
  call  
 KSPGMRESSetCGSRefinementType 
 (ksp,KSP_GMRES_CGS_REFINEMENT_IFNEEDED,pterr)
 --

 If I comment these two lines, the code is OK. Any idea for this?  
 Thanks.

 Bests,
 Hui



suprelu

2009-11-20 Thread Lukasz Kaczmarczyk
Thanks, for help. Strange enough It is working now. 

Regards,
Lukasz

On 20 Nov 2009, at 17:33, Barry Smith wrote:

 
  Please send configure.log to petsc-maint at mcs.anl.gov
 
  Looks like you have a corrupt tar.gz file; recommend deleting it and getting 
 it again.
 
   Barry
 
 On Nov 20, 2009, at 11:22 AM, Lukasz Kaczmarczyk wrote:
 
 Hello
 I try to configure petsc (p9) with SuperLU and I get error message
 
 *
UNABLE to CONFIGURE with GIVEN OPTIONS(see configure.log for 
 details):
 ---
 Error unzipping _d_SuperLU_DIST.tar.gz: Could not execute 'cd 
 /Users/likask/MyBuild/src/petsc-3.0.0-p9/externalpackages; gunzip 
 _d_SuperLU_DIST.tar.gz':
 
 gzip: _d_SuperLU_DIST.tar.gz: not in gzip format
 *
 
 Regards,
 Lukasz
 
 

Lukasz Kaczmarczyk
Lecturer
Department of Civil Engineering,
University of Glasgow,
GLASGOW, G12 8LT
Tel: +44 141 3305348
email: likask at civil.gla.ac.uk
web: http://www.civil.gla.ac.uk/~kaczmarczyk/
web: http://code.google.com/p/yaffems/