/[escript]/trunk/paso/src/SystemMatrixPattern_unrollBlocks.c
ViewVC logotype

Contents of /trunk/paso/src/SystemMatrixPattern_unrollBlocks.c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2608 - (show annotations)
Tue Aug 18 01:25:18 2009 UTC (10 years, 1 month ago) by jfenwick
File MIME type: text/plain
File size: 4081 byte(s)
Updating Lutz' email

1
2 /*******************************************************
3 *
4 * Copyright (c) 2003-2009 by University of Queensland
5 * 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
14
15 /**************************************************************/
16
17 /* Paso: SystemMatrixPattern_unrollBlocks */
18
19 /**************************************************************/
20
21 /* Author: Lutz Gross, l.gross@uq.edu.au */
22
23 /**************************************************************/
24
25 #include "SystemMatrixPattern.h"
26
27 /**************************************************************/
28
29 /* creates SystemMatrixPattern */
30
31
32 Paso_SystemMatrixPattern* Paso_SystemMatrixPattern_unrollBlocks(Paso_SystemMatrixPattern* pattern,
33 int type, dim_t output_block_size,dim_t input_block_size) {
34 Paso_SystemMatrixPattern*out=NULL;
35 Paso_Pattern *new_mainPattern=NULL, *new_col_couplePattern=NULL, *new_row_couplePattern=NULL;
36 Paso_Distribution* new_output_distribution=NULL, *new_input_distribution=NULL;
37 Paso_Connector *new_col_connector=NULL, *new_row_connector=NULL;
38
39 if ( ( output_block_size == 1 ) && (input_block_size == 1) && ((pattern->type & PATTERN_FORMAT_OFFSET1) == (type & PATTERN_FORMAT_OFFSET1) ) ) {
40 out = Paso_SystemMatrixPattern_getReference(pattern);
41 } else {
42 new_mainPattern=Paso_Pattern_unrollBlocks(pattern->mainPattern,type,output_block_size,input_block_size);
43 new_col_couplePattern=Paso_Pattern_unrollBlocks(pattern->col_couplePattern,type,output_block_size,input_block_size);
44 new_row_couplePattern=Paso_Pattern_unrollBlocks(pattern->row_couplePattern,type,output_block_size,input_block_size);
45 if (output_block_size>1) {
46 new_output_distribution=Paso_Distribution_alloc(pattern->output_distribution->mpi_info,
47 pattern->output_distribution->first_component,
48 output_block_size,0);
49 new_row_connector=Paso_Connector_unroll(pattern->row_connector,output_block_size);
50 } else {
51 new_output_distribution=Paso_Distribution_getReference(pattern->output_distribution);
52 new_row_connector= Paso_Connector_getReference(pattern->row_connector);
53 }
54 if (input_block_size>1) {
55 new_input_distribution=Paso_Distribution_alloc(pattern->input_distribution->mpi_info,
56 pattern->input_distribution->first_component,
57 input_block_size,0);
58 new_col_connector=Paso_Connector_unroll(pattern->col_connector,input_block_size);
59 } else {
60 new_input_distribution=Paso_Distribution_getReference(pattern->input_distribution);
61 new_col_connector=Paso_Connector_getReference(pattern->col_connector);
62 }
63 if (Paso_noError()) {
64 out=Paso_SystemMatrixPattern_alloc(type,
65 new_output_distribution,
66 new_input_distribution,
67 new_mainPattern,
68 new_col_couplePattern,
69 new_row_couplePattern,
70 new_col_connector,
71 new_row_connector);
72 }
73 Paso_Pattern_free(new_mainPattern);
74 Paso_Pattern_free(new_col_couplePattern);
75 Paso_Pattern_free(new_row_couplePattern);
76 Paso_Distribution_free(new_output_distribution);
77 Paso_Distribution_free(new_input_distribution);
78 Paso_Connector_free(new_row_connector);
79 Paso_Connector_free(new_col_connector);
80 }
81 if (Paso_noError()) {
82 return out;
83 } else {
84 Paso_SystemMatrixPattern_free(out);
85 return NULL;
86 }
87 }

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26