# Annotation of /trunk/paso/src/Solver_applyBlockDiagonalMatrix.c

Revision 3642 - (hide annotations)
Thu Oct 27 03:41:51 2011 UTC (7 years, 7 months ago) by caltinay
File MIME type: text/plain
File size: 2538 byte(s)
```Assorted spelling/comment fixes in paso.

```
 1 ksteube 1312 2 /******************************************************* 3 ksteube 1811 * 4 jfenwick 2881 * Copyright (c) 2003-2010 by University of Queensland 5 ksteube 1811 * Earth Systems Science Computational Center (ESSCC) 6 * http://www.uq.edu.au/esscc 7 * 8 * Primary Business: Queensland, Australia 9 * Licensed under the Open Software License version 3.0 10 * http://www.opensource.org/licenses/osl-3.0.php 11 * 12 *******************************************************/ 13 dhawcroft 631 14 ksteube 1811 15 jgs 150 /**************************************************************/ 16 17 caltinay 3642 /* Paso: apply block diagonal matrix D: x=D*b */ 18 jgs 150 19 caltinay 3642 /* should be called within a parallel region */ 20 /* barrier synchronization should be performed to make sure */ 21 /* that the input vector available */ 22 jgs 150 23 /**************************************************************/ 24 25 caltinay 3642 /* Copyrights by ACcESS Australia 2003, 2004, 2005 */ 26 /* Author: Lutz Gross, l.gross@uq.edu.au */ 27 jgs 150 28 /**************************************************************/ 29 30 gross 700 #include "Paso.h" 31 jgs 150 32 /**************************************************************/ 33 34 35 void Paso_Solver_applyBlockDiagonalMatrix(dim_t n_block,dim_t n,double* D,index_t* pivot,double* x,double* b) { 36 dim_t i; 37 gross 495 register dim_t i3,i9; 38 register double b0,b1,b2,D00,D10,D20,D01,D11,D21,D02,D12,D22; 39 40 jgs 150 if (n_block==1) { 41 gross 1556 #pragma omp parallel for private(i) schedule(static) 42 gross 495 for (i=0;i

## Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision