- (NSDate *)dateFromString:(NSString *)string {
    //Wed Mar 14 16:40:08 +0800 2012
    if (!string) return nil;
    NSDateFormatter *dateformatter=dateformatter = [[NSDateFormatter alloc] init];
    NSTimeZone *tz = [NSTimeZone timeZoneWithAbbreviation:@"GMT"];
    [dateformatter setTimeZone:tz];
    [dateformatter setDateFormat:@"EEE MMM dd HH:mm:ss Z yyyy"];;
    return [dateformatter dateFromString:string];
}

- (NSDate*)mydateFromString:(NSString *)string;
{
    //Wed Mar 14 16:40:08 +0800 2012
    if (!string) return nil;
    static NSCalendar *gregorian=nil;
    static NSDateComponents *comps=nil;
    static NSDictionary *month=nil;
    if (gregorian==nil) {
        gregorian = [[NSCalendar alloc] initWithCalendarIdentifier:NSGregorianCalendar];
        comps = [[NSDateComponents alloc] init];
        month = [NSDictionary dictionaryWithObjectsAndKeys:
                 [NSNumber numberWithInt:1], @"Jan",
                 [NSNumber numberWithInt:2], @"Feb",
                 [NSNumber numberWithInt:3], @"Mar",
                 [NSNumber numberWithInt:4], @"Apr",
                 [NSNumber numberWithInt:5], @"May",
                 [NSNumber numberWithInt:6], @"Jun",
                 [NSNumber numberWithInt:7], @"Jul",
                 [NSNumber numberWithInt:8], @"Aug",
                 [NSNumber numberWithInt:9], @"Sep",
                 [NSNumber numberWithInt:10], @"Oct",
                 [NSNumber numberWithInt:11], @"Nov",
                 [NSNumber numberWithInt:12], @"Dec",
                 nil];
    }

    @try {
        NSString *t=[string substringWithRange:NSMakeRange(26, 4)];
        [comps setYear:[t intValue]];
        t=[string substringWithRange:NSMakeRange(4, 3)];
        [comps setMonth:[[month objectForKey:t] intValue]];
        t=[string substringWithRange:NSMakeRange(8, 2)];
        [comps setDay:[t intValue]];
        t=[string substringWithRange:NSMakeRange(11, 2)];
        [comps setHour:[t intValue]];
        t=[string substringWithRange:NSMakeRange(14, 2)];
        [comps setMinute:[t intValue]];
        t=[string substringWithRange:NSMakeRange(17, 2)];
        [comps setSecond:[t intValue]];
        t=[string substringWithRange:NSMakeRange(20, 5)];
        //全球共24个标准时区，每个时区为1个小时，下面计算该时区offset秒数
        [comps setTimeZone:[NSTimeZone timeZoneForSecondsFromGMT:[t intValue]/100*3600]];
    }
    @catch (NSException *exception) {
    }
    @finally {
    }

    return [gregorian dateFromComponents:comps];
}

- (NSDate *)dateFromString:(NSString *)string {
    //Wed Mar 14 16:40:08 +0800 2012
    if (!string) return nil;
    struct tm tm;
    time_t t;
    string=[string substringFromIndex:4];
    strptime([string cStringUsingEncoding:NSUTF8StringEncoding], "%b %d %H:%M:%S %z %Y", &tm);
    tm.tm_isdst = -1;
    t = mktime(&tm);
    return [NSDate dateWithTimeIntervalSince1970:t];
}

2012-05-05 15:57:51.942 timetest[18488:707] 1.991521 //第一种
2012-05-05 15:57:52.096 timetest[18488:707] 0.921144　//第二种
2012-05-05 15:57:54.088 timetest[18488:707] 0.153897　//第三种

function levenshtein(a, b) {
    var i;
    var j;
    var cost;
    var d = new Array();

    if (a.length == 0) {
        return b.length;
    }

    if (b.length == 0) {
        return a.length;
    }

    for (i = 0; i <= a.length; i++) {
        d[i] = new Array();
        d[i][0] = i;
    }

    for (j = 0; j <= b.length; j++) {
        d[0][j] = j;
    }

    for (i = 1; i <= a.length; i++) {
        for (j = 1; j <= b.length; j++) {

          cost = (a.charAt(i - 1) == b.charAt(j - 1)) ? 0 : 1;
          
            d[i][j] = Math.min(d[i - 1][j] + 1, d[i][j - 1] + 1, d[i - 1][j - 1] + cost);

            if (i > 1 && j > 1 &&
              a.charAt(i - 1) == b.charAt(j - 2) &&
              a.charAt(i - 2) == b.charAt(j - 1)) {
                d[i][j] = Math.min(d[i][j], d[i - 2][j - 2] + cost)
            }
        }
    }
    return d[a.length][b.length];
}

function checkCommonMail(id,errorid) {
    var email = document.getElementById(id).value;
    var useremailDomain=email.substring(email.indexOf("@")+1).toLowerCase();
      var error=document.getElementById(errorid);
      
      var commonDomain = [
          "hotmail.com",
          "gmail.com",
          ["yahoo.cn","yahoo.hk","yahoo.com.cn","yahoo.com.hk","yahoo.com"],
          "163.com",
          "sina.com",
          "live.com",
          "qq.com",
          "vip.qq.com"];

  for(var i=0;i<commonDomain.length;i++){
      if (commonDomain[i] instanceof Array) {
          var pardomain=commonDomain[i];
          var dis=new Array();
          for(var j=0;j<pardomain.length;j++){
              dis[j]=levenshtein(useremailDomain,pardomain[j]);
          }
          var mindis=100000,minj=0;
          for(var j=0;j<dis.length;j++){
              if(dis[j]==0) return;
              if(dis[j]<mindis) {
                  mindis=dis[j];
                  minj=j;
              }
          }
          if(mindis==0) return;
          if(mindis<3){
              error.innerHTML = '您的邮件地址很像'+pardomain[minj]+'？';
          }
          delete dis;
          return;
      }else{
          var dis=levenshtein(useremailDomain,commonDomain[i]);
          if(dis==0) return;
          if(dis<2){
              error.innerHTML = '您的邮件地址很像'+commonDomain[i]+'？';
              return;
          }
          delete dis;
      }
  }
}

function checkMail(id) {
  var email = document.getElementById(id).value;
    var pattern = /^\w+([-+.]\w+)*@\w+([-.]\w+)*\.\w+([-.]\w+)*$/;
    if (pattern.test(email)) {
        return true;
    } else {
        return false;
    }
}

$('#content a.button')
  .bind('click.button', function (e) {
    $(this).data('clicked', true).unbind()
    e.preventDefault()
  })
  .css({
      opacity: 1
    , color: 'red'
  })
  .fadeOut(250)
$.map(['a', 'b', 'c'], function (letter) {
  return letter.toUpperCase()
})
$.ajax('/data', function (resp) {
  $('#content').html(resp)
})

/*!
  * =============================================================
  * Ender: open module JavaScript framework (https://ender.no.de)
  * Build: ender build jeesh reqwest
  * =============================================================
  */

$.domReady(function () {...})

$('#boosh a[rel~="bookmark"]').each(function (el) { ... })

$('#boosh p a[rel~="bookmark"]').hide().html('hello').css({
  color: 'red',
  'text-decoration': 'none'
}).addClass('blamo').after('✓').show();

$('#content a').bind('keydown input', handler)
$('#content a').emit('customEvent')
$('#content a').remove('click.myClick')

reqwest({
    url: 'path/to/html'
  , method: 'post'
  , data: { foo: 'bar', baz: 100 }
  , success: function (resp) {
      qwery('#content').html(resp)
    }
})

$.ajax({
    url: 'path/to/html'
  , method: 'post'
  , data: { foo: 'bar', baz: 100 }
  , success: function (resp) {
      qwery('#content').html(resp)
    }
})

{% instagram media_id %}

rake new_page['photos']

rake new_page['photos/instagram.php']

---
layout: page
title: "我的照片"
date: 2011-12-21 05:22
comments: true
sharing: true
footer: true
description: My Photos
---

<?php
include("./instagram.php");
?>

<?php

  $access_token='19120281.12eb222.1614fb6d1334221222d66fcd5f258fb';
  $count_images=90;
  $cache_time=60;
  $image_size='thumbnail';
  
  /**
  Output each image with HTML markup
 */
  function echoimage($val, $imagesize) {
      $image = $val["images"][$imagesize]["url"];
      $link = $val["link"];
      $tag=md5($link);
      return "<a href=\"$link\" target=\"_blank\"><img src=\"$image\"/></a>";
      }
      
  /**
  Getting the Data from the API
 */
  function getDataFromApi($token, $number){
      // Instagram API bearbeiten
      $url="https://api.instagram.com/v1/users/self/media/recent/?access_token=$token&count=$number";
      $contents = file_get_contents($url);
      return $contents;
  }
  
  /**
  Gets the Data from either the Cache or the API
 */
  function getData($token, $number, $cache_time){

      $cache_folder = "your cache path";
  
      if(!is_dir($cache_folder)){
        if(mkdir($cache_folder, 0777) == false){
          return getDataFromApi($token, $number);
        }
      }
  
      $cachefile = $cache_folder."user.cache";
  
      if (file_exists($cachefile) && time()-filemtime($cachefile)<$cache_time) {
        try{
          $contents = file_get_contents($cachefile);
        }catch(Exception $e){
          $contents = getDataFromApi($token, $number);
          file_put_contents($cachefile, $contents);
        }
  
      } else {
        $contents = getDataFromApi($token, $number);
        file_put_contents($cachefile, $contents);
      }
  
      return $contents;
  }
  
  $json = json_decode(getData($access_token, $count_images, $cache_time), true);

  echo '<div class="photosdiv">';
  foreach ($json["data"] as $value)
      echo echoimage($value, $image_size);
  echo '</div>';
  
?>

id obj=self;
if ([obj respondsToSelector:@selector(function1:)) {
}
if ([obj isKindOfClass:[NSArray class]] ) {
}
if ([obj conformsToProtocol:@protocol(myProtocol)]) {
}
if ([[obj class] isSubclassOfClass:[NSArray class]]) {
}
[obj someNonExistFunction];

id knife=[Knife grateKnife];
NSArray *monsterList=[NSArray array];
[monsterList makeObjectsPerformSelector:@selector(killMonster:) withObject:knife];

/********************************************************************
 * id        objc_msgSend(id self,
 *           SEL op,
 *           ...)
 *
 * On entry: a1 is the message receiver,
 *           a2 is the selector
 ********************************************************************/

  ENTRY objc_msgSend
# check whether receiver is nil
  teq     a1, #0
  moveq   a2, #0
  bxeq    lr
 
# save registers and load receiver's class for CacheLookup
  stmfd   sp!, {a4,v1-v3}
  ldr     v1, [a1, #ISA]

# receiver is non-nil: search the cache
  CacheLookup a2, LMsgSendCacheMiss

# cache hit (imp in ip) - prep for forwarding, restore registers and call
  teq  v1, v1     /* set nonstret (eq) */
  ldmfd   sp!, {a4,v1-v3}
  bx      ip

# cache miss: go search the method lists
LMsgSendCacheMiss:
  ldmfd    sp!, {a4,v1-v3}
  b    _objc_msgSend_uncached

LMsgSendExit:
  END_ENTRY objc_msgSend


  .text
  .align 2
_objc_msgSend_uncached:

# Push stack frame
  stmfd    sp!, {a1-a4,r7,lr}
  add     r7, sp, #16
  SAVE_VFP

# Load class and selector
  ldr  a1, [a1, #ISA]        /* class = receiver->isa  */
 # MOVE  a2, a2          /* selector already in a2 */

# Do the lookup
  MI_CALL_EXTERNAL(__class_lookupMethodAndLoadCache)
  MOVE    ip, a1

# Prep for forwarding, Pop stack frame and call imp
  teq  v1, v1     /* set nonstret (eq) */
  RESTORE_VFP
  ldmfd    sp!, {a1-a4,r7,lr}
  bx   ip

- (void)myThread:(id)sender
{
  NSAutoreleasePool *pool=[[NSAutoreleasePool alloc] init];
  while (TRUE) {
      //do some jobs
      //break in some condition
      usleep(10000);
      [pool drain];
  }
  [pool release];
}

NSMutableArray *targetQueue;
NSMutableArray *actionQueue;
- (void)myThread:(id)sender
{
    NSAutoreleasePool *pool=[[NSAutoreleasePool alloc] init];
    while (TRUE) {

        //do some jobs
        //break in some condition
        int n=[targetQueue count];
        assert(n==[actionQueue count]);
        for(int i=0;i<n;i++){
            id target=[targetQueue objectAtIndex:i];
            SEL action=NSSelectorFromString([actionQueue objectAtIndex:i]);
            if ([target respondsToSelector:action]) {
                [target performSelector:action withObject:nil];
            }
        }

        usleep(10000);

        [pool drain];
    }

    [pool release];
}

while (TRUE) {
  //break in some condition
}

@interface MyNSTimer : NSObject
{
  id target;
  SEL action;
  float interval;
  CFAbsoluteTime lasttime;
}
- (void)invoke;
@end

@implementation MyNSTimer
- (void)invoke;
{
  if ([target respondsToSelector:action]) {
        [target performSelector:action withObject:nil];
    }
}
@end

@interface MyNSRunloop : NSObject
{
  NSMutableArray *timerQueue;
}
- (void)addTimer:(MyNSTimer*)t;
- (void)executeOnce;
@end

@implementation MyNSRunloop
- (void)addTimer:(MyNSTimer*)t;
{
  @synchronized(timerQueue){
      [timerQueue addObject:t];
    }
}
- (void)executeOnce;
{
  CFAbsoluteTime currentTime=CFAbsoluteTimeGetCurrent();
  @synchronized(timerQueue){
      for(MyNSTimer *t in timerQueue){
          if(currentTime-t.lasttime>t.interval){
              t.lasttime=currentTime;
                  [t invoke];
          }
      }
  }
}
@end

@interface MyNSThread : NSObject
{
  MyNSRunloop *runloop;
}
- (void)main:(id)sender;
@end

@implementation MyNSThread
- (void)main:(id)sender
{
    NSAutoreleasePool *pool=[[NSAutoreleasePool alloc] init];
    while (TRUE) {
        //do some jobs
        //break in some condition
        [runloop executeOnce];
        usleep(10000);
        [pool drain];
    }
    [pool release];
}
@end

MyNSTimer *timer1=[MyNSTimer scheduledTimerWithTimeInterval:1
                                                   target:obj1
                                                 selector:@selector(download1:)];
[myThread1.runloop addTimer:timer1];

MyNSTimer *timer2=[MyNSTimer scheduledTimerWithTimeInterval:2
                                                   target:obj2
                                                 selector:@selector(download2:)];
[myThread1.runloop addTimer:timer2];

//
//  LXTextField.h
//  lxtextfield
//
//  Created by xu lian on 12-03-09.
//  Copyright (c) 2012 Beyondcow. All rights reserved.
//

#import <AppKit/AppKit.h>

@interface LXTextField : NSView
{
    NSString *stringValue;
    BOOL selectAll;
}
@property(retain,nonatomic) NSString *stringValue;

@end

//
//  LXTextField.m
//  lxtextfield
//
//  Created by xu lian on 12-03-09.
//  Copyright (c) 2012 Beyondcow. All rights reserved.
//

#import "LXTextField.h"

@implementation LXTextField
@synthesize stringValue;

- (void)awakeFromNib
{
    selectAll = NO;
}

- (id)initWithFrame:(NSRect)frameRect
{
    if( self = [super initWithFrame:frameRect] ){
        selectAll = NO;
    }
    return self;
}

- (BOOL)acceptsFirstResponder
{
    return YES;
}

- (BOOL)becomeFirstResponder
{
    return YES;
}

- (BOOL)resignFirstResponder
{
    selectAll=NO;
    [self setNeedsDisplay:YES];
    return YES;
}


- (void)setStringValue:(NSString *)string{
    stringValue = string;
    [self setNeedsDisplay:YES];
}

- (void)drawRect:(NSRect)dirtyRect
{
    if (selectAll) {
        NSRect r = NSZeroRect;
        r.size = [stringValue sizeWithAttributes:nil];
        [[NSColor selectedControlColor] set];
        NSRectFill(r);
    }
    [stringValue drawAtPoint:NSZeroPoint withAttributes:nil];
}

- (IBAction)selectAll:(id)sender;
{
    selectAll=YES;
    [self setNeedsDisplay:YES];
}

- (IBAction)copy:(id)sender;
{
    NSPasteboard *pasteBoard = [NSPasteboard generalPasteboard];
    [pasteBoard declareTypes:[NSArray arrayWithObjects:NSStringPboardType, nil] owner:nil];
    [pasteBoard setString:stringValue forType:NSStringPboardType];
}

- (void)mouseDown:(NSEvent *)theEvent
{
    if ([theEvent clickCount]>=2) {
        selectAll=YES;
    }
    [self setNeedsDisplay:YES];
}

- (void)keyDown:(NSEvent *)theEvent
{
}

@end

- (void)keyDown:(NSEvent *)theEvent
{
  if(condition){
      do something;
  }else{
      [super keyDown:theEvent];
  }
}

NSLog(@"%d",[object retainCount]);

while([object retainCount]>0){
      [object release];
}

/**************************************************
*   VSM by keefo 2008.05.26
**************************************************/
     
#include <stdio.h>
#include <math.h>
#include <map>
#include <string>
#include <vector>
using namespace std;
#define FILENUM                50 
#define isDigit(ch) ((ch)>='0' && (ch)<='9')
#define isAlpha(ch) (((ch)>='a' && (ch)<='z') || ((ch)>='A' && (ch)<='Z'))
#define isSign(ch) (((ch)<127 && (ch)>=0)?signindex[ch]:0)
map<string, int> dic; //全局词典。文本集里出现过的词全部保存在这里
map<string, int> filedic[FILENUM];
map<string, int> stopword;//停用词表
vector<string> diclist;//单篇文档词典
int signindex[200]={0};//符号hash表用来去除文章里的符号


int loadstopword() //加载停用词到stopword中
{
    FILE *fp;
    char str[32];
    stopword.clear();
    fp=fopen("stopwords.txt","r");
    if(fp==NULL)
    {
        printf("Can't open stopword file.\n");
        system("pause");
        return 1;
    }
    printf("Loading stopword...");
    while(fscanf(fp,"%s",str)!=EOF)
        stopword[str]=1;
    fclose(fp);
    printf("OK!\n");
    return 0;
}

int trimfile()
{
    int i;
    char ch,name[128],word[128];
    FILE *fp1,*fp2;
    printf("Create trim file waiting...\n");
    //去除符号
    for(i=0;i<FILENUM;i++)
    {
        sprintf(name,"file/file%d.txt",i+1);//训练数据保存在file目录下
        fp1=fopen(name,"r");
        if(fp1==NULL){
            printf("Open file %s error!\n",name);
            system("pause");
            return 1;
        }
        sprintf(name,"trim%d.txt",i+1);
        fp2=fopen(name,"w+");
        if(fp2==NULL){
            printf("Create file %s error!\n",name);
            system("pause");
            return 1;
        }
        ch=fgetc(fp1);
        while(ch!=EOF)
        {
            if(!(isDigit(ch) || isAlpha(ch) || isSign(ch)))
            {
                fputc(ch,fp2);
            }
            ch=fgetc(fp1);
        }
        fclose(fp1);
        fclose(fp2);
    }
    //去除停用词
    for(i=0;i<FILENUM;i++)
    {
        sprintf(name,"trim%d.txt",i+1);
        fp1=fopen(name,"r+");
        if(fp1==NULL){
            printf("Can't open file %s!\n",name);
            system("pause");
            return 1;
        }
        sprintf(name,"file%d.txt",i+1);
        printf("%d ",i+1);
        fp2=fopen(name,"w");
        if(fp2==NULL)
        {
            printf("Create file %s error!\n",name);
            system("pause");
            return 1;
        }
        while(fscanf(fp1,"%s",word)!=EOF)
        {
            if(!stopword[word])
                fprintf(fp2,"%s ",word);
        }
        fclose(fp1);
        fclose(fp2);
    }
    printf("\n");
    system("del trim*.txt");//删除临时文件
    return 0;
}


void initsign()//初始化符号表
{
    memset(signindex,0,sizeof(signindex));
    signindex[33]=1;
    signindex[34]=1;
    signindex[37]=1;
    signindex[38]=1;
    signindex[40]=1;
    signindex[41]=1;
    signindex[42]=1;
    signindex[43]=1;
    signindex[44]=1;
    signindex[45]=1;
    signindex[46]=1;
    signindex[47]=1;
    signindex[58]=1;
    signindex[59]=1;
    signindex[60]=1;
    signindex[61]=1;
    signindex[62]=1;
    signindex[63]=1;
    signindex[91]=1;
    signindex[93]=1;
    signindex[95]=1;
    signindex[123]=1;
    signindex[125]=1;
    signindex[126]=1;
}

int allwordnum,process[11]={0};
void processing(int k)//显示计算进度
{
    float s=(float)k/allwordnum*100;
    int ca=(int)s/10;
    if(!process[ca])
    {
        if(ca==10)
            printf("%d%% OK!",ca*10);
        else
            printf("%d%% ",ca*10);
        process[ca]=1;
    }
}


int main(int argc,char *argv[])
{
    FILE *fp;
    char name[128],word[128];
    int i,j,filewordnum[FILENUM];
    initsign();
    if(loadstopword())
        return 1;
    if(trimfile())
        return 1;
    //计算tf-idf建立向量空间模型
    printf("Create VSM waiting ...\n");
    diclist.clear();
    dic.clear();
    for(i=0;i<FILENUM;i++) //按顺序打开处理好的文本，统计词
    {
        filedic[i].clear();
        sprintf(name,"file%d.txt",i+1);
        printf("%d ",i+1);
        fp=fopen(name,"r");
        if(fp == NULL){
            printf("Processing file %d error!",i+1);
            system("pause");
            continue;
        }
        filewordnum[i]=0;
        while(fscanf(fp,"%s",word) != EOF) //循环查找该文件内所有单词
        {
            if(filedic[i][word])
            {
                //如果在filedic[i]里发现这个单词,出现个数增加，总字典内个数增加
                filedic[i][word]++;
                dic[word]++;
            }
            else
            {
                filedic[i][word]=1;//没出现过的新单词加入filedic[i]
                if(dic[word]==0)
                {
                    //如果也不在dic里，同样加入dic和diclist内。
                    diclist.push_back(word);
                    dic[word]=1;
                }
                else
                {
                    //在总字典里出现过，出现个数增加
                    dic[word]++;
                }
            }
            filewordnum[i]++;
        }
        fclose(fp);
    }
    system("del file*.txt");//删除临时文件
    fp=fopen("T_v.txt","w");
    if(fp == NULL){
        printf("Create v file error.");
        system("pause");
        return 1;
    }
    allwordnum=diclist.size();

    printf("\n计算&打印VSM T_v.txt\n");
    printf("一共有 %d 个词需要计算 waiting ... \n",allwordnum);
    fprintf(fp,"%d %d\n",FILENUM,allwordnum);
    int knum;
    string key;
    double tf,idf;
    FILE *fpword;
    fpword=fopen("T_word.txt","w");
    fprintf(fpword,"%d\n",allwordnum);
    for(i=0;i<allwordnum;i++)
    {
        key=diclist[i];
        fprintf(fpword,"%s\n",key.c_str());
        for(knum=j=0;j<FILENUM;j++)
            if(filedic[j][key])
                knum++;//knum出现该词的文档数

        idf=log((double)FILENUM/knum);//逆文频率指数,这里注意要排除全局都出现的关键词。log(1)=0 会导致后期迭代失败
        for(j=0;j<FILENUM;j++)
        {
            tf=(double)filedic[j][key]/filewordnum[j];//词频
            fprintf(fp,"%.8lf ", tf*idf);
        }
        fprintf(fp,"\n");
        processing(i+1);//显示计算进度
    }
    fclose(fp);
    fclose(fpword);
    return 0;
}

/**************************************************
*   NMF by keefo 2008.05.26
**************************************************/
 
#include <stdio.h>
#include <string.h>
#include <memory.h>
#include <stdlib.h>
#include <time.h>

#define ARRAY2D(VTYPE,X,Y) (VTYPE **) malloc2d(X, Y, sizeof(VTYPE))
#define MATRIX(NAME,X,Y) matrix NAME;NAME.p=ARRAY2D(double,X,Y);NAME.x=X;NAME.y=Y;
#define FREE(X) free2d(&X)

typedef struct{
  double **p;
  int x,y;
} matrix;

void** malloc2d(int W, int H, int size)
{
  int j,rowSize=W*size,colSize=H*sizeof(void *);
  void **a = (void **)malloc(colSize+H*rowSize);
  char *begin=(char *)a+colSize;
  for(j=0;j<H;j++)
        a[j]=begin+j*rowSize;
  return a;
}

void free2d(matrix *M)
{
  free(M->p);
  M->p=0;
  M->x=M->y=0;
}

matrix RandomMatrix(int x,int y) //生成随机矩阵
{
  int i,j;
  MATRIX(A,x,y);
  srand((int)time(NULL));
  for(i=0;i<y;i++)
  {
        for(j=0;j<x;j++)
        {
            A.p[i][j]=((double)rand()/((double)RAND_MAX+1))*10;
        }
  }
  return A;
}

matrix ONES(int x,int y) //生成全为1的矩阵
{
  int i,j;
  MATRIX(A,x,y);
  for(i=0;i<A.y;i++)
  {
        for(j=0;j<A.x;j++)
            A.p[i][j]=1.0;
  }
  return A;
}

void SUM(matrix A,matrix *B)
{
  int i,j;
  double sum;
  for(i=0;i<A.x;i++)
  {
        sum=0.0;
        for(j=0;j<A.y;j++)
            sum+=A.p[j][i];
        B->p[0][i]=sum;
  }
}

void DOTMUL(matrix A,matrix B,matrix *C) //矩阵点乘
{
  int i,j;
  if(A.x!=B.x || A.y!=B.y)
        return;
  
  for(i=0;i<A.y;i++)
        for(j=0;j<A.x;j++)
            C->p[i][j]=A.p[i][j]*B.p[i][j];
}

void MUL(matrix A,matrix B,matrix *C) //矩阵相乘
{
  int i,j,k;
  double sum;
  if(A.x!=B.y)
        return;
  for(i=0;i<A.y;i++)
  {
        for(j=0;j<B.x;j++)
        {
            sum=0.0;
            for(k=0;k<A.x;k++)
                sum+=A.p[i][k]*B.p[k][j];
            C->p[i][j]=sum;
        }
  }
}

void DOTDIV(matrix A,matrix B,matrix *C) //矩阵点除
{
  int i,j;
  if(A.x!=B.x || A.y!=B.y)
        return;
  
  for(i=0;i<A.y;i++)
        for(j=0;j<A.x;j++)
            C->p[i][j]=A.p[i][j]/B.p[i][j];
}

void T(matrix S,matrix *D) //矩阵转置
{
  int i,j;
  for(i=0;i<S.y;i++)
  {
        for(j=0;j<S.x;j++)
        {
            D->p[j][i]=S.p[i][j];
        }
  }
}

void fprint(char name[128],matrix m)
{
  int i,j;
  FILE *fp;
  fp=fopen(name,"w");
  fprintf(fp,"%d %d\n",m.x,m.y);
  if(m.p){
        for(i=0;i<m.y;i++){
            for(j=0;j<m.x;j++)
                fprintf(fp,"%.8f ",m.p[i][j]);
            fprintf(fp,"\n");
        }
  }
  fclose(fp);
}

int main(int argc, char * argv[])
{
  int i,j,x,y,r,iteratenum;
  char word[128];
  FILE *fp;
  
  fp=fopen("config.ini","r");
  if(fp==NULL){
      printf("Can't open file config.ini!\n");
      system("pause");
      return 1;
  }
  fscanf(fp,"%s %d",word,&r);
  fscanf(fp,"%s %d",word,&iteratenum);
  fclose(fp);
  
  fp=fopen("T_v.txt","r");
  if(fp==NULL){
        printf("Can't open file T_v.txt!\n");
        system("pause");
        return 1;
  }
  
  fscanf(fp,"%d %d",&x,&y);
  
  MATRIX(X,x,y);
  
  for(i=0;i<y;i++)
        for(j=0;j<x;j++)
            fscanf(fp,"%lf",&X.p[i][j]);
  fclose(fp);
  
  printf("r=%d\n",r);
  
  matrix B=RandomMatrix(r,X.y); //初始化B，为随机非负数
  matrix H=RandomMatrix(X.x,r); //初始化H，为随机非负数
  /*************定义中间变量*************/
  matrix one=ONES(1,X.y);
  MATRIX(XX,X.x,X.y);
  MATRIX(TB,B.y,B.x);
  MATRIX(BB,B.x,B.y);
  MATRIX(TH,H.y,H.x);
  MATRIX(HH,H.x,H.y);
  MATRIX(BH,H.x,B.y);
  MATRIX(SB,B.x,1);
  
  /*************归一化B的每一列*************/
  SUM(B,&SB);
  MUL(one,SB,&BB);
  DOTDIV(B,BB,&B);
  
  /*************开始迭代计算*************/
  printf("共%d次迭代计算\n",iteratenum);
  for(i=0;i<iteratenum;i++)
  {
      printf("%d ",i+1);
      // a：H = H .* ( B' * ( X ./ ( B * H ) ) )
      MUL(B,H,&BH);
      DOTDIV(X,BH,&XX);
      T(B,&TB);
      MUL(TB,XX,&HH);
      DOTMUL(H,HH,&H);
      // b：B = B .* ( ( X ./ ( B * H ) ) * H' )
      T(H,&TH);
      MUL(B,H,&BH);
      DOTDIV(X,BH,&XX);
      MUL(XX,TH,&BB);
      DOTMUL(B,BB,&B);
      // c：B 归一化
      SUM(B,&SB);
      MUL(one,SB,&BB);
      DOTDIV(B,BB,&B);
  }
  
  //MUL(B,H,&XX);//测试计算出来的BH相乘是否接近X
  //fprint("BxH.txt",XX);
  
  T(B,&TB);
  MUL(TB,X,&HH);
  
  printf("\n计算完成!\n");
  printf("写入数据 T_eigenvector.txt...\n");
  fprint("T_eigenvector.txt",HH);
  fprint("T_B.txt",TB);
  
  FREE(B);
  FREE(H);
  FREE(X);
  FREE(one);
  FREE(BB);
  FREE(HH);
  FREE(TB);
  FREE(TH);
  FREE(SB);
  FREE(XX);
  return 0;
}