# Contents of /trunk/esys2/finley/src/finleyC/Assemble_PDEMatrix_System2.c

Revision 82 - (show annotations)
Tue Oct 26 06:53:54 2004 UTC (15 years, 5 months ago) by jgs
File MIME type: text/plain
File size: 6852 byte(s)
```Initial revision

```
 1 /* \$Id\$ */ 2 3 /**************************************************************/ 4 5 /* assembles the system of numEq PDEs into the stiffness matrix S: */ 6 7 /* -div(A*grad u)-div(B*u)+C*grad u + D*u */ 8 9 /* -(A_{k,i,m,j} u_m,j)_i-(B_{k,i,m} u_m)_i+C_{k,m,j} u_m,j-D_{k,m} u_m */ 10 11 /* u has numComp components. */ 12 13 /* Shape of the coefficients: */ 14 15 /* A = numEqu x numDim x numComp x numDim */ 16 /* B = numDim x numEqu x numComp */ 17 /* C = numEqu x numDim x numComp */ 18 /* D = numEqu x numComp */ 19 20 21 /**************************************************************/ 22 23 /* Copyrights by ACcESS Australia, 2003 */ 24 /* author: gross@access.edu.au */ 25 /* Version: \$Id\$ */ 26 27 /**************************************************************/ 28 29 30 #include "Common.h" 31 #include "Assemble.h" 32 33 /**************************************************************/ 34 35 void Finley_Assemble_PDEMatrix_System2(int NS,int numDim,int numQuad,int numEqu,int numComp, 36 double* S,double* DSDX, double* Vol, 37 int NN, double* EM_S, 38 double* A, int extendedA, 39 double* B, int extendedB, 40 double* C, int extendedC, 41 double* D, int extendedD ) { 42 int s,r,k,m,i,j,q; 43 double rtmp; 44 45 46 /**************************************************************/ 47 /* process A: */ 48 /**************************************************************/ 49 if (NULL!=A) { 50 if (extendedA) { 51 for (s=0;s

## Properties

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