/[escript]/branches/domexper/dudley/src/Assemble_getAssembleParameters.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Assemble_getAssembleParameters.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3205 by jfenwick, Fri Sep 24 00:30:43 2010 UTC revision 3206 by jfenwick, Fri Sep 24 03:20:22 2010 UTC
# Line 74  void Assemble_getAssembleParameters(Dudl Line 74  void Assemble_getAssembleParameters(Dudl
74        if ( Paso_Distribution_getMyNumComponents(S->row_distribution)*S->row_block_size==parm->numEqu* Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution)) {        if ( Paso_Distribution_getMyNumComponents(S->row_distribution)*S->row_block_size==parm->numEqu* Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution)) {
75             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);
76             parm->row_DOF=nodes->degreesOfFreedomMapping->target;             parm->row_DOF=nodes->degreesOfFreedomMapping->target;
77             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,FALSE,reducedIntegrationOrder);             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
78        }        }
79        else if ( Paso_Distribution_getMyNumComponents(S->row_distribution)*S->row_block_size==parm->numEqu* Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution)) {        else if ( Paso_Distribution_getMyNumComponents(S->row_distribution)*S->row_block_size==parm->numEqu* Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution)) {
80             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);
81             parm->row_DOF=nodes->reducedDegreesOfFreedomMapping->target;             parm->row_DOF=nodes->reducedDegreesOfFreedomMapping->target;
82         parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,TRUE,reducedIntegrationOrder);         parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
83        } else {        } else {
84             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: number of rows in matrix does not match the number of degrees of freedom in mesh");             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: number of rows in matrix does not match the number of degrees of freedom in mesh");
85        }        }
# Line 87  void Assemble_getAssembleParameters(Dudl Line 87  void Assemble_getAssembleParameters(Dudl
87        if ( Paso_Distribution_getMyNumComponents(S->col_distribution)*S->col_block_size==parm->numComp* Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution)) {        if ( Paso_Distribution_getMyNumComponents(S->col_distribution)*S->col_block_size==parm->numComp* Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution)) {
88             parm->col_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);             parm->col_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);
89             parm->col_DOF=nodes->degreesOfFreedomMapping->target;             parm->col_DOF=nodes->degreesOfFreedomMapping->target;
90             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,FALSE,reducedIntegrationOrder);             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
91        } else if ( Paso_Distribution_getMyNumComponents(S->col_distribution)*S->col_block_size==parm->numComp* Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution)) {        } else if ( Paso_Distribution_getMyNumComponents(S->col_distribution)*S->col_block_size==parm->numComp* Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution)) {
92             parm->col_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);             parm->col_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);
93             parm->col_DOF=nodes->reducedDegreesOfFreedomMapping->target;             parm->col_DOF=nodes->reducedDegreesOfFreedomMapping->target;
94             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,TRUE,reducedIntegrationOrder);             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
95        } else {        } else {
96             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: number of columns in matrix does not match the number of degrees of freedom in mesh");             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: number of columns in matrix does not match the number of degrees of freedom in mesh");
97        }        }
# Line 102  void Assemble_getAssembleParameters(Dudl Line 102  void Assemble_getAssembleParameters(Dudl
102        if (numSamplesEqual(F,1, Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution))) {        if (numSamplesEqual(F,1, Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution))) {
103             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->degreesOfFreedomDistribution);
104             parm->row_DOF=nodes->degreesOfFreedomMapping->target;             parm->row_DOF=nodes->degreesOfFreedomMapping->target;
105             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,FALSE,reducedIntegrationOrder);             parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
106        }        }
107        else if (numSamplesEqual(F,1, Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution))) {        else if (numSamplesEqual(F,1, Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution))) {
108             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);             parm->row_DOF_UpperBound =  Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);
109             parm->row_DOF=nodes->reducedDegreesOfFreedomMapping->target;             parm->row_DOF=nodes->reducedDegreesOfFreedomMapping->target;
110         parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,TRUE,reducedIntegrationOrder);         parm->row_jac=Dudley_ElementFile_borrowJacobeans(elements,nodes,reducedIntegrationOrder);
111        } else {        } else {
112             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: length of RHS vector does not match the number of degrees of freedom in mesh");             Dudley_setError(TYPE_ERROR,"Assemble_getAssembleParameters: length of RHS vector does not match the number of degrees of freedom in mesh");
113        }        }

Legend:
Removed from v.3205  
changed lines
  Added in v.3206

  ViewVC Help
Powered by ViewVC 1.1.26