PageRenderTime 24ms CodeModel.GetById 34ms RepoModel.GetById 0ms app.codeStats 0ms

/hpc_kernel_samples/sparse_linear_algebra/casestudies/pagerank-petsc/petsc-2.3.2-p10/src/mat/impls/rowbs/mpi/iccbs.c.html

https://gitlab.com/steinret/CodeVault
HTML | 101 lines | 86 code | 15 blank | 0 comment | 0 complexity | cd1ad49459a16d8c6e1afe5c2a4d7402 MD5 | raw file
    

  1. <center><a href="iccbs.c">Actual source code: iccbs.c</a></center><br>
    2. 

  2. 

  3. <html>
    4. 

  4. <head>
    5. 

  5. <title></title>
    6. 

  6. <meta name="generator" content="c2html 0.9.5">
    7. 

  7. <meta name="date" content="2007-03-29T00:23:47+00:00">
    8. 

  8. </head>
    9. 

  9. 

  10. <body bgcolor="#FFFFFF">
    11. 

  11. <pre width="80"><a name="line1">  1: </a><strong><font color="#228B22">#define PETSCMAT_DLL</font></strong>
    12. 

  12. 

  13. <a name="line3">  3: </a><font color="#B22222">/*</font>
    14. 

  14. <a name="line4">  4: </a><font color="#B22222">   Defines a Cholesky factorization preconditioner with BlockSolve95 interface.</font>
    15. 

    16. 

  16. <a name="line6">  6: </a><font color="#B22222">   Note that BlockSolve95 works with a scaled and permuted preconditioning matrix.</font>
    17. 

  17. <a name="line7">  7: </a><font color="#B22222">   If the linear system matrix and preconditioning matrix are the same, we then</font>
    18. 

  18. <a name="line8">  8: </a><font color="#B22222">   work directly with the permuted and scaled linear system:</font>
    19. 

  19. <a name="line9">  9: </a><font color="#B22222">      - original system:  Ax = b</font>
    20. 

  20. <a name="line10"> 10: </a><font color="#B22222">      - permuted and scaled system:   Cz = f, where</font>
    21. 

  21. <a name="line11"> 11: </a><font color="#B22222">             C = P D^{-1/2} A D^{-1/2}</font>
    22. 

  22. <a name="line12"> 12: </a><font color="#B22222">             z = P D^{1/2} x</font>
    23. 

  23. <a name="line13"> 13: </a><font color="#B22222">             f = P D^{-1/2} b</font>
    24. 

  24. <a name="line14"> 14: </a><font color="#B22222">             D = diagonal of A</font>
    25. 

  25. <a name="line15"> 15: </a><font color="#B22222">             P = permutation matrix determined by coloring</font>
    26. 

  26. <a name="line16"> 16: </a><font color="#B22222">   In this case, we use pre-solve and post-solve phases to handle scaling and</font>
    27. 

  27. <a name="line17"> 17: </a><font color="#B22222">   permutation, and by default the scaled residual norm is monitored for the</font>
    28. 

  28. <a name="line18"> 18: </a><font color="#B22222">   ILU/ICC preconditioners.  Use the option</font>
    29. 

  29. <a name="line19"> 19: </a><font color="#B22222">     -ksp_truemonitor</font>
    30. 

  30. <a name="line20"> 20: </a><font color="#B22222">   to print both the scaled and unscaled residual norms.</font>
    31. 

  31. <a name="line21"> 21: </a><font color="#B22222">*/</font>
    32. 

  32. 

  33. <a name="line23"> 23: </a> #include <A href="../../../../../include/petsc.h.html">petsc.h</A>
    34. 

  34. 

  35. <a name="line25"> 25: </a> #include <A href="../../../../../src/mat/impls/rowbs/mpi/mpirowbs.h.html">src/mat/impls/rowbs/mpi/mpirowbs.h</A>
    36. 

  36. 

  37. <a name="line29"> 29: </a><strong><font color="#4169E1"><a name="MatScaleSystem_MPIRowbs"></a><A href="../../../../../docs/manualpages/Sys/PetscErrorCode.html#PetscErrorCode">PetscErrorCode</A> MatScaleSystem_MPIRowbs(<A href="../../../../../docs/manualpages/Mat/Mat.html#Mat">Mat</A> mat,<A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A> rhs,<A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A> x)</font></strong>
    38. 

  38. <a name="line30"> 30: </a>{
    39. 

  39. <a name="line31"> 31: </a>  Mat_MPIRowbs *bsif  = (Mat_MPIRowbs*)mat-&gt;data;
    40. 

  40. <a name="line32"> 32: </a>  <A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A>          v = bsif-&gt;xwork;
    41. 

  41. <a name="line33"> 33: </a>  <A href="../../../../../docs/manualpages/Sys/PetscScalar.html#PetscScalar">PetscScalar</A>  *xa,*rhsa,*va;
    42. 

  42. 

  43. <a name="line37"> 37: </a>  <font color="#B22222">/* Permute and scale RHS and solution vectors */</font>
    44. 

  44. <a name="line38"> 38: </a>  <font color="#4169E1">if</font> (x) {
    45. 

  45. <a name="line39"> 39: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(x,&amp;xa);
    46. 

  46. <a name="line40"> 40: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(v,&amp;va);
    47. 

  47. <a name="line41"> 41: </a>    BSperm_dvec(xa,va,bsif-&gt;pA-&gt;perm);CHKERRBS(0);
    48. 

  48. <a name="line42"> 42: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(x,&amp;xa);
    49. 

  49. <a name="line43"> 43: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(v,&amp;va);
    50. 

  50. <a name="line44"> 44: </a>    <A href="../../../../../docs/manualpages/Vec/VecPointwiseDivide.html#VecPointwiseDivide">VecPointwiseDivide</A>(x,v,bsif-&gt;diag);
    51. 

  51. <a name="line45"> 45: </a>  }
    52. 

  52. 

  53. <a name="line47"> 47: </a>  <font color="#4169E1">if</font> (rhs) {
    54. 

  54. <a name="line48"> 48: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(rhs,&amp;rhsa);
    55. 

  55. <a name="line49"> 49: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(v,&amp;va);
    56. 

  56. <a name="line50"> 50: </a>    BSperm_dvec(rhsa,va,bsif-&gt;pA-&gt;perm);CHKERRBS(0);
    57. 

  57. <a name="line51"> 51: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(rhs,&amp;rhsa);
    58. 

  58. <a name="line52"> 52: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(v,&amp;va);
    59. 

  59. <a name="line53"> 53: </a>    <A href="../../../../../docs/manualpages/Vec/VecPointwiseMult.html#VecPointwiseMult">VecPointwiseMult</A>(rhs,v,bsif-&gt;diag);
    60. 

  60. <a name="line54"> 54: </a>  }
    61. 

  61. <a name="line55"> 55: </a>  <font color="#4169E1">return</font>(0);
    62. 

  62. <a name="line56"> 56: </a>}
    63. 

  63. 

  64. <a name="line60"> 60: </a><strong><font color="#4169E1"><a name="MatUnScaleSystem_MPIRowbs"></a><A href="../../../../../docs/manualpages/Sys/PetscErrorCode.html#PetscErrorCode">PetscErrorCode</A> MatUnScaleSystem_MPIRowbs(<A href="../../../../../docs/manualpages/Mat/Mat.html#Mat">Mat</A> mat,<A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A> rhs,<A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A> x)</font></strong>
    65. 

  65. <a name="line61"> 61: </a>{
    66. 

  66. <a name="line62"> 62: </a>  Mat_MPIRowbs *bsif  = (Mat_MPIRowbs*)mat-&gt;data;
    67. 

  67. <a name="line63"> 63: </a>  <A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A>          v = bsif-&gt;xwork;
    68. 

  68. <a name="line64"> 64: </a>  <A href="../../../../../docs/manualpages/Sys/PetscScalar.html#PetscScalar">PetscScalar</A>  *xa,*va,*rhsa;
    69. 

  69. 

  70. <a name="line68"> 68: </a>  <font color="#B22222">/* Unpermute and unscale the solution and RHS vectors */</font>
    71. 

  71. <a name="line69"> 69: </a>  <font color="#4169E1">if</font> (x) {
    72. 

  72. <a name="line70"> 70: </a>    <A href="../../../../../docs/manualpages/Vec/VecPointwiseMult.html#VecPointwiseMult">VecPointwiseMult</A>(v,x,bsif-&gt;diag);
    73. 

  73. <a name="line71"> 71: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(v,&amp;va);
    74. 

  74. <a name="line72"> 72: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(x,&amp;xa);
    75. 

  75. <a name="line73"> 73: </a>    BSiperm_dvec(va,xa,bsif-&gt;pA-&gt;perm);CHKERRBS(0);
    76. 

  76. <a name="line74"> 74: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(x,&amp;xa);
    77. 

  77. <a name="line75"> 75: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(v,&amp;va);
    78. 

  78. <a name="line76"> 76: </a>  }
    79. 

  79. <a name="line77"> 77: </a>  <font color="#4169E1">if</font> (rhs) {
    80. 

  80. <a name="line78"> 78: </a>    <A href="../../../../../docs/manualpages/Vec/VecPointwiseDivide.html#VecPointwiseDivide">VecPointwiseDivide</A>(v,rhs,bsif-&gt;diag);
    81. 

  81. <a name="line79"> 79: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(rhs,&amp;rhsa);
    82. 

  82. <a name="line80"> 80: </a>    <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(v,&amp;va);
    83. 

  83. <a name="line81"> 81: </a>    BSiperm_dvec(va,rhsa,bsif-&gt;pA-&gt;perm);CHKERRBS(0);
    84. 

  84. <a name="line82"> 82: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(rhs,&amp;rhsa);
    85. 

  85. <a name="line83"> 83: </a>    <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(v,&amp;va);
    86. 

  86. <a name="line84"> 84: </a>  }
    87. 

  87. <a name="line85"> 85: </a>  <font color="#4169E1">return</font>(0);
    88. 

  88. <a name="line86"> 86: </a>}
    89. 

  89. 

  90. <a name="line90"> 90: </a><strong><font color="#4169E1"><a name="MatUseScaledForm_MPIRowbs"></a><A href="../../../../../docs/manualpages/Sys/PetscErrorCode.html#PetscErrorCode">PetscErrorCode</A> MatUseScaledForm_MPIRowbs(<A href="../../../../../docs/manualpages/Mat/Mat.html#Mat">Mat</A> mat,<A href="../../../../../docs/manualpages/Sys/PetscTruth.html#PetscTruth">PetscTruth</A> scale)</font></strong>
    91. 

  91. <a name="line91"> 91: </a>{
    92. 

  92. <a name="line92"> 92: </a>  Mat_MPIRowbs *bsif  = (Mat_MPIRowbs*)mat-&gt;data;
    93. 

  93. 

  94. <a name="line95"> 95: </a>  bsif-&gt;vecs_permscale = scale;
    95. 

  95. <a name="line96"> 96: </a>  <font color="#4169E1">return</font>(0);
    96. 

  96. <a name="line97"> 97: </a>}
    97. 

  97. 

  98. </pre>
    99. 

  99. </body>
    100. 

  100. 

  101. </html>
    102. 

  102.