00001 // -*- C++ -*- 00002 00003 // PLearn (A C++ Machine Learning Library) 00004 // Copyright (C) 1998 Pascal Vincent 00005 // Copyright (C) 1999-2002 Pascal Vincent, Yoshua Bengio, Rejean Ducharme and University of Montreal 00006 // Copyright (C) 2001-2002 Nicolas Chapados, Ichiro Takeuchi, Jean-Sebastien Senecal 00007 // Copyright (C) 2002 Xiangdong Wang, Christian Dorion 00008 00009 // Redistribution and use in source and binary forms, with or without 00010 // modification, are permitted provided that the following conditions are met: 00011 // 00012 // 1. Redistributions of source code must retain the above copyright 00013 // notice, this list of conditions and the following disclaimer. 00014 // 00015 // 2. Redistributions in binary form must reproduce the above copyright 00016 // notice, this list of conditions and the following disclaimer in the 00017 // documentation and/or other materials provided with the distribution. 00018 // 00019 // 3. The name of the authors may not be used to endorse or promote 00020 // products derived from this software without specific prior written 00021 // permission. 00022 // 00023 // THIS SOFTWARE IS PROVIDED BY THE AUTHORS ``AS IS'' AND ANY EXPRESS OR 00024 // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 00025 // OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN 00026 // NO EVENT SHALL THE AUTHORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, 00027 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED 00028 // TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 00029 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 00030 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 00031 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 00032 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00033 // 00034 // This file is part of the PLearn library. For more information on the PLearn 00035 // library, go to the PLearn Web site at www.plearn.org 00036 00037 00038 /* ******************************************************* 00039 * $Id: KernelVMatrix.cc,v 1.2 2004/02/20 21:11:45 chrish42 Exp $ 00040 * This file is part of the PLearn library. 00041 ******************************************************* */ 00042 00043 #include "KernelVMatrix.h" 00044 00045 // From Old Kernel.cc: all includes are putted in every file. 00046 // To be revised manually 00047 #include <cmath> 00048 #include "stringutils.h" 00049 #include "Kernel.h" 00050 #include "TMat_maths.h" 00051 #include "PLMPI.h" 00053 namespace PLearn { 00054 using namespace std; 00055 00056 00057 // ***************** 00058 // * KernelVMatrix * 00059 // ***************** 00060 00061 KernelVMatrix::KernelVMatrix(VMat data1, VMat data2, Ker the_ker) 00062 : VMatrix(data1->length(), data2->length()), 00063 d1(data1), d2(data2), ker(the_ker), 00064 input1(data1->width()), input2(data2->width()) 00065 {} 00066 00067 00068 /* 00069 IMPLEMENT_NAME_AND_DEEPCOPY(KernelVMatrix); 00070 void KernelVMatrix::makeDeepCopyFromShallowCopy(map<const void*, void*>& copies) 00071 { 00072 Kernel::makeDeepCopyFromShallowCopy(copies); 00073 deepCopyField(d1, copies); 00074 deepCopyField(d2, copies); 00075 deepCopyField(ker, copies); 00076 deepCopyField(input1, copies); 00077 deepCopyField(input2, copies); 00078 } 00079 */ 00080 00081 real KernelVMatrix::get(int i, int j) const 00082 { 00083 #ifdef BOUNDCHECK 00084 if(i<0 || i>=length() || j<0 || j>=width()) 00085 PLERROR("In KernelVMatrix::get OUT OF BOUNDS"); 00086 #endif 00087 00088 d1->getRow(i,input1); 00089 d2->getRow(j,input2); 00090 return ker(input1,input2); 00091 } 00092 00093 00094 void KernelVMatrix::getSubRow(int i, int j, Vec v) const 00095 { 00096 #ifdef BOUNDCHECK 00097 if(i<0 || i>=length() || j<0 || j+v.length()>width()) 00098 PLERROR("In KernelVMatrix::getRow OUT OF BOUNDS"); 00099 #endif 00100 00101 d1->getRow(i,input1); 00102 for(int jj=0; jj<v.length(); jj++) 00103 { 00104 d2->getRow(j+jj,input2); 00105 v[jj] = ker(input1,input2); 00106 } 00107 } 00108 00109 00110 00111 } // end of namespace PLearn 00112