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

/hpc_kernel_samples/sparse_linear_algebra/casestudies/pagerank-petsc/petsc-2.3.2-p10/src/dm/da/examples/tutorials/ex5.c.html

https://gitlab.com/steinret/CodeVault
HTML | 77 lines | 64 code | 13 blank | 0 comment | 0 complexity | 799b35100787d1caf1eb38386b140ab4 MD5 | raw file
    

  1. <center><a href="ex5.c">Actual source code: ex5.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:33:34+00:00">
    8. 

  8. </head>
    9. 

  9. 

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

  11. <pre width="80">
    12. 

  12. <a name="line2">  2: </a>static char help[] = <font color="#666666">"Tests <A href="../../../../../docs/manualpages/DA/DAGetElements.html#DAGetElements">DAGetElements</A>() and <A href="../../../../../docs/manualpages/Vec/VecView.html#VecView">VecView</A>() contour plotting for 2d DAs.\n\n"</font>;
    13. 

  13. 

  14. <a name="line4"> 4: </a> #include <A href="../../../../../include/petscda.h.html">petscda.h</A>
    15. 

  15. <a name="line5"> 5: </a> #include <A href="../../../../../include/petscsys.h.html">petscsys.h</A>
    16. 

  16. 

  17. <a name="line9">  9: </a><strong><font color="#4169E1"><a name="main"></a>int main(int argc,char **argv)</font></strong>
    18. 

  18. <a name="line10"> 10: </a>{
    19. 

  19. <a name="line11"> 11: </a>  <A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A>       M = 10,N = 8,m = <A href="../../../../../docs/manualpages/Sys/PETSC_DECIDE.html#PETSC_DECIDE">PETSC_DECIDE</A>,n = <A href="../../../../../docs/manualpages/Sys/PETSC_DECIDE.html#PETSC_DECIDE">PETSC_DECIDE</A>,ne,i;
    20. 

  20. <a name="line12"> 12: </a>  const <A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A> *e;
    21. 

  21. <a name="line14"> 14: </a>  <A href="../../../../../docs/manualpages/Sys/PetscTruth.html#PetscTruth">PetscTruth</A>     flg;
    22. 

  22. <a name="line15"> 15: </a>  <A href="../../../../../docs/manualpages/DA/DA.html#DA">DA</A>             da;
    23. 

  23. <a name="line16"> 16: </a>  <A href="../../../../../docs/manualpages/Viewer/PetscViewer.html#PetscViewer">PetscViewer</A>    viewer;
    24. 

  24. <a name="line17"> 17: </a>  <A href="../../../../../docs/manualpages/Vec/Vec.html#Vec">Vec</A>            local,global;
    25. 

  25. <a name="line18"> 18: </a>  <A href="../../../../../docs/manualpages/Sys/PetscScalar.html#PetscScalar">PetscScalar</A>    value;
    26. 

  26. <a name="line19"> 19: </a>  <A href="../../../../../docs/manualpages/DA/DAPeriodicType.html#DAPeriodicType">DAPeriodicType</A> ptype = DA_NONPERIODIC;
    27. 

  27. <a name="line20"> 20: </a>  <A href="../../../../../docs/manualpages/DA/DAStencilType.html#DAStencilType">DAStencilType</A>  stype = DA_STENCIL_BOX;
    28. 

  28. <a name="line21"> 21: </a>  <A href="../../../../../docs/manualpages/Sys/PetscScalar.html#PetscScalar">PetscScalar</A>    *lv;
    29. 

  29. 

  30. <a name="line23"> 23: </a>  <A href="../../../../../docs/manualpages/Sys/PetscInitialize.html#PetscInitialize">PetscInitialize</A>(&amp;argc,&amp;argv,(char*)0,help);
    31. 

  31. <a name="line24"> 24: </a>  <A href="../../../../../docs/manualpages/Viewer/PetscViewerDrawOpen.html#PetscViewerDrawOpen">PetscViewerDrawOpen</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,0,<font color="#666666">""</font>,300,0,300,300,&amp;viewer);
    32. 

  32. 

  33. <a name="line26"> 26: </a>  <font color="#B22222">/* Read options */</font>
    34. 

  34. <a name="line27"> 27: </a>  <A href="../../../../../docs/manualpages/Sys/PetscOptionsGetInt.html#PetscOptionsGetInt">PetscOptionsGetInt</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<font color="#666666">"-M"</font>,&amp;M,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>);
    35. 

  35. <a name="line28"> 28: </a>  <A href="../../../../../docs/manualpages/Sys/PetscOptionsGetInt.html#PetscOptionsGetInt">PetscOptionsGetInt</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<font color="#666666">"-N"</font>,&amp;N,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>);
    36. 

  36. <a name="line29"> 29: </a>  <A href="../../../../../docs/manualpages/Sys/PetscOptionsGetInt.html#PetscOptionsGetInt">PetscOptionsGetInt</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<font color="#666666">"-m"</font>,&amp;m,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>);
    37. 

  37. <a name="line30"> 30: </a>  <A href="../../../../../docs/manualpages/Sys/PetscOptionsGetInt.html#PetscOptionsGetInt">PetscOptionsGetInt</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<font color="#666666">"-n"</font>,&amp;n,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>);
    38. 

  38. <a name="line31"> 31: </a>  <A href="../../../../../docs/manualpages/Sys/PetscOptionsHasName.html#PetscOptionsHasName">PetscOptionsHasName</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<font color="#666666">"-star_stencil"</font>,&amp;flg);
    39. 

  39. <a name="line32"> 32: </a>  <font color="#4169E1">if</font> (flg) stype = <A href="../../../../../docs/manualpages/DA/DA_STENCIL_STAR.html#DA_STENCIL_STAR">DA_STENCIL_STAR</A>;
    40. 

  40. 

  41. <a name="line34"> 34: </a>  <font color="#B22222">/* Create distributed array and get vectors */</font>
    42. 

  42. <a name="line35"> 35: </a>  <A href="../../../../../docs/manualpages/DA/DACreate2d.html#DACreate2d">DACreate2d</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,ptype,stype,M,N,m,n,1,1,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,&amp;da);
    43. 

  43. <a name="line36"> 36: </a>  <A href="../../../../../docs/manualpages/DA/DACreateGlobalVector.html#DACreateGlobalVector">DACreateGlobalVector</A>(da,&amp;global);
    44. 

  44. <a name="line37"> 37: </a>  <A href="../../../../../docs/manualpages/DA/DACreateLocalVector.html#DACreateLocalVector">DACreateLocalVector</A>(da,&amp;local);
    45. 

  45. 

  46. <a name="line39"> 39: </a>  value = -3.0;
    47. 

  47. <a name="line40"> 40: </a>  <A href="../../../../../docs/manualpages/Vec/VecSet.html#VecSet">VecSet</A>(global,value);
    48. 

  48. <a name="line41"> 41: </a>  <A href="../../../../../docs/manualpages/DA/DAGlobalToLocalBegin.html#DAGlobalToLocalBegin">DAGlobalToLocalBegin</A>(da,global,<A href="../../../../../docs/manualpages/Sys/INSERT_VALUES.html#INSERT_VALUES">INSERT_VALUES</A>,local);
    49. 

  49. <a name="line42"> 42: </a>  <A href="../../../../../docs/manualpages/DA/DAGlobalToLocalEnd.html#DAGlobalToLocalEnd">DAGlobalToLocalEnd</A>(da,global,<A href="../../../../../docs/manualpages/Sys/INSERT_VALUES.html#INSERT_VALUES">INSERT_VALUES</A>,local);
    50. 

  50. 

  51. <a name="line44"> 44: </a>  <A href="../../../../../docs/manualpages/DA/DAGetElements.html#DAGetElements">DAGetElements</A>(da,&amp;ne,&amp;e);
    52. 

  52. <a name="line45"> 45: </a>  <A href="../../../../../docs/manualpages/Vec/VecGetArray.html#VecGetArray">VecGetArray</A>(local,&amp;lv);
    53. 

  53. <a name="line46"> 46: </a>  <font color="#4169E1">for</font> (i=0; i&lt;ne; i++) {
    54. 

  54. <a name="line47"> 47: </a>    <A href="../../../../../docs/manualpages/Sys/PetscPrintf.html#PetscPrintf">PetscPrintf</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,<font color="#666666">"i %D e[3*i] %D %D %D\n"</font>,i,e[3*i],e[3*i+1],e[3*i+2]);
    55. 

  55. <a name="line48"> 48: </a>    lv[e[3*i]] = i;
    56. 

  56. <a name="line49"> 49: </a>  }
    57. 

  57. <a name="line50"> 50: </a>  <A href="../../../../../docs/manualpages/Vec/VecRestoreArray.html#VecRestoreArray">VecRestoreArray</A>(local,&amp;lv);
    58. 

  58. <a name="line51"> 51: </a>  <A href="../../../../../docs/manualpages/DA/DARestoreElements.html#DARestoreElements">DARestoreElements</A>(da,&amp;ne,&amp;e);
    59. 

  59. 

  60. <a name="line53"> 53: </a>  <A href="../../../../../docs/manualpages/DA/DALocalToGlobal.html#DALocalToGlobal">DALocalToGlobal</A>(da,local,<A href="../../../../../docs/manualpages/Sys/ADD_VALUES.html#ADD_VALUES">ADD_VALUES</A>,global);
    61. 

  61. 

  62. <a name="line55"> 55: </a>  <A href="../../../../../docs/manualpages/DA/DAView.html#DAView">DAView</A>(da,viewer);
    63. 

  63. <a name="line56"> 56: </a>  <A href="../../../../../docs/manualpages/Vec/VecView.html#VecView">VecView</A>(global,viewer);
    64. 

  64. 

  65. <a name="line58"> 58: </a>  <font color="#B22222">/* Free memory */</font>
    66. 

  66. <a name="line59"> 59: </a>  <A href="../../../../../docs/manualpages/Viewer/PetscViewerDestroy.html#PetscViewerDestroy">PetscViewerDestroy</A>(viewer);
    67. 

  67. <a name="line60"> 60: </a>  <A href="../../../../../docs/manualpages/Vec/VecDestroy.html#VecDestroy">VecDestroy</A>(local);
    68. 

  68. <a name="line61"> 61: </a>  <A href="../../../../../docs/manualpages/Vec/VecDestroy.html#VecDestroy">VecDestroy</A>(global);
    69. 

  69. <a name="line62"> 62: </a>  <A href="../../../../../docs/manualpages/DA/DADestroy.html#DADestroy">DADestroy</A>(da);
    70. 

  70. <a name="line63"> 63: </a>  <A href="../../../../../docs/manualpages/Sys/PetscFinalize.html#PetscFinalize">PetscFinalize</A>();
    71. 

  71. <a name="line64"> 64: </a>  <font color="#4169E1">return</font> 0;
    72. 

  72. <a name="line65"> 65: </a>}
    73. 

  73. <a name="line66"> 66: </a>
    74. 

  74. </pre>
    75. 

  75. </body>
    76. 

  76. 

  77. </html>
    78. 

  78.