| version 1.133, 2006/06/05 08:11:10 |
version 1.134, 2006/06/06 07:14:16 |
|
|
| /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.132 2006/06/05 01:29:24 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.133 2006/06/05 08:11:10 noro Exp $ */ |
| |
|
| #include "nd.h" |
#include "nd.h" |
| |
|
| Line 4185 int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr |
|
| Line 4185 int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr |
|
| |
|
| v.id = O_VECT; v.len = col; v.body = (pointer *)svect; |
v.id = O_VECT; v.len = col; v.body = (pointer *)svect; |
| maxrs = 0; |
maxrs = 0; |
| hmag = p_mag((P)svect[0])*nd_scale; |
for ( i = 0; i < col && !svect[i]; i++ ); |
| |
if ( i == col ) return maxrs; |
| |
hmag = p_mag((P)svect[i])*nd_scale; |
| for ( i = 0; i < nred; i++ ) { |
for ( i = 0; i < nred; i++ ) { |
| ivect = imat[i]; |
ivect = imat[i]; |
| k = ivect->head; |
k = ivect->head; |
| Line 4224 int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr |
|
| Line 4226 int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr |
|
| } |
} |
| break; |
break; |
| } |
} |
| |
for ( j = k+1; j < col && !svect[j]; j++ ); |
| |
if ( j == col ) break; |
| |
if ( hmag && ((double)p_mag((P)svect[j]) > hmag) ) { |
| |
igcdv(&v,&t); |
| |
hmag = ((double)p_mag((P)svect[j]))*nd_scale; |
| |
} |
| } |
} |
| if ( hmag && ((double)p_mag((P)svect[0]) > hmag) ) |
|
| igcdv(&v,&t); |
|
| } |
} |
| igcdv(&v,&t); |
igcdv(&v,&t); |
| if ( DP_Print ) { |
if ( DP_Print ) { |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to nd.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / engine