NetSolve and PETSc

>> [mv, idx, jdx] = find(A);

>> nrows = size(A, 1);

% case: default

>> [y, its, ierr] = netsolve('itmeth’, mv, x, idx, jdx, nrows, tol, maxit);

% case: Jacobi preconditioner

>> [y, its, ierr] = netsolve('itmeth', mv, x, idx, jdx, nrows, tol, maxit,

'-pc_type jacobi');

% case: incomplete ILU preconditioner

>> [y, its, ierr] = netsolve('itmeth', mv, x, idx, jdx, nrows, tol, maxit,

'-pc_type ilu -pc_ilu_levels 4');

• Example of using packages remotely

Previous slide

Next slide

Back to first slide

View graphic version