<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font face="Courier New, Courier, monospace">Hi,<br>
I am about to use PETSc linear algebra. Two questions.<br>
<br>
Q1: </font><font face="Courier New, Courier, monospace">I have
the stiffness matrix and the rhs vector assembled (outside) and in
PETSc order. Do I still need to create a scatter:</font><font
face="Courier New, Courier, monospace"> "ierr =
VecScatterCreate(bb,isglobal,localX,islocal,&scatter)"?<br>
<br>
Q2: I have PETSc vectors and matrix set:</font><br>
<br>
<font face="Courier New, Courier, monospace">ierr =
VecCreate(MPI_COMM_WORLD,&bb);CHKERRQ(ierr);<br>
ierr = VecSetFromOptions(bb);CHKERRQ(ierr);CHKPT<br>
ierr = VecSetSizes(bb,Nvlocal,Nvglobal);CHKERRQ(ierr); <br>
ierr = MatCreate(PETSC_COMM_WORLD,&A);CHKERRQ(ierr);<br>
ierr = MatSetSizes(A,Nvlocal,</font><font face="Courier New,
Courier, monospace"><font face="Courier New, Courier, monospace">Nvlocal</font>,Nvglobal,Nvglobal);CHKERRQ(ierr);</font><font
face="Courier New, Courier, monospace"><br>
<br>
I passed the rhs myData to bb using </font><font face="Courier
New, Courier, monospace"><font face="Courier New, Courier,
monospace"><font face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace">VecSetValues</font></font></font>:</font><br>
<font face="Courier New, Courier, monospace"><font face="Courier
New, Courier, monospace"><font face="Courier New, Courier,
monospace">ierr = VecSetValues(bb,Nvlocal,vertices,myData,INSERT_VALUES);CHKERRQ(ierr);<br>
<br>
T</font>he matrix is in HYPRE sparse format. What is the
recommended way to pass it to A? <br>
<br>
</font></font><font face="Courier New, Courier, monospace">Thanks
in advance,<br>
--Ted<br>
</font>
</body>
</html>