Re: [petsc-users] Integrate PETSC with existing Fortran

2019-01-23 Thread Yaxiong Chen via petsc-users
Hi Matt, I changed optimal_mechanical_part.f90 to optimal_mechanical_part.F90. But I still received the following error: HiDACX/obj/debug/optimal_mechanical_part.o mpicc -Dqh_QHpointer -ansi -Wall -I/usr/local/include/libqhull -c /Users/yaxiong/Documents/HiDAC_Code_Git/HiDAC_Tao_NURBS_Pets

Re: [petsc-users] Integrate PETSC with existing Fortran

2019-01-23 Thread Matthew Knepley via petsc-users
On Wed, Jan 23, 2019 at 1:27 PM Yaxiong Chen wrote: > Hi Matt, > > > I tried to modify the structure of present code and use KSP in the main > program(optimal_mechanical_part). Now the makefile is as following : > > > #=

Re: [petsc-users] Integrate PETSC with existing Fortran

2019-01-23 Thread Yaxiong Chen via petsc-users
Hi Matt, I tried to modify the structure of present code and use KSP in the main program(optimal_mechanical_part). Now the makefile is as following : # OUT_DIR := HiDACX/ ## Define mode ifndef MODE .PHONY: debug

Re: [petsc-users] Printing parallel matrix

2019-01-23 Thread Jed Brown via petsc-users
Evan Um via petsc-users writes: > Dear PETSC users, > > I try to verify a matrix by printing a PETSC matrix and comparing its > elements with reference values. Below is my test code. > > It works well when a single process is used. The output file is created > quickly. In contrast, when multiple