00001 #include "gsurf.h"
00002 #include "gnode.h"
00003 #include "gtopology.h"
00004 #include <math.h>
00005
00006 gsurf::gsurf (void)
00007 {
00008
00009 nn=0;
00010
00011 ndofn=NULL;
00012
00013
00014 re=-1;
00015
00016 adjel=NULL;
00017
00018
00019 nlistfn = NULL;
00020 nlistln = NULL;
00021
00022 cnmult = NULL;
00023
00024
00025 nv = NULL;
00026
00027
00028 threshold=1.0e-6;
00029 }
00030
00031
00032 gsurf::~gsurf (void)
00033 {
00034 long i;
00035
00036 delete [] ndofn;
00037 delete [] adjel;
00038 delete [] nlistfn;
00039 delete [] nlistln;
00040 delete [] nv;
00041
00042 for (i=0;i<nn;i++){
00043 delete [] cnmult[i];
00044 }
00045 delete [] cnmult;
00046 }
00047
00048
00049
00050
00051
00052
00053
00054
00055 void gsurf::normal_vector (gnode *gnodes)
00056 {
00057 double norm;
00058
00059 nv = new double [3];
00060
00061
00062 norm = nv[0]*nv[0] + nv[1]*nv[1] + nv[2]*nv[2];
00063
00064 if (fabs(norm)<threshold){
00065 print_err("zero length of the normal vector", __FILE__, __LINE__, __func__);
00066 }
00067
00068 norm = sqrt(norm);
00069
00070 nv[0]/=norm;
00071 nv[1]/=norm;
00072 nv[2]/=norm;
00073 }
00074
00075
00076
00077
00078
00079
00080
00081
00082 void gsurf::print (FILE *out)
00083 {
00084 long i;
00085
00086 fprintf (out,"\n");
00087 fprintf (out,"\n number of nodes %ld",nn);
00088 fprintf (out,"\n node multiplicity %ld",nm);
00089
00090 if (nlistfn!=NULL){
00091 fprintf (out,"\n surface nodes\n");
00092 for (i=0;i<nn;i++){
00093 fprintf (out," %ld",nlistfn[i]);
00094 }
00095 fprintf (out,"\n");
00096 for (i=0;i<nn;i++){
00097 fprintf (out," %ld",nlistln[i]);
00098 }
00099 }
00100 if (nv!=NULL){
00101 fprintf (out,"\n normal vector");
00102 fprintf (out,"\n %lf %lf %lf",nv[0],nv[1],nv[2]);
00103 }
00104
00105 }
00106
00107
00108
00109
00110
00111
00112 void gsurf::give_norvect (double *v)
00113 {
00114 v[0]=nv[0];
00115 v[1]=nv[1];
00116 v[2]=nv[2];
00117 }
00118
00119
00120
00121
00122
00123
00124 void gsurf::check_normal (vector &x,vector &y,ivector &nod)
00125 {
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167 }
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177 void gsurf::alloc_cn (long nccnfn,long nccnln)
00178 {
00179 long i,ndofnf,ndofnl;
00180
00181
00182 ndofnf=nccnfn;
00183
00184 ndofnl=nccnln;
00185
00186 if (ndofnf!=ndofnl){
00187
00188 }
00189 else{}
00190
00191 if (cnmult!=NULL){
00192 for (i=0;i<nn;i++){
00193 delete [] cnmult[i];
00194 }
00195 delete [] cnmult;
00196 }
00197
00198 cnmult = new long* [nn];
00199 for (i=0;i<nn;i++){
00200 cnmult[i] = new long [ndofnf];
00201 }
00202 }
00203
00204
00205
00206
00207
00208
00209
00210
00211 void gsurf::give_first_node_numbers (long *fnn)
00212 {
00213 long i;
00214
00215 for (i=0;i<nn;i++){
00216 fnn[i]=nlistfn[i];
00217 }
00218 }
00219
00220
00221
00222
00223
00224
00225
00226
00227 void gsurf::give_last_node_numbers (long *lnn)
00228 {
00229 long i;
00230
00231 for (i=0;i<nn;i++){
00232 lnn[i]=nlistln[i+nn];
00233 }
00234 }
00235
00236
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248 void gsurf::give_mult_code_numbers (long fln,long *mcn)
00249 {
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264 }
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275
00276
00277
00278
00279
00280
00281
00282
00283
00284