PLearn 0.1
|
00001 // -*- C++ -*- 00002 00003 // VMat_Func_utils.cc 00004 // 00005 // Copyright (C) 2004 Pascal Vincent 00006 // 00007 // Redistribution and use in source and binary forms, with or without 00008 // modification, are permitted provided that the following conditions are met: 00009 // 00010 // 1. Redistributions of source code must retain the above copyright 00011 // notice, this list of conditions and the following disclaimer. 00012 // 00013 // 2. Redistributions in binary form must reproduce the above copyright 00014 // notice, this list of conditions and the following disclaimer in the 00015 // documentation and/or other materials provided with the distribution. 00016 // 00017 // 3. The name of the authors may not be used to endorse or promote 00018 // products derived from this software without specific prior written 00019 // permission. 00020 // 00021 // THIS SOFTWARE IS PROVIDED BY THE AUTHORS ``AS IS'' AND ANY EXPRESS OR 00022 // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 00023 // OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN 00024 // NO EVENT SHALL THE AUTHORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, 00025 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED 00026 // TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 00027 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 00028 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 00029 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 00030 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00031 // 00032 // This file is part of the PLearn library. For more information on the PLearn 00033 // library, go to the PLearn Web site at www.plearn.org 00034 00035 /* ******************************************************* 00036 * $Id: VMat_Func_utils.cc 5557 2006-05-10 20:36:58Z lamblin $ 00037 ******************************************************* */ 00038 00039 // Authors: Pascal Vincent 00040 00044 #include "VMat_Func_utils.h" 00045 #include <plearn/var/Func.h> 00046 #include <plearn/vmat/VMat.h> 00047 00048 namespace PLearn { 00049 using namespace std; 00050 00051 // Put function implementations here. 00052 void evaluateSumOfFprop(VMat vm, Func f, Vec& output_result, int nsamples) 00053 { 00054 //if (f->outputs.size()!=1) 00055 // PLERROR("In evaluateSumOfFprop: function must have a single variable output (maybe you can concat the vars into a single one, if this is really what you want)"); 00056 00057 static int curpos = 0; 00058 int l = vm->length(); 00059 int w = vm->width(); 00060 if (nsamples == -1) nsamples = l; 00061 Vec input_value(w); 00062 Vec output_value(output_result.length()); 00063 00064 f->recomputeParents(); 00065 output_result.clear(); 00066 00067 for(int i=0; i<nsamples; i++) 00068 { 00069 vm->getRow(curpos++, input_value); 00070 f->fprop(input_value, output_value); 00071 output_result += output_value; 00072 if(curpos == l) curpos = 0; 00073 } 00074 } 00075 00076 void evaluateSumOfFbprop(VMat vm, Func f, Vec& output_result, Vec& output_gradient, int nsamples) 00077 { 00078 // if(f->outputs.size()!=1) 00079 // PLERROR("In evaluateSumOfFprop: function must have a single variable output (maybe you can concat the vars into a single one, if this is really what you want)"); 00080 00081 static int curpos = 0; 00082 int l = vm->length(); 00083 int w = vm->width(); 00084 if (nsamples == -1) nsamples = l; 00085 Vec input_value(w); 00086 Vec input_gradient(w); 00087 Vec output_value(output_result.length()); 00088 00089 f->recomputeParents(); 00090 output_result.clear(); 00091 00092 for(int i=0; i<nsamples; i++) 00093 { 00094 vm->getRow(curpos++, input_value); 00095 f->fbprop(input_value, output_value, input_gradient, output_gradient); 00096 //displayFunction(f, true); 00097 output_result += output_value; 00098 if(curpos == l) curpos = 0; 00099 } 00100 } 00101 00102 00103 } // end of namespace PLearn 00104 00105 00106 /* 00107 Local Variables: 00108 mode:c++ 00109 c-basic-offset:4 00110 c-file-style:"stroustrup" 00111 c-file-offsets:((innamespace . 0)(inline-open . 0)) 00112 indent-tabs-mode:nil 00113 fill-column:79 00114 End: 00115 */ 00116 // vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:encoding=utf-8:textwidth=79 :