Lines Matching refs:nc
325 PetscInt nc = 0; in IJacobian() local
332 col[nc].i = i; in IJacobian()
333 col[nc].j = j; in IJacobian()
334 val[nc++] = 1.0; in IJacobian()
340 col[nc].i = i + 1; in IJacobian()
341 col[nc].j = j; in IJacobian()
342 val[nc++] = fwc / user->dx; in IJacobian()
343 col[nc].i = i; in IJacobian()
344 col[nc].j = j + 1; in IJacobian()
345 val[nc++] = -user->disper_coe / user->dy; in IJacobian()
346 col[nc].i = i; in IJacobian()
347 col[nc].j = j; in IJacobian()
348 val[nc++] = -fwc / user->dx + fthetac + user->disper_coe / user->dy; in IJacobian()
350 col[nc].i = i + 1; in IJacobian()
351 col[nc].j = j; in IJacobian()
352 val[nc++] = fwc / user->dx; in IJacobian()
353 col[nc].i = i; in IJacobian()
354 col[nc].j = j - 1; in IJacobian()
355 val[nc++] = user->disper_coe / user->dy; in IJacobian()
356 col[nc].i = i; in IJacobian()
357 col[nc].j = j; in IJacobian()
358 val[nc++] = -fwc / user->dx + fthetac - user->disper_coe / user->dy; in IJacobian()
360 col[nc].i = i - 1; in IJacobian()
361 col[nc].j = j; in IJacobian()
362 val[nc++] = -fwc / user->dx; in IJacobian()
363 col[nc].i = i; in IJacobian()
364 col[nc].j = j + 1; in IJacobian()
365 val[nc++] = -user->disper_coe / user->dy; in IJacobian()
366 col[nc].i = i; in IJacobian()
367 col[nc].j = j; in IJacobian()
368 val[nc++] = fwc / user->dx + fthetac + user->disper_coe / user->dy; in IJacobian()
370 col[nc].i = i - 1; in IJacobian()
371 col[nc].j = j; in IJacobian()
372 val[nc++] = -fwc / user->dx; in IJacobian()
373 col[nc].i = i; in IJacobian()
374 col[nc].j = j - 1; in IJacobian()
375 val[nc++] = user->disper_coe / user->dy; in IJacobian()
376 col[nc].i = i; in IJacobian()
377 col[nc].j = j; in IJacobian()
378 val[nc++] = fwc / user->dx + fthetac - user->disper_coe / user->dy; in IJacobian()
380 col[nc].i = i + 1; in IJacobian()
381 col[nc].j = j; in IJacobian()
382 val[nc++] = fwc / user->dx; in IJacobian()
383 col[nc].i = i; in IJacobian()
384 col[nc].j = j + 1; in IJacobian()
385 val[nc++] = -user->disper_coe / (2 * user->dy); in IJacobian()
386 col[nc].i = i; in IJacobian()
387 col[nc].j = j - 1; in IJacobian()
388 val[nc++] = user->disper_coe / (2 * user->dy); in IJacobian()
389 col[nc].i = i; in IJacobian()
390 col[nc].j = j; in IJacobian()
391 val[nc++] = -fwc / user->dx + fthetac; in IJacobian()
393 col[nc].i = i - 1; in IJacobian()
394 col[nc].j = j; in IJacobian()
395 val[nc++] = -fwc / user->dx; in IJacobian()
396 col[nc].i = i; in IJacobian()
397 col[nc].j = j + 1; in IJacobian()
398 val[nc++] = -user->disper_coe / (2 * user->dy); in IJacobian()
399 col[nc].i = i; in IJacobian()
400 col[nc].j = j - 1; in IJacobian()
401 val[nc++] = user->disper_coe / (2 * user->dy); in IJacobian()
402 col[nc].i = i; in IJacobian()
403 col[nc].j = j; in IJacobian()
404 val[nc++] = fwc / user->dx + fthetac; in IJacobian()
406 col[nc].i = i + 1; in IJacobian()
407 col[nc].j = j; in IJacobian()
408 val[nc++] = fwc / (2 * user->dx); in IJacobian()
409 col[nc].i = i - 1; in IJacobian()
410 col[nc].j = j; in IJacobian()
411 val[nc++] = -fwc / (2 * user->dx); in IJacobian()
412 col[nc].i = i; in IJacobian()
413 col[nc].j = j + 1; in IJacobian()
414 val[nc++] = -user->disper_coe / user->dy; in IJacobian()
415 col[nc].i = i; in IJacobian()
416 col[nc].j = j; in IJacobian()
417 val[nc++] = user->disper_coe / user->dy + fthetac; in IJacobian()
419 col[nc].i = i + 1; in IJacobian()
420 col[nc].j = j; in IJacobian()
421 val[nc++] = fwc / (2 * user->dx); in IJacobian()
422 col[nc].i = i - 1; in IJacobian()
423 col[nc].j = j; in IJacobian()
424 val[nc++] = -fwc / (2 * user->dx); in IJacobian()
425 col[nc].i = i; in IJacobian()
426 col[nc].j = j - 1; in IJacobian()
427 val[nc++] = user->disper_coe / user->dy; in IJacobian()
428 col[nc].i = i; in IJacobian()
429 col[nc].j = j; in IJacobian()
430 val[nc++] = -user->disper_coe / user->dy + fthetac; in IJacobian()
437 col[nc].i = i - 1; in IJacobian()
438 col[nc].j = j; in IJacobian()
439 val[nc++] = c1; in IJacobian()
440 col[nc].i = i + 1; in IJacobian()
441 col[nc].j = j; in IJacobian()
442 val[nc++] = -c1; in IJacobian()
443 col[nc].i = i; in IJacobian()
444 col[nc].j = j - 1; in IJacobian()
445 val[nc++] = c3 + c5; in IJacobian()
446 col[nc].i = i; in IJacobian()
447 col[nc].j = j + 1; in IJacobian()
448 val[nc++] = -c3 + c5; in IJacobian()
449 col[nc].i = i; in IJacobian()
450 col[nc].j = j; in IJacobian()
451 val[nc++] = -2 * c5 - a; in IJacobian()
453 PetscCall(MatSetValuesStencil(Jpre, 1, &row, nc, col, val, INSERT_VALUES)); in IJacobian()