PageRenderTime 26ms CodeModel.GetById 26ms RepoModel.GetById 0ms app.codeStats 1ms

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

https://gitlab.com/steinret/CodeVault
HTML | 94 lines | 77 code | 17 blank | 0 comment | 0 complexity | f4ac80750a19e20bffb3cb29d8772395 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:32:48+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 AODataRemap(). \n\n"</font>;
    13. 

  13. 

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

  15. 

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

  17. <a name="line9">  9: </a>{
    18. 

  18. <a name="line10"> 10: </a>  <A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A>       n,nglobal,bs = 1,*keys,*data,i,start,*news;
    19. 

  19. <a name="line12"> 12: </a>  <A href="../../../../../docs/manualpages/Sys/PetscMPIInt.html#PetscMPIInt">PetscMPIInt</A>    <A href="../../../../../docs/manualpages/Sys/rank.html#rank">rank</A>,<A href="../../../../../docs/manualpages/Sys/size.html#size">size</A>;
    20. 

  20. <a name="line13"> 13: </a>  <A href="../../../../../docs/manualpages/AO/AOData.html#AOData">AOData</A>         aodata;
    21. 

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

  22. 

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

  24. <a name="line17"> 17: </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>);
    25. 

  25. 

  26. <a name="line19"> 19: </a>  <A href="http://www.mcs.anl.gov/mpi/www/www3/MPI_Comm_rank.html#MPI_Comm_rank">MPI_Comm_rank</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,&amp;<A href="../../../../../docs/manualpages/Sys/rank.html#rank">rank</A>); n = <A href="../../../../../docs/manualpages/Sys/rank.html#rank">rank</A> + 2;
    27. 

  27. <a name="line20"> 20: </a>  <A href="http://www.mcs.anl.gov/mpi/www/www3/MPI_Allreduce.html#MPI_Allreduce">MPI_Allreduce</A>(&amp;n,&amp;nglobal,1,MPI_INT,MPI_SUM,<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>);
    28. 

  28. <a name="line21"> 21: </a>  <A href="http://www.mcs.anl.gov/mpi/www/www3/MPI_Comm_size.html#MPI_Comm_size">MPI_Comm_size</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,&amp;<A href="../../../../../docs/manualpages/Sys/size.html#size">size</A>);
    29. 

  29. 

  30. <a name="line23"> 23: </a>  <font color="#B22222">/*</font>
    31. 

  31. <a name="line24"> 24: </a><font color="#B22222">       Create a database with one  key and one segment</font>
    32. 

  32. <a name="line25"> 25: </a><font color="#B22222">  */</font>
    33. 

  33. <a name="line26"> 26: </a>  <A href="../../../../../docs/manualpages/AO/AODataCreateBasic.html#AODataCreateBasic">AODataCreateBasic</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,&amp;aodata);
    34. 

  34. 

  35. <a name="line28"> 28: </a>  <font color="#B22222">/*</font>
    36. 

  36. <a name="line29"> 29: </a><font color="#B22222">       Put one segment in the key</font>
    37. 

  37. <a name="line30"> 30: </a><font color="#B22222">  */</font>
    38. 

  38. <a name="line31"> 31: </a>  <A href="../../../../../docs/manualpages/AO/AODataKeyAdd.html#AODataKeyAdd">AODataKeyAdd</A>(aodata,<font color="#666666">"key1"</font>,<A href="../../../../../docs/manualpages/Sys/PETSC_DECIDE.html#PETSC_DECIDE">PETSC_DECIDE</A>,nglobal);
    39. 

  39. 

  40. <a name="line33"> 33: </a>  <font color="#B22222">/* allocate space for the keys each processor will provide */</font>
    41. 

  41. <a name="line34"> 34: </a>  <A href="../../../../../docs/manualpages/Sys/PetscMalloc.html#PetscMalloc">PetscMalloc</A>(n*<font color="#4169E1">sizeof</font>(<A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A>),&amp;keys);
    42. 

  42. 

  43. <a name="line36"> 36: </a>  <font color="#B22222">/*</font>
    44. 

  44. <a name="line37"> 37: </a><font color="#B22222">     We assign the first set of keys (0 to 2) to processor 0, etc.</font>
    45. 

  45. <a name="line38"> 38: </a><font color="#B22222">     This computes the first local key on each processor</font>
    46. 

  46. <a name="line39"> 39: </a><font color="#B22222">  */</font>
    47. 

  47. <a name="line40"> 40: </a>  <A href="http://www.mcs.anl.gov/mpi/www/www3/MPI_Scan.html#MPI_Scan">MPI_Scan</A>(&amp;n,&amp;start,1,MPI_INT,MPI_SUM,<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>);
    48. 

  48. <a name="line41"> 41: </a>  start -= n;
    49. 

  49. 

  50. <a name="line43"> 43: </a>  <font color="#4169E1">for</font> (i=0; i&lt;n; i++) {
    51. 

  51. <a name="line44"> 44: </a>    keys[i]     = start + i;
    52. 

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

  53. 

  54. <a name="line47"> 47: </a>  <font color="#B22222">/* </font>
    55. 

  55. <a name="line48"> 48: </a><font color="#B22222">      Allocate data for the first key and first segment </font>
    56. 

  56. <a name="line49"> 49: </a><font color="#B22222">  */</font>
    57. 

  57. <a name="line50"> 50: </a>  <A href="../../../../../docs/manualpages/Sys/PetscMalloc.html#PetscMalloc">PetscMalloc</A>(bs*n*<font color="#4169E1">sizeof</font>(<A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A>),&amp;data);
    58. 

  58. <a name="line51"> 51: </a>  <font color="#4169E1">for</font> (i=0; i&lt;n; i++) {
    59. 

  59. <a name="line52"> 52: </a>    data[i]   = start + i + 1; <font color="#B22222">/* the data is the neighbor to the right */</font>
    60. 

  60. <a name="line53"> 53: </a>  }
    61. 

  61. <a name="line54"> 54: </a>  data[n-1] = 0; <font color="#B22222">/* make it periodic */</font>
    62. 

  62. <a name="line55"> 55: </a>  <A href="../../../../../docs/manualpages/AO/AODataSegmentAdd.html#AODataSegmentAdd">AODataSegmentAdd</A>(aodata,<font color="#666666">"key1"</font>,<font color="#666666">"key1"</font>,bs,n,keys,data,PETSC_INT);
    63. 

  63. <a name="line56"> 56: </a>  <A href="../../../../../docs/manualpages/Sys/PetscFree.html#PetscFree">PetscFree</A>(data);
    64. 

  64. <a name="line57"> 57: </a>  <A href="../../../../../docs/manualpages/Sys/PetscFree.html#PetscFree">PetscFree</A>(keys);
    65. 

  65. 

  66. <a name="line59"> 59: </a>  <font color="#B22222">/*</font>
    67. 

  67. <a name="line60"> 60: </a><font color="#B22222">        View the database</font>
    68. 

  68. <a name="line61"> 61: </a><font color="#B22222">  */</font>
    69. 

  69. <a name="line62"> 62: </a>  <A href="../../../../../docs/manualpages/AO/AODataView.html#AODataView">AODataView</A>(aodata,<A href="../../../../../docs/manualpages/Viewer/PETSC_VIEWER_STDOUT_WORLD.html#PETSC_VIEWER_STDOUT_WORLD">PETSC_VIEWER_STDOUT_WORLD</A>);
    70. 

  70. <a name="line63"> 63: </a>
    71. 

  71. <a name="line64"> 64: </a>  <font color="#B22222">/*</font>
    72. 

  72. <a name="line65"> 65: </a><font color="#B22222">         Remap the database so that i -&gt; nglobal - i - 1</font>
    73. 

  73. <a name="line66"> 66: </a><font color="#B22222">  */</font>
    74. 

  74. <a name="line67"> 67: </a>  <A href="../../../../../docs/manualpages/Sys/PetscMalloc.html#PetscMalloc">PetscMalloc</A>(n*<font color="#4169E1">sizeof</font>(<A href="../../../../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</A>),&amp;news);
    75. 

  75. <a name="line68"> 68: </a>  <font color="#4169E1">for</font> (i=0; i&lt;n; i++) {
    76. 

  76. <a name="line69"> 69: </a>    news[i] = nglobal - i - start - 1;
    77. 

  77. <a name="line70"> 70: </a>  }
    78. 

  78. <a name="line71"> 71: </a>  <A href="../../../../../docs/manualpages/AO/AOCreateBasic.html#AOCreateBasic">AOCreateBasic</A>(<A href="../../../../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</A>,n,news,<A href="../../../../../docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A>,&amp;ao);
    79. 

  79. <a name="line72"> 72: </a>  <A href="../../../../../docs/manualpages/Sys/PetscFree.html#PetscFree">PetscFree</A>(news);
    80. 

  80. <a name="line73"> 73: </a>  <A href="../../../../../docs/manualpages/AO/AODataKeyRemap.html#AODataKeyRemap">AODataKeyRemap</A>(aodata,<font color="#666666">"key1"</font>,ao);
    81. 

  81. <a name="line74"> 74: </a>  <A href="../../../../../docs/manualpages/AO/AODestroy.html#AODestroy">AODestroy</A>(ao);
    82. 

  82. <a name="line75"> 75: </a>  <A href="../../../../../docs/manualpages/AO/AODataView.html#AODataView">AODataView</A>(aodata,<A href="../../../../../docs/manualpages/Viewer/PETSC_VIEWER_STDOUT_WORLD.html#PETSC_VIEWER_STDOUT_WORLD">PETSC_VIEWER_STDOUT_WORLD</A>);
    83. 

  83. <a name="line76"> 76: </a>  <A href="../../../../../docs/manualpages/AO/AODataDestroy.html#AODataDestroy">AODataDestroy</A>(aodata);
    84. 

  84. 

  85. <a name="line78"> 78: </a>  <A href="../../../../../docs/manualpages/Sys/PetscFinalize.html#PetscFinalize">PetscFinalize</A>();
    86. 

  86. <a name="line79"> 79: </a>  <font color="#4169E1">return</font> 0;
    87. 

  87. <a name="line80"> 80: </a>}
    88. 

  88. <a name="line81"> 81: </a>
    89. 

  89. 

  90. 

  91. </pre>
    92. 

  92. </body>
    93. 

  93. 

  94. </html>
    95. 

  95.