MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 10 MatAssemblyEnd_SeqAIJ(): Matrix size: 25495389 X 10 storage space: 0 unneeded,254953890 used PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374784 Total number of mallocs used during MatSetValues calls =0 MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 504677 MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0 PetscFOpen(): Opening file /group/ml/data/tensor/nell/ PetscInitialize(): PETSc successfully started: number of processors = 1 Results and Errors with -info -mat-view-info: Ierr = MatMatMult(X1, CC, MAT_INITIAL_MATRIX, PETSC_DEFAULT, &M) CHKERRQ(ierr) Ierr = MatConvert(MC, MATBAIJ, MAT_INITIAL_MATRIX, &CC) CHKERRQ(ierr) Ierr = MatCreateMAIJ(C, J, &MC) CHKERRQ(ierr) Ierr = MatAssemblyEnd(C, MAT_FINAL_ASSEMBLY) CHKERRQ(ierr) Ierr = MatAssemblyBegin(C, MAT_FINAL_ASSEMBLY) CHKERRQ(ierr) Ierr = MatSetValues(C, 1, &k, 1, &r, &one, INSERT_VALUES) CHKERRQ(ierr) Ierr = MatCreateSeqAIJ(PETSC_COMM_SELF, K, R, R, NULL, &C) CHKERRQ(ierr) Create a matrix C (K x R) with all values 1 Ierr = PetscPrintf(PETSC_COMM_WORLD, "Setup Time: %2.1e \n", v1-v) CHKERRQ(ierr) Ierr = MatAssemblyEnd(X1, MAT_FINAL_ASSEMBLY) Ierr = MatAssemblyBegin(X1, MAT_FINAL_ASSEMBLY) Ierr = MatSetValues(X1, 1, &i, 1, &j, &val, INSERT_VALUES) CHKERRQ(ierr) Ierr = MatSeqAIJSetPreallocation(X1, 0, nnz) CHKERRQ(ierr) Ierr = MatSetType(X1, MATSEQAIJ) CHKERRQ(ierr) Ierr = MatSetBlockSizes(X1, I, J) CHKERRQ(ierr)
Ierr = MatSetSizes(X1, PETSC_DECIDE, PETSC_DECIDE, I, J*K) CHKERRQ(ierr) Ierr = MatCreate(PETSC_COMM_SELF, &X1) CHKERRQ(ierr) Is this memory problem, and do I have to change seqaij into mpiaij and use multi processors? Or do I have another methods fixing it? If you know the method, then please let me know it. However, I got errors such as out of memory and Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range. I need to multiply a large seqaij matrix(X1) and a maij(or baij) matrix(CC). Next message: Errors from large matrices.Previous message: How to write the dmma mesh into tecplot format.Errors from large matrices Joon hee Choi choi240 at
0 kommentar(er)
