Changes in sasmodels/kernelpy.py [e44432d:3199b17] in sasmodels


Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/kernelpy.py

    re44432d r3199b17  
    3333logger = logging.getLogger(__name__) 
    3434 
     35 
    3536class PyModel(KernelModel): 
    3637    """ 
     
    3839    """ 
    3940    def __init__(self, model_info): 
    40         # Make sure Iq is available and vectorized 
     41        # Make sure Iq is available and vectorized. 
    4142        _create_default_functions(model_info) 
    4243        self.info = model_info 
    4344        self.dtype = np.dtype('d') 
     45        logger.info("make python model " + self.info.name) 
    4446 
    4547    def make_kernel(self, q_vectors): 
     
    5254        """ 
    5355        pass 
     56 
    5457 
    5558class PyInput(object): 
     
    9093        self.q = None 
    9194 
     95 
    9296class PyKernel(Kernel): 
    9397    """ 
     
    130134        parameter_vector = np.empty(len(partable.call_parameters)-2, 'd') 
    131135 
    132         # Create views into the array to hold the arguments 
     136        # Create views into the array to hold the arguments. 
    133137        offset = 0 
    134138        kernel_args, volume_args = [], [] 
     
    173177                        else (lambda mode: 1.0)) 
    174178 
    175  
    176  
    177179    def _call_kernel(self, call_details, values, cutoff, magnetic, effective_radius_type): 
    178180        # type: (CallDetails, np.ndarray, np.ndarray, float, bool) -> np.ndarray 
     
    194196        self.q_input.release() 
    195197        self.q_input = None 
     198 
    196199 
    197200def _loops(parameters,    # type: np.ndarray 
     
    253256        total = np.zeros(nq, 'd') 
    254257        for loop_index in range(call_details.num_eval): 
    255             # update polydispersity parameter values 
     258            # Update polydispersity parameter values. 
    256259            if p0_index == p0_length: 
    257260                pd_index = (loop_index//pd_stride)%pd_length 
     
    264267            p0_index += 1 
    265268            if weight > cutoff: 
    266                 # Call the scattering function 
     269                # Call the scattering function. 
    267270                # Assume that NaNs are only generated if the parameters are bad; 
    268271                # exclude all q for that NaN.  Even better would be to have an 
     
    272275                    continue 
    273276 
    274                 # update value and norm 
     277                # Update value and norm. 
    275278                total += weight * Iq 
    276279                weight_norm += weight 
     
    292295    any functions that are not already marked as vectorized. 
    293296    """ 
    294     # Note: must call create_vector_Iq before create_vector_Iqxy 
     297    # Note: Must call create_vector_Iq before create_vector_Iqxy. 
    295298    _create_vector_Iq(model_info) 
    296299    _create_vector_Iqxy(model_info) 
Note: See TracChangeset for help on using the changeset viewer.