static DaoArray* DaoArray_CopyX( DaoArray *self, DaoType *tp )
{
DaoArray *copy = DaoArray_New( self->etype );
if( tp && tp->tid == DAO_ARRAY && tp->nested->size ){
int nt = tp->nested->items.pType[0]->tid;
if( nt >= DAO_INTEGER && nt <= DAO_COMPLEX ) copy->etype = nt;
}
DaoArray_ResizeArray( copy, self->dims, self->ndim );
DaoArray_CopyArray( copy, self );
return copy;
}
static void STD_Array( DaoProcess *proc, DaoValue *p[], int N )
{
DaoInteger idint = {DAO_INTEGER,0,0,0,0,0};
DaoValue *res, *index = (DaoValue*)(void*)&idint;
DaoVmCode *sect = DaoGetSectionCode( proc->activeCode );
DaoArray *array = DaoProcess_PutArray( proc );
DaoArray *first = NULL;
DaoArray *sub = NULL;
daoint i, j, k, entry, size = 1;
/* if multi-dimensional array is disabled, DaoProcess_PutArray() will raise exception. */
#ifdef DAO_WITH_NUMARRAY
for(i=0; i<N; i++){
daoint d = p[i]->xInteger.value;
if( d < 0 ){
DaoProcess_RaiseException( proc, DAO_ERROR_PARAM, NULL );
break;
}
size *= d;
}
if( size == 0 ) return;
if( sect == NULL ) return; // TODO exception
if( DaoProcess_PushSectionFrame( proc ) == NULL ) return;
entry = proc->topFrame->entry;
DaoProcess_AcquireCV( proc );
for(i=0; i<size; i++){
idint.value = i;
if( sect->b >0 ) DaoProcess_SetValue( proc, sect->a, index );
proc->topFrame->entry = entry;
DaoProcess_Execute( proc );
if( proc->status == DAO_PROCESS_ABORTED ) break;
res = proc->stackValues[0];
if( i == 0 ){
int D = N;
DaoArray_SetDimCount( array, N + (res->type == DAO_ARRAY ? res->xArray.ndim : 0) );
for(j=0; j<N; j++) array->dims[j] = p[j]->xInteger.value;
if( res->type == DAO_ARRAY ){
first = DaoArray_Copy( (DaoArray*) res );
if( first->ndim == 2 && (first->dims[0] == 1 || first->dims[1] == 1) ){
D += 1;
array->dims[N] = first->dims[ first->dims[0] == 1 ];
}else{
D += first->ndim;
memmove( array->dims + N, first->dims, first->ndim*sizeof(daoint) );
}
}
DaoArray_ResizeArray( array, array->dims, D );
}
if( res->type == DAO_ARRAY ){
sub = (DaoArray*) res;
if( first == NULL || DaoArray_AlignShape( sub, NULL, first->dims, first->ndim ) ==0 ){
DaoProcess_RaiseException( proc, DAO_ERROR, "inconsistent elements or subarrays" );
break;
}
k = i * sub->size;
for(j=0; j<sub->size; j++){
switch( array->etype ){
case DAO_INTEGER : array->data.i[k+j] = DaoArray_GetInteger( sub, j ); break;
case DAO_FLOAT : array->data.f[k+j] = DaoArray_GetFloat( sub, j ); break;
case DAO_DOUBLE : array->data.d[k+j] = DaoArray_GetDouble( sub, j ); break;
case DAO_COMPLEX : array->data.c[k+j] = DaoArray_GetComplex( sub, j ); break;
}
}
}else{
switch( array->etype ){
case DAO_INTEGER : array->data.i[i] = DaoValue_GetInteger( res ); break;
case DAO_FLOAT : array->data.f[i] = DaoValue_GetFloat( res ); break;
case DAO_DOUBLE : array->data.d[i] = DaoValue_GetDouble( res ); break;
case DAO_COMPLEX : array->data.c[i] = DaoValue_GetComplex( res ); break;
}
}
}
DaoProcess_ReleaseCV( proc );
DaoProcess_PopFrame( proc );
if( first ) DaoArray_Delete( first );
#endif
}
static void DaoMT_Functional( DaoProcess *proc, DaoValue *P[], int N, int F )
{
DMutex mutex;
DCondVar condv;
DaoTaskData *tasks;
DaoValue *param = P[0];
DaoValue *result = NULL;
DaoList *list = NULL;
DaoArray *array = NULL;
DaoVmCode *sect = DaoGetSectionCode( proc->activeCode );
int i, entry, threads = P[1]->xInteger.value;
daoint index = -1, status = 0, joined = 0;
DNode *node = NULL;
switch( F ){
case DVM_FUNCT_MAP :
if( param->type == DAO_ARRAY ){
array = DaoProcess_PutArray( proc );
result = (DaoValue*) array;
}else{
list = DaoProcess_PutList( proc );
result = (DaoValue*) list;
}
break;
case DVM_FUNCT_APPLY : DaoProcess_PutValue( proc, param ); break;
case DVM_FUNCT_FIND : DaoProcess_PutValue( proc, dao_none_value ); break;
}
if( threads <= 0 ) threads = 2;
if( sect == NULL || DaoMT_PushSectionFrame( proc ) == 0 ) return;
if( list ){
DArray_Clear( & list->items );
if( param->type == DAO_LIST ) DArray_Resize( & list->items, param->xList.items.size, NULL );
if( param->type == DAO_MAP ) DArray_Resize( & list->items, param->xMap.items->size, NULL );
#ifdef DAO_WITH_NUMARRAY
}else if( array && F == DVM_FUNCT_MAP ){
DaoArray_GetSliceShape( (DaoArray*) param, & array->dims, & array->ndim );
DaoArray_ResizeArray( array, array->dims, array->ndim );
#endif
}
DMutex_Init( & mutex );
DCondVar_Init( & condv );
entry = proc->topFrame->entry;
tasks = (DaoTaskData*) dao_calloc( threads, sizeof(DaoTaskData) );
DaoProcess_PopFrame( proc );
for(i=0; i<threads; i++){
DaoTaskData *task = tasks + i;
task->param = param;
task->result = result;
task->proto = proc;
task->sect = sect;
task->funct = F;
task->entry = entry;
task->first = i;
task->step = threads;
task->index = & index;
task->node = & node;
task->joined = & joined;
task->condv = & condv;
task->mutex = & mutex;
task->clone = DaoVmSpace_AcquireProcess( proc->vmSpace );
task->clone->mutex = & mutex;
if( i ) DaoCallServer_AddTask( DaoMT_RunFunctional, task, 1 );
}
DaoMT_RunFunctional( tasks );
DMutex_Lock( & mutex );
while( joined < threads ) DCondVar_TimedWait( & condv, & mutex, 0.01 );
DMutex_Unlock( & mutex );
for(i=0; i<threads; i++){
DaoTaskData *task = tasks + i;
DaoVmSpace_ReleaseProcess( proc->vmSpace, task->clone );
status |= task->status;
}
if( F == DVM_FUNCT_FIND ){
DaoTuple *tuple = DaoProcess_PutTuple( proc, 0 );
if( param->type == DAO_LIST && index != -1 ){
DaoValue **items = param->xList.items.items.pValue;
GC_ShiftRC( items[index], tuple->items[1] );
tuple->items[1] = items[index];
tuple->items[0]->xInteger.value = index;
}else if( param->type == DAO_MAP && node ){
GC_ShiftRC( node->key.pValue, tuple->items[0] );
GC_ShiftRC( node->value.pValue, tuple->items[1] );
tuple->items[0] = node->key.pValue;
tuple->items[1] = node->value.pValue;
}
}
if( status ) DaoProcess_RaiseException( proc, DAO_ERROR, "code section execution failed!" );
DMutex_Destroy( & mutex );
DCondVar_Destroy( & condv );
dao_free( tasks );
}
/*
// Note: reference count is handled for "value2"!
//
// Item of list/tuple etc. can be directly passed as parameter "value2",
// to avoid creating unnecessary intermediate objects.
*/
static int DaoParser_Deserialize( DaoParser *self, int start, int end, DaoValue **value2, DArray *types, DaoNamespace *ns, DaoProcess *proc, DMap *omap )
{
DaoToken **tokens = self->tokens->items.pToken;
DaoType *it1 = NULL, *it2 = NULL, *type = NULL;
DaoValue *value = *value2;
DaoValue *tmp = NULL;
DaoValue *tmp2 = NULL;
DaoObject *object;
DaoCdata *cdata;
DaoArray *array;
DaoTuple *tuple;
DaoList *list;
DaoMap *map;
DArray *dims;
DNode *node;
void *key = NULL;
char *str;
int i, j, k, n;
int minus = 0;
int next = start + 1;
int tok2 = start < end ? tokens[start+1]->type : 0;
int maybetype = tok2 == DTOK_COLON2 || tok2 == DTOK_LT || tok2 == DTOK_LCB;
if( tokens[start]->type == DTOK_AT && tok2 == DTOK_LCB ){
int rb = DaoParser_FindPairToken( self, DTOK_LCB, DTOK_RCB, start, end );
if( rb < 0 ) return next;
sscanf( tokens[start+2]->string.mbs, "%p", & key );
node = DMap_Find( omap, key );
if( node ) DaoValue_Copy( node->value.pValue, value2 );
return rb + 1;
}
if( tokens[start]->name == DTOK_ID_SYMBOL ){
DString *mbs = DString_New(1);
while( tokens[start]->name == DTOK_ID_SYMBOL ){
DString_Append( mbs, & tokens[start]->string );
start += 1;
}
type = DaoNamespace_MakeType( ns, mbs->mbs, DAO_ENUM, NULL, NULL, 0 );
DString_Delete( mbs );
if( type == NULL ) return start;
if( tokens[start]->name != DTOK_LCB ) return start;
end = DaoParser_FindPairToken( self, DTOK_LCB, DTOK_RCB, start, end );
if( end < 0 ) return start;
next = end + 1;
start += 1;
end -= 1;
}else if( tokens[start]->type == DTOK_IDENTIFIER && maybetype ){
type = DaoParser_ParseType( self, start, end, & start, NULL );
if( type == NULL ) return next;
if( tokens[start]->name != DTOK_LCB ) return start;
end = DaoParser_FindPairToken( self, DTOK_LCB, DTOK_RCB, start, end );
if( end < 0 ) return start;
next = end + 1;
start += 1;
end -= 1;
}
if( type == NULL ){
type = types->items.pType[0];
if( type && type->tid >= DAO_ARRAY ){
if( tokens[start]->name != DTOK_LCB ) return start;
end = DaoParser_FindPairToken( self, DTOK_LCB, DTOK_RCB, start, end );
if( end < 0 ) return start;
next = end + 1;
start += 1;
end -= 1;
}
}
if( type == NULL ) return next;
DaoValue_Copy( type->value, value2 );
if( start > end ) return next;
if( tokens[start]->name == DTOK_SUB ){
minus = 1;
start += 1;
if( start > end ) return next;
}
if( type->nested && type->nested->size >0 ) it1 = type->nested->items.pType[0];
if( type->nested && type->nested->size >1 ) it2 = type->nested->items.pType[1];
if( tokens[start]->name == DTOK_LB ){
int rb = DaoParser_FindPairToken( self, DTOK_LB, DTOK_RB, start, end );
if( rb < 0 ) return next;
sscanf( tokens[start+1]->string.mbs, "%p", & key );
DMap_Insert( omap, key, *value2 );
start = rb + 1;
}
str = tokens[start]->string.mbs;
#if 0
printf( "type: %s %s\n", type->name->mbs, str );
for(i=start; i<=end; i++) printf( "%s ", tokens[i]->string.mbs ); printf( "\n" );
#endif
value = *value2;
switch( type->tid ){
case DAO_NONE :
break;
case DAO_INTEGER :
value->xInteger.value = DaoDecodeInteger( str );
if( minus ) value->xInteger.value = - value->xInteger.value;
break;
case DAO_FLOAT :
value->xFloat.value = DaoDecodeDouble( str );
if( minus ) value->xFloat.value = - value->xFloat.value;
break;
case DAO_DOUBLE :
value->xDouble.value = DaoDecodeDouble( str );
if( minus ) value->xDouble.value = - value->xDouble.value;
break;
case DAO_COMPLEX :
value->xComplex.value.real = DaoDecodeDouble( str );
if( minus ) value->xComplex.value.real = - value->xComplex.value.real;
if( start + 1 > end ) return start+1;
minus = 0;
if( tokens[start + 1]->name == DTOK_SUB ){
minus = 1;
start += 1;
if( start + 1 > end ) return start+1;
}
value->xComplex.value.imag = DaoDecodeDouble( tokens[start+1]->string.mbs );
if( minus ) value->xComplex.value.imag = - value->xComplex.value.imag;
next = start + 2;
break;
#ifdef DAO_WITH_LONGINT
case DAO_LONG :
value->xLong.value->base = DaoDecodeInteger( str );
start += 1;
if( tokens[start]->name == DTOK_ADD ){
value->xLong.value->sign = 1;
start += 1;
}else if( tokens[start]->name == DTOK_SUB ){
value->xLong.value->sign = -1;
start += 1;
}
for(i=start; i<=end; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
DLong_PushBack( value->xLong.value, DaoDecodeInteger( tokens[i]->string.mbs ) );
}
break;
#endif
case DAO_STRING :
n = tokens[start]->string.size - 1;
for(i=1; i<n; i++){
char c1 = str[i];
char c2 = str[i+1];
if( c1 < 'A' || c1 > 'P' ) continue;
DString_AppendChar( value->xString.data, (char)((c1-'A')*16 + (c2-'A')) );
i += 1;
}
if( str[0] == '\"' ) DString_ToWCS( value->xString.data );
break;
case DAO_ENUM :
value->xEnum.value = DaoDecodeInteger( str );
break;
case DAO_ARRAY :
#ifdef DAO_WITH_NUMARRAY
if( tokens[start]->name != DTOK_LSB ) return next;
k = DaoParser_FindPairToken( self, DTOK_LSB, DTOK_RSB, start, end );
if( k < 0 ) return next;
n = 1;
for(i=start+1; i<k; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
n *= strtol( tokens[i]->string.mbs, 0, 0 );
}
if( n < 0 ) return next;
if( it1 == NULL || it1->tid == 0 || it1->tid > DAO_COMPLEX ) return next;
array = & value->xArray;
dims = DArray_New(0);
for(i=start+1; i<k; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
j = strtol( tokens[i]->string.mbs, 0, 0 );
DArray_Append( dims, (size_t) j );
}
n = dims->size;
DaoArray_ResizeArray( array, dims->items.pInt, n );
DArray_PushFront( types, it1 );
DArray_Delete( dims );
n = 0;
for(i=k+1; i<=end; i++){
j = i + 1;
while( j <= end && tokens[j]->name != DTOK_COMMA ) j += 1;
DaoParser_Deserialize( self, i, j-1, & tmp, types, ns, proc, omap );
switch( it1->tid ){
case DAO_INTEGER : array->data.i[n] = tmp->xInteger.value; break;
case DAO_FLOAT : array->data.f[n] = tmp->xFloat.value; break;
case DAO_DOUBLE : array->data.d[n] = tmp->xDouble.value; break;
}
i = j;
n += 1;
}
DArray_PopFront( types );
#endif
break;
case DAO_LIST :
list = & value->xList;
DArray_PushFront( types, it1 );
n = 0;
for(i=start; i<=end; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
DArray_Append( & list->items, NULL );
k = DaoParser_Deserialize( self, i, end, list->items.items.pValue + n, types, ns, proc, omap );
i = k - 1;
n += 1;
}
DArray_PopFront( types );
break;
case DAO_MAP :
map = & value->xMap;
n = 0;
for(i=start; i<=end; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
DaoValue_Clear( & tmp );
DaoValue_Clear( & tmp2 );
DArray_PushFront( types, it1 );
i = DaoParser_Deserialize( self, i, end, &tmp, types, ns, proc, omap );
DArray_PopFront( types );
if( tokens[i]->name == DTOK_COMMA ) continue;
if( map->items->size == 0 ){
if( tokens[i]->name == DTOK_COLON ){
DMap_Delete( map->items );
map->items = DHash_New( D_VALUE, D_VALUE );
}
}
if( tokens[i]->name == DTOK_COLON || tokens[i]->name == DTOK_FIELD ) i += 1;
DArray_PushFront( types, it2 );
i = DaoParser_Deserialize( self, i, end, &tmp2, types, ns, proc, omap );
DArray_PopFront( types );
node = DMap_Insert( map->items, (void*) tmp, (void*) tmp2 );
i -= 1;
n += 1;
}
break;
case DAO_TUPLE :
tuple = & value->xTuple;
n = 0;
for(i=start; i<=end; i++){
if( tokens[i]->name == DTOK_COMMA ) continue;
it1 = NULL;
if( type->nested && type->nested->size > n ){
it1 = type->nested->items.pType[n];
if( it1 && it1->tid == DAO_PAR_NAMED ) it1 = & it1->aux->xType;
}
DArray_PushFront( types, it1 );
i = DaoParser_Deserialize( self, i, end, tuple->items + n, types, ns, proc, omap );
DArray_PopFront( types );
i -= 1;
n += 1;
}
break;
case DAO_OBJECT :
DArray_PushFront( types, NULL );
DaoParser_Deserialize( self, start, end, & tmp, types, ns, proc, omap );
DArray_PopFront( types );
if( tmp == NULL ) break;
object = DaoClass_MakeObject( & type->aux->xClass, tmp, proc );
if( object ) DaoValue_Copy( (DaoValue*) object, value2 );
break;
case DAO_CDATA :
case DAO_CSTRUCT :
DArray_PushFront( types, NULL );
DaoParser_Deserialize( self, start, end, & tmp, types, ns, proc, omap );
DArray_PopFront( types );
if( tmp == NULL ) break;
cdata = DaoCdata_MakeObject( & type->aux->xCdata, tmp, proc );
if( cdata ) DaoValue_Copy( (DaoValue*) cdata, value2 );
break;
}
DaoValue_Clear( & tmp );
DaoValue_Clear( & tmp2 );
return next;
}
