|
RTOpPack: Extra C/C++ Code for Vector Reduction/Transformation Operators
Version of the Day
|
00001 /* 00002 // @HEADER 00003 // *********************************************************************** 00004 // 00005 // Moocho: Multi-functional Object-Oriented arCHitecture for Optimization 00006 // Copyright (2003) Sandia Corporation 00007 // 00008 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00009 // license for use of this work by or on behalf of the U.S. Government. 00010 // 00011 // Redistribution and use in source and binary forms, with or without 00012 // modification, are permitted provided that the following conditions are 00013 // met: 00014 // 00015 // 1. Redistributions of source code must retain the above copyright 00016 // notice, this list of conditions and the following disclaimer. 00017 // 00018 // 2. Redistributions in binary form must reproduce the above copyright 00019 // notice, this list of conditions and the following disclaimer in the 00020 // documentation and/or other materials provided with the distribution. 00021 // 00022 // 3. Neither the name of the Corporation nor the names of the 00023 // contributors may be used to endorse or promote products derived from 00024 // this software without specific prior written permission. 00025 // 00026 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY 00027 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 00028 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 00029 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE 00030 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 00031 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 00032 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 00033 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 00034 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 00035 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 00036 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00037 // 00038 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 00039 // 00040 // *********************************************************************** 00041 // @HEADER 00042 */ 00043 00044 #include "RTOp_obj_values_vtbl.h" 00045 #include "RTOp_obj_free_free.h" 00046 00047 #include <stdlib.h> 00048 00049 /* Local function definitions */ 00050 00051 static int get_obj_type_num_entries( 00052 const struct RTOp_obj_type_vtbl_t* vtbl 00053 ,const void* instance_data 00054 ,int* num_values 00055 ,int* num_indexes 00056 ,int* num_chars 00057 ) 00058 { 00059 *num_values = *(RTOp_index_type*)instance_data; 00060 *num_indexes = 0; 00061 *num_chars = 0; 00062 return 0; 00063 } 00064 00065 static int obj_create( 00066 const struct RTOp_obj_type_vtbl_t* vtbl, const void* instance_data, void** obj ) 00067 { 00068 int k; 00069 const int 00070 num_values = *(RTOp_index_type*)instance_data, 00071 mem_size = sizeof(RTOp_value_type) * num_values; 00072 *obj = malloc( mem_size ); 00073 for( k = 0; k < num_values; ++k ) 00074 ((RTOp_value_type*)*obj)[k] = 0.0; 00075 return 0; 00076 } 00077 00078 static int obj_reinit( 00079 const struct RTOp_obj_type_vtbl_t* vtbl, const void* instance_data, void* obj ) 00080 { 00081 const int num_values = *(RTOp_index_type*)instance_data; 00082 int k; 00083 for( k = 0; k < num_values; ++k ) 00084 ((RTOp_value_type*)obj)[k] = 0.0; 00085 return 0; 00086 } 00087 00088 static int extract_state( 00089 const struct RTOp_obj_type_vtbl_t* vtbl 00090 ,const void * instance_data 00091 ,void * obj 00092 ,int num_values 00093 ,RTOp_value_type value_data[] 00094 ,int num_indexes 00095 ,RTOp_index_type index_data[] 00096 ,int num_chars 00097 ,RTOp_char_type char_data[] 00098 ) 00099 { 00100 int num_values_state; 00101 int k; 00102 assert(instance_data); 00103 num_values_state = *(RTOp_index_type*)instance_data; 00104 assert(obj); 00105 assert( num_values == num_values_state ); 00106 assert( num_indexes == 0 ); 00107 assert( num_chars == 0 ); 00108 for( k = 0; k < num_values; ++k ) 00109 value_data[k] = ((RTOp_value_type*)obj)[k]; 00110 return 0; 00111 } 00112 00113 static int load_state( 00114 const struct RTOp_obj_type_vtbl_t* vtbl 00115 ,const void * instance_data 00116 ,int num_values 00117 ,const RTOp_value_type value_data[] 00118 ,int num_indexes 00119 ,const RTOp_index_type index_data[] 00120 ,int num_chars 00121 ,const RTOp_char_type char_data[] 00122 ,void ** obj 00123 ) 00124 { 00125 int num_values_state; 00126 int k; 00127 assert(instance_data); 00128 num_values_state = *(RTOp_index_type*)instance_data; 00129 assert( obj ); 00130 assert( num_values == num_values_state ); 00131 assert( num_indexes == 0 ); 00132 assert( num_chars == 0 ); 00133 if(*obj == NULL) 00134 obj_create(vtbl,instance_data,obj); 00135 for( k = 0; k < num_values; ++k ) 00136 ((RTOp_value_type*)*obj)[k] = value_data[k]; 00137 return 0; 00138 } 00139 00140 const struct RTOp_obj_type_vtbl_t RTOp_obj_values_vtbl = 00141 { 00142 get_obj_type_num_entries 00143 ,obj_create 00144 ,obj_reinit 00145 ,RTOp_obj_free_free 00146 ,extract_state 00147 ,load_state 00148 };
1.7.6.1