00001 #include <string.h> 00002 #include "lhsrhsc.h" 00003 #include "globalc.h" 00004 #include "gtopology.h" 00005 00006 00007 00008 lhsrhsc::lhsrhsc () 00009 { 00010 nlc=Cb->nlc; lhs=NULL; lhsi=NULL; tdlhs=NULL; rhs=NULL; 00011 } 00012 00013 lhsrhsc::~lhsrhsc () 00014 { 00015 delete [] lhs; delete [] lhsi; delete [] tdlhs; delete [] rhs; 00016 } 00017 00018 /** 00019 function allocates array for left hand side, time derivative of 00020 left hand side, initial values and right sides of problem 00021 00022 6.4.2003, JK 00023 */ 00024 void lhsrhsc::alloc () 00025 { 00026 lhs = new double [nlc*Ndofc]; 00027 memset (lhs,0,(nlc*Ndofc)*sizeof(double)); 00028 lhsi = new double [nlc*Ndofc]; 00029 memset (lhsi,0,(nlc*Ndofc)*sizeof(double)); 00030 tdlhs = new double [nlc*Ndofc]; 00031 memset (tdlhs,0,(nlc*Ndofc)*sizeof(double)); 00032 rhs = new double [nlc*Ndofc]; 00033 memset (rhs,0,(nlc*Ndofc)*sizeof(double)); 00034 } 00035 00036 double *lhsrhsc::give_lhs (long lcid) 00037 { 00038 return (lhs+lcid*Ndofc); 00039 } 00040 00041 double *lhsrhsc::give_lhsi (long lcid) 00042 { 00043 return (lhsi+lcid*Ndofc); 00044 } 00045 00046 double *lhsrhsc::give_tdlhs (long lcid) 00047 { 00048 return (tdlhs+lcid*Ndofc); 00049 } 00050 00051 double *lhsrhsc::give_rhs (long lcid) 00052 { 00053 return (rhs+lcid*Ndofc); 00054 } 00055 00056 void lhsrhsc::output (FILE *out,long lcid) 00057 { 00058 } 00059 00060 /** 00061 function reads initial conditions 00062 00063 @param in - input data stream 00064 17.1.2002 00065 */ 00066 void lhsrhsc::initcond (FILE *in) 00067 { 00068 00069 } 00070